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This work is designed to overview our present knowledge about universality classes occurring in 
nonequilibrium systems defined on regular lattices. In the first section I summarize the most 
important critical exponents, relations and the field theoretical formalism used in the text. In the 
second section I briefly address the question of scaling behavior at first order phase transitions. 
In section three I review dynamical extensions of basic static classes, show the efi'ect of mixing 
dynamics and the percolation behavior. The main body of this work is given in section four where 
genuine, dynamical universality classes specific to nonequilibrium systems are introduced. In 
section five I continue overviewing such nonequilibrium classes but in coupled, multi-component 
systems. Most of the known nonequilibrium transition classes are explored in low dimensions 
between active and absorbing states of reaction-diffusion type of systems. However by mapping 
they can be related to universal behavior of interface growth models, which I overview in section 
six. Finally in section seven I summarize families of absorbing state system classes, mean-field 
classes and give an outlook for further directions of research. 
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I. INTRODUCTION 

Universal scaling behavior is an attractive feature in 
statistical physics because a wide range of models can 
be classified purely in terms of their collective behav- 
ior. Scaling phenomenon has been observed in many 
branches of physics, chemistry, biology, economy ... etc., 
most frequently by critical phase transitions. Nonequi- 
librium phase transitions may appear in models of pop- 
ulation (Albano, 1994), epidemics (Liggct, 1985; MoUi- 
son, 1977), catalysis (Ziff et al, 1986), cooperative trans- 
port (Chowdhury et al., 2000; Havlin and ben Avra- 
ham, 1987), enzyme biology (Berry, 2003) and markets 
(Bouchaud and Georges, 1990) for example. 

Dynamical extensions of static universality classes — 
established in equilibrium — are the simplest nonequilib- 
rium models systems, but beyond that critical phenom- 
ena, with new classes have been explored so far (Grass- 
berger, 1996; Hinrichsen, 2000a; Marro and Dickman, 
1999). While the theory of phase transitions is quite 
well understood in thermodynamic equilibrium its re- 
search in nonequilibrium is rather new. In general phase 
transitions, scaling and universality retain most of the 
fundamental concepts of equilibrium models. The ba- 
sic ingredients affecting universality classes arc again the 
collective behavior of systems, the symmetries, the con- 
servation laws and the spatial dimensions as described by 
renormalization group theory. Besides them several new 
factors have also been identified recently. Low dimen- 
sional systems are of primary interest because the fluc- 
tuation effects are relevant, hence the mean-field type of 
description is not valid. In the past decades this field of 
research was growing very quickly and now we are faced 
with a zoo of models, universality classes, strange nota- 
tions and abbreviations. This article aims to help new- 
comers as well as researchers to navigate in the literature 



by systematically reviewing most of the explored univer- 
sality classes. I define models by their field theory (when 
it is available), show their symmetries or other impor- 
tant features and list the critical exponents and scaling 
relations. 

Nonequilibrium systems can be classified into two cat- 
egories: (a) Systems which do have a hermitian Hamilto- 
nian and whose stationary states are given by the proper 
Gibbs-Boltzmann distribution. However, they are pre- 
pared in an initial condition which is far from the station- 
ary state and sometimes, in the thermodynamic limit, 
the system may never reach the true equilibrium. These 
nonequilibrium systems include, for example, phase or- 
dering systems, spin glasses, glasses etc. I begin the re- 
view of classes by showing the scaling behavior of the 
simplest prototypes of such systems in section III These 
are defined by adding simple dynamics to static models, 
(b) Systems without a hermitian Hamiltonian defined by 
transition rates, which do not satisfy the detailed balance 
condition (the local time reversal symmetry is broken). 
They may or may not have a steady state and even if they 
have one, it is not a Gibbs state. Such models can be 
created by combining different dynamics or by generat- 
ing currents in them externally. The critical phenomena 
of these systems are referred here as "Out of equilibrium 
classes" and discussed in section III. There are also sys- 
tems, which are not related to equilibrium models, in the 
simplest case these are lattice Markov processes of inter- 
acting particle systems (Liggct, 1985). These arc referred 
here as "Genuinely non-equilibrium systems" and are dis- 
cussed in the rest of the work. The discussion of the 
latter type of systems is splitted to three parts: In sec- 
tion IV phase transition classes of simplest such models 
is presented. These are usually reaction-diffusion (RD) 
type of models exhibiting phase transition to absorbing 
states. In section V I list the known classes which occur 
by combinations of basic genuine class processes. These 
models are coupled multi-component RD systems. While 
the former two sections are related to critical phenom- 
ena near extinction in section VI I discuss universality 
classes in systems, where site variables are non- vanishing, 
in surface growth models. The bosonic field theoretical 
description is applicable for them as well. I point out 
mapping between growth and RD systems when it is pos- 
sible. In section II I briefly touch the point of discontin- 
uous nonequilibrium phase transitions especially because 
dynamical scaling may occur at such points. 

I define a critical universality class by the complete 
set of exponents at the phase transition. Therefore dif- 
ferent dynamics split up the basic static classes of ho- 
mogeneous systems. I emphasize the role of symmetries, 
boundary conditions which affect these classes. I also 
point out very recent evidence according to which in low 
dimensional systems symmetries arc not necessarily the 
most relevant factors of universality classes. Although 
the systems covered here might prove to be artificial to 
experimentalists or to application oriented people they 
constitute the fundamental blocks of understanding of 
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nonequilibrium critical phenomena. Note that the under- 
standing of even so simple models runs into tremendous 
difficulties very often. 

I shall not discuss the critical behavior of quantum 
systems (Racz, 2002), self-organized critical phenomena 
(Bak et al., 1987) and neither show experimental real- 
izations. The discussion of the applied methods is also 
omitted due to the lack of space, although in section 
(I.E) I give a brief introduction of the field theoretical 
approach. This section shows the formalism for defin- 
ing nonequilibrium models. This is necessary to express 
the symmetry relations affecting critical behavior. Re- 
searchers from other branches of science are provided a 
kind of catalog of classes in which they can identify their 
models and find corresponding theories. To help navigat- 
ing in the text and in the literature I provided a list of 
the most common abbreviations at the end of the text. 

Besides scaling exponents and scaling relations there 
are many other interesting features of universality classes 
like scaling distribution, extremal statistics, finite size ef- 
fects, statistics of fluctuations in surface growth models 
etc., which I do not discuss in this review. Still I believe 
the shown material provides a useful frame for orientation 
in this huge field. There is no general theory of nonequi- 
librium phase transitions, hence a widespread overview 
of known classes can help theorists deducing the relevant 
factors determining universality classes. 

There are two recent, similar reviews available. One 
of them is by (Marro and Dickman, 1999), which gives 
a pedagogical introduction to driven lattice gas sys- 
tems and to fundamental particle systems with absorb- 
ing states. The other one (Hinrichsen, 2000a) focuses 
more on basic absorbing state phase transitions, methods 
and experimental realizations. However the field evolves 
rapidly and since the publication of this two remarkable 
introductory works a series of new developments have 
come up. The present work aims to give a comprehen- 
sive overview of known nonequilibrium dynamical classes, 
incorporating surface growth classes, classes of spin mod- 
els, percolation and multi-component system classes and 
damage spreading transitions. The relations and map- 
pings of the corresponding models are pointed out. The 
effects of boundary conditions, long-range interactions 
and disorder are shown systematically for each class were 
it is known. Since a debate on the conditions of the par- 
ity conserving class has not been settled yet I provide a 
discussion on it through some surface catalytic model ex- 
amples. Naturally this review can not be complete and I 
apologize for the omitted references. 



A. Critical exponents of equilibrium systems 

In this section I briefly summarize the definition of 
well known critical exponents of homogeneous equilib- 
rium systems and show some scaling relations (Amit, 
1984; Fisher, 1967; Kadanoff and al., 1967; Ma, 1976; 
Stanley, 1971). The basic exponents are deflned via the 



scaling laws: 

CH OC a^l((|r-Te|/Te)-"--l) , (1) 



m OC {Tc-Tf , (2) 

X cx |T-Te|-T , (3) 

m a H^/^" , (4) 

G^^^r) OC r^-'^-"" , (5) 

^ OC IT-Tcl-"^ . (6) 



Here ch denotes the specific heat, m the order param- 
eter, X the susceptibility and £, the correlation length. 
The presence of another degree of freedom besides the 
temperature T, like a (small) external field (labeled by 
H), leads to other interesting power laws when H^O. 
The d present in the expression of two-point correlation 
function Gc {r) is the space dimension of the system. 

Some laws are valid both to the right and to the left of 
the critical point; the values of the relative proportional- 
ity constants, or amplitudes, arc in general different for 
the two branches of the functions, whereas the exponent 
is the same. However there are universal amplitude rela- 
tions among them. We can see that there are altogether 
six basic exponents Nevertheless they are not indepen- 
dent of each other, but related by some simple scaling 
relations: 

an + 2/3 + 7 = 2, an + /3(fe + 1) = 2, (7) 
(2 - r]a)i'± = J, i'±d =2- an ■ 

The last relation is a so-called hyper-scaling law, which 
depends on the spatial dimension d and is not valid above 
the upper critical dimension dc, for example by the Gaus- 
sian thcxjry. Therefore below d^ there arc only two inde- 
pendent exponents in equilibrium. One of the most in- 
teresting aspects of second order phase transitions is the 
so-called universality, i.e., the fact that systems which 
can be very different from each other share the same set 
of critical indices (exponents and some amplitude ratios). 
One can therefore hope to assign all systems to classes 
each of them being identified by a set of critical indices. 



B. Static percolation cluster exponents 

Universal behavior may occur at percolation (Grim- 
mett, 1999; Stauffer and Aharony, 1994), which can be 
considered a purely geometrical phenomenon describing 
the occurrence of infinitely large connected clusters on 
lattices. On the other hand such clusters emerge at 
critical phase transition of lattice models indeed. The 
definition of connected clusters is not unambiguous. It 
may mean the set of sites or bonds with variables in 
the same state, or connected by bonds with probability 
6 = 1 - exp(-2J/fcT). 

By changing the system control parameters {p pp) 
(that usually is the temperature in equilibrium systems) 



4 



the coherence length between sites may diverge as 

m^\p-ppr^ , (8) 

hence percolation at Pp like standard critical phenomena 

exhibits rcnornializibility and universality of critical ex- 
ponents. At Pp the cluster size (s) distribution follows 
the scaling law: 

n.cxs-V(b-Pp|s") . (9) 

while moments of this distribution exhibit singular be- 
havior the exponents: 

^sn,(p) oc Ip-PpI^" , (10) 
Y,^^ris{p) cx \p-pp\-^- . (11) 

s 

Further critical exponents and scaling relations among 
them are shown in (Stauffer and Aharony, 1994). In case 
of completely random placement of (sites, bonds, etc) 
variables (with probability p) on lattices we find random 
isotropic (ordinary) percolation (see Sect.IV.B.l). 
Percolating clusters may arise at critical, thermal transi- 
tions or by nonequilibrium processes. If the critical point 
{Pc) of the order parameter does not coincide with pp 
than at the percolation transition the order parameter 
coherence length is finite and does not influence the per- 
colation properties. We observe random percolation in 
that case. In contrast if Pp = Pc percolation is influ- 
enced by the order parameter behavior and we find dif- 
ferent, correlated percolation universality (Coniglio 
and Klein, 1980; Fortuin and Kasteleyn, 1972; Stauffer 
and Aharony, 1994) whose exponents may coincide with 
those of the order parameter. 

According to the Fortuin-Kasteleyn construction of 
clusters (Fortuin and Kasteleyn, 1972) two nearest- 
neighbor spins of the same state belong to the same 
cluster with probability b — 1 ~ exp(— 2J/A:T). It was 
shown that using this prescription for Zn and 0{n) and 
symmetric models (Bialas and al., 2000; Blanchard and 
al., 2000; Coniglio and Klein, 1980; Fortunato and Satz, 
2001) the thermal phase transition point coincides 
with the percolation limits of such clusters. On the other 
hand in case of "pure-site clusters" {b= 1) different, uni- 
versal cluster exponents arc reported in two dimensional 
models (Fortunato, 2002; Fortunato and Satz, 2001) (see 
Sects. III.A.l, III.B.l, III.D.l). 

C. Dynamical critical exponents 

Nonequilibrium systems were first introduced to study 

relaxation to equilibrium states (Halperin and Hohen- 
berg, 1977) and phase ordering kinetics (Binder and 
Stauffer, 1974; Marro et al, 1979). Power-law time de- 
pendences were investigated away from the critical point 
as well, example by the domain growth in a quench to 



T = 0. Later the combination of different heat-baths, dif- 
ferent dynamics, external currents became popular inves- 
tigation tools of fully nonequilibrium models. To describe 
the dynamical behavior of a critical system additional 
exponents were introduced. For example the relation of 
the divergences of the relaxation time r and correlation 
length ^ is described by the dynamical exponent Z 

TOC^^ . (12) 

Systems out of equilibrium may show anisotropic scaling 
of two (and n) point functions 

G{br,b'^t) =b-^='G{r,t) (13) 

where r and t denote spatial and temporal coordinates, x 
is the scaling dimension and C is the anisotropy exponent. 
As a consequence the temporal (z/||) and spatial {i'±) cor- 
relation length exponents may be different, described by 

C = z. 

Z = C = ^ . (14) 

For some years it was believed that dynamical critical 
phenomena are characterized by a set of three critical 
exponents, comprising two independent static exponents 
(other static exponents being related to these by scaling 
laws) and the dynamical exponent Z. Recently, it was 
discovered that there is another dynamical exponent, the 
'non-equilibrium' or short-time exponent A, needed to 
describe two-time correlations in a spin system ({si}) of 
size L relaxing to the critical state from a disordered 
initial condition (Huse, 1989; Janssen et al., 1989). 

A{t, 0) = ^<J2 s^{0)si{t) >oc i-^/^ (15) 

i 

More recently the persistence exponents 9i and 
9g were introduced by (Derrida et al., 1994; Majumdar 
et al., 1996). These are associated with the probability, 
p{t), that the local or global order parameter has not 
changed sign in time t following a quench to the critical 
point. In many systems of physical interest these expo- 
nents decay algebraically as 

p{t) oc t-^ (16) 

(see however example Sect. V.A). It turned out that in 
systems where the scaling relation 

OgZ = X-d+l-r]a/2 (17) 

is satisfied the dynamics of the global order parame- 
ter is a Markov process. In contrast in systems with 
non-Markovian global order parameter dg is in general 
a new, non-trivial critical exponent (Majumdar et al, 
1996). For example it was shown that while in the d = 1 
Glauber Ising model the magnetization is Markovian and 
the scaling relation (17) is fulfilled, at the critical point 
of the d = 1 NEKIM Ising model this is not satisfied 
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and the persistence behavior there is characterized by a 
different, non-trivial 6g exponent (Menyhard and Odor, 
1997) (see discussion in section IV. D. 2). As we can see 
the universality classes of static models are split by the 
dynamical exponents. 

D. Critical exponents of spreading processes 

In the previous section I defined quantities describing 
dynamical properties of the bulk of a material. In a dual 
way to this one may also consider cluster properties aris- 
ing by initiating a process from an ordered (correlated) 
state with a small cluster of activity. Here I define a 
basic set of critical exponents that occur in spreading 
processes and show the scaling relations among them. In 
such processes phase transition may exist to absorbing 
state(s) where the density of spreading entity (particle, 
agent, epidemic etc.) disappears. The order parameter 
is usually the density of active sites {sj} 

i 

which in the supercritical phase vanishes as 

p°°cx\p-pcf, (19) 

as the control parameter p is varied. The "dual" quan- 
tity is the ultimate survival probability Poo of an infinite 
cluster of active sites that scales in the active phase as 

Poo(x\p-pcf (20) 

with some critical exponent /?' (Grassberger and de la 
Torre, 1979). In field theoretical description of such pro- 
cesses /3 is associated with the particle annihilation, /3' 
with the particle creation operator and in case of time re- 
versal symmetry (see Eq. (88)) they are equal. The crit- 
ical long-time behavior of these quantities are described 
by 

p{t) oc t-^fiAt^/^ii) , P{t) cx t-^giAt^/^ii) , (21) 

where a and 6 are the critical exponents for decay and 
survival, A = \p — pc\, f and g are universal scaling func- 
tions (Grassberger and de la Torre, 1979; Jansscn, 2003; 
Muhoz et al., 1997). The obvious scaling relations among 
them are 

a = /3/z/||, S = p'/uii. (22) 

For finite systems (of size N = L'^) these quantities scale 
as 

p{t) oc r'^/'^ii /'(AtV-ii^ t<i/z/^) ^ (23) 
P{t) cx t-^'/^w g^At^/^w^f^'^ IN) . (24) 

For "relatively short times" or for initial conditions with 
a single active seed the the number of active sites N{t) 



and its mean square of spreading distance [xi) from the 
origin 

R\t) = ^^{Y,xl{t)) (25) 

follow the "initial slip" scaling laws (Grassberger and 
de la Torre, 1979) 

N{t)(xt'' , (26) 

R'^it) cx t' , (27) 
and usually the z = 2/Z relation holds. 

1. Damage spreading exponents 

Phase transitions between chaotic and non-chaotic 
states may be described by damage spreading (DS). 
While DS was first introduced in biology (Kauffman, 
1969) it has become an interesting topic in physics as 
well (Creutz, 1986; Derrida and Weisbuch, 1987; Stanley 
et al., 1986). The main question is if a damage introduced 
in a dynamical system survives or disappears. To investi- 
gate this the usual technique is to make replica(s) of the 
original system and let them evolve with the same dy- 
namics and external noise. This method has been found 
to bo very useful to measure accurately dynamical ex- 
ponents of equilibrium systems (Grassberger, 1995a). It 
has turned out however, that DS properties do depend 
on the applied dynamics. An example is the case of the 
two-dimensional Ising model with heat-bath algorithm 
versus Glauber dynamics (Grassberger, 1995b; Jan and 
de Arcangelis, 1994; Mariz et al, 1990). 

To avoid dependency on the dynamics a definition of 
"physical" family of DS dynamics was suggested by (Hin- 
richsen et al., 1997) according to the active phase may be 
divided to a sub-phase in which DS occurs for every mem- 
ber of the family, another sub-phase where the damage 
heals for every member of the family and a third possible 
sub-phase, where DS is possible for some members and 
the damage disappears for other members. The family 
of possible DS dynamics is defined to be consistent with 
the physics of the single replicas (symmetries, interaction 
ranges etc.). 

Usually the order parameter of the damage is the Ham- 
ming distance between replicas 

where s{i) and s'{i) denote variables of the replicas. At 
continuous DS transitions D exhibits power-law singu- 
larities as physical quantities at the critical point. For 
example one can follow the fate of a single difference be- 
tween two (or more) replicas and measure the spreading 
exponents: 

D{t) cx , (29) 
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Similarly the survival probability of damage variables be- 
haves as: 



(30) 



and similarly to eq.(25) the average mean square spread- 
ing distance of damage variables from the center scales 
as: 



(31) 



Grassbcrger conjectured, that DS transitions should be- 
long to DP class (sec Sect. IV. A) unless they coincide 
with other transition points and provided the probabil- 
ity for a locally damaged state to become healed is not 
zero (Grassberger, 1995c). This hypothesis has been con- 
firmed by simulations of many different systems. 



E. Field theoretical approach to reaction-difTusion 
systems 

In this review I define noncquilibrium systems for- 
mally by their field theoretical action where it is pos- 
sible. Therefore in this subsection I give a brief introduc- 
tion to (bosonic) field theoretical formalism. This will be 
through the simplest example of reaction-diffusion sys- 
tems, via the: A + A ^ (d annihilating random walk 
(ARW) (see Sect.IV.C.l). Similar stochastic differential 
equation formalism can also be set up for growth pro- 
cesses in most cases. For a more complete introduction 
see (Cardy, 1996, 1997; Tauber, 2003). 

A proper field theoretical treatment should start from 
the master equation for the microscopic time evolution 
of probabilities p{a; t) of states a 



dp{a; t) 
dt 







.f3p{a;t) , (32) 



where Ra^0 denotes the transition matrix from state 

a to [3. In field theory this can be expressed in Fock 
space formalism with annihilation (aj) and creation (cj) 
operators satisfying the commutation relation 



(33) 



The states are built up from the vacuum |0 > as the 
linear superposition 

^{t)=Y,P{nun2,...;t)cTcT.:\^> , (34) 



with occupation number coefficients p(ni, n2, t). The 
evolution of states can be described by a Schrodinger-like 
equation 



dt 



(35) 



with a generally non-hermitian Hamiltonian, which in 
case of the ARW process looks like 

H = D J2{ci - Cj){ai - aj) - A ^(a| - c|a|) , (36) 



here D denotes the diffusion strength and A the annihi- 
lation rate. By going to the continuum limit this turns 
into 



H 



d'*x[£)(VV')(V(/))-A((/>2-V'' (/>')] , (37) 



and in the path integral formalism over fields (j){x,t), 
ip{x,t) with weight e"'^^'^''^) one can define an action, 
that in case of ARW is 

S = j dtd'^x [ii}dt(t> + DVV'Vf/) - A(^2 - V^<^^)] • (38) 

The action is analyzed by renormalization group (RG) 

methods at criticality (Amit, 1984; Ma, 1976), usually 
by perturbative epsilon expansion below the upper crit- 
ical dimension dc — that is the lower limit of the valid- 
ity of the mean-field (MF) behavior of the system. The 
symmetries of the model can be expressed in terms of 
the <f>{x,t) field and ip{x,t) response field variables and 
the corresponding hyper-scaling relations can be derived 
(Janssen, 2003; Munoz et al., 1997). 

By a Gaussian transformation one may set up an al- 
ternative formalism — integrating out the response field 
— the Langevin equation, that in case of ARW is 



dt(j) = DV^cj) - 2X(j)^ + ri{x,t) 



(39) 



with a Gaussian noise, exhibiting correlations 

< rjix, t)ri{x', t') >= -X(t>'^6{x - x')6{t - t') . (40) 

Here 5 denotes the Dirac delta functional and A is the 
noise amplitude. From the Langevin equation - if it ex- 
ists - one can deduce a naive upper critical dimension 
{dc) by power counting. However this estimate may be 
modified by fiuctuations, which can be analyzed by the 
application of the RG method. 

II. SCALING AT FIRST ORDER PHASE TRANSITIONS 

In noncquilibrium systems dynamical scaling of vari- 
ables may occur even when the order parameter jumps 
at the transition. We call such a transition first order, 
although the free energy is not defined. First order phase 
transitions have rarely been seen in low dimensions. This 
is due to the fact that in lower dimension the fiuctuations 
are more relevant and may destabilize the ordered phase. 
Therefore fiuctuation induced second ordered phase tran- 
sitions are likely to appear. Hinrichsen advanced the 
hypothesis (Hinrichsen, 2000b) that first-order transi- 
tions do not exist in 1-1-1 dimensional systems without 
extra symmetries, conservation laws, special boundary 
conditions or long-range interactions (which can be gen- 
erated by macroscopic currents or anomalous diffusion in 
noncquilibrium systems for instance). Examples are the 
Glauber and the NEKIM Ising spin systems (see sections 

III. A, IV. D. 2) possessing Z2 symmetry in one dimension 



7 









0s' 


A 


Tjs 


5, 


Zs 


Glauber 





1/2 


.99(2) 


1/2 


.0006(4) 


.500(5) 


1 


NEKIM 


.00(2) 


.444 


.45(1) 


.49(1) 


.288(4) 


.287(3) 


1.14 



TABLE I Critical Id Ising spin exponents at the Glauber and 
NEKIM transition points (Menyhard and Odor, 1998) 



(Glauber, 1963; Menyhard and Odor, 1998), where the 
introduction of a "temperature like" flip inside of a do- 
main or an external field {h) causes discontinuous jump 
in the magnetization order parameter (m) . Interestingly 
enough the correlation length diverges at the transition 
point: ^ oc p^^^ and static 



as well as cluster critical exponents can be defined: 

P,{t,h) oc t-^' 
Rl{t,h) at^» 
\m{t,h) -m{0)\ cxt''' 



(41) 



lim Ps{t,h) oc h/^'" 



(42) 
(43) 

(44) 
(45) 



Here s refers to the spin variables. The table II summa- 
rizes the results obtained for these transitions. Other ex- 
amples for first order transition are known in driven diffu- 
sive systems (Jansscn and Schmittman, 1986), in the Id 
asymmetric exclusion process (Derrida, 1998), in bosonic 
annihilation-fission models (Sect.V.F), in asymmetric 
triplet and quadruplet models (Odor, 2003a) (Sect. IV. F) 
and DCF models for Da > Db and d > 1 (Oerding et ai, 
2000) (Sect.V.H). By simulations it is quite difficult to 
decide whether a transition is really discontinuous. The 
order parameter of weak first order transitions - where 
the jump is small - may look very similar to continuous 
transitions. The hysteresis of the order parameter that is 
considered to be the indication of a first order transition 
is a demanding task to measure. There are some exam- 
ples with debates over the order of the transition (see for 
example (Dickman and Tome, 1991; Hinrichsen, 2000b, 
2001a; Lipowski, 1999; Lipowski and Lopata, 1999; Szol- 
noki, 2000; Tome and de Oliviera, 1989)). In some cases 
the mean-field solution results in first order transition 
(Menyhard and Odor, 1995; Odor et ai, 1993). In two 
dimensions there are certain stochastic cellular automata 
for which systematic cluster mean-field techniques com- 
bined with simulations made it possible to prove first 
oder transitions (Odor and Szolnoki, 1996) firmly (see 
table II). 



III. OUT OF EQUILIBRIUM CLASSES 

In this section I begin introducing basic nonequilibrium 

classes starting with the simplest dynamical extensions 
of equilibrium models. These dynamical systems exhibit 



n 


y = l 

Pc 


y = 

Pc 


= 2 

P(Pc) 


y = 

Pc 


= 3 

p{Pc) 


1 


0.111 


0.354 


0.216 


0.534 


0.372 


2 


0.113 


0.326 


0.240 


0.455 


0.400 


4 


0.131 


0.388 


0.244 


0.647 


0.410 


simulation 


0.163 


0.404 


0.245 


0.661 


0.418 



TABLE 11 Convergence of the critical point estimates in 
various (y = 1, 2, 3) two dimensional SCA calculated by n- 
cluster (GMF) approximation (see section 1V.A.2). First or- 
der transitions are denoted by boldface numbers. The gap 
sizes (p(pc)) of the order parameter shown for y = 2,3 in- 
crease with n, approximating the simulation value (Odor and 
Szolnoki, 1996). 



hermitian Hamiltonian and starting from a nonequilib- 
rium state they evolve into a Gibbs state. Such nonequi- 
librium models include, for example, phase ordering sys- 
tems, spin glasses, glasses etc. In these cases one is usu- 
ally interested in the 'nonequilibrium dynamics' at the 
equilibrium critical point. 

It is an important and new universal phenomenon that 
scaling behavior can be observed far away from critical- 
ity as well. In quenches to zero temperature of model 
A systems (which do not conserve the order parameter) 
the characteristic length in the late time regime grows 
with an universal power-law ^ oc t^^^, while in case of 
model B systems (which conserve the order parameter) 
^ oc t^/^. In model C systems a conserved secondary 
density is coupled to the non-conserved order parame- 
ter. Such models may exhibit model A behavior or 
^ oc depending on the model parameters. 

The effects of such conservation laws in critical systems 
without hermitian Hamiltonian have also been investi- 
gated and will be discussed in later sections. 

Since percolation is a central topic in reaction-diffusion 
systems discussed in Sects. IV and V, for the sake of 
completeness I show recent percolation results obtained 
for systems of Section III too. Another novel feature of 
dynamic phase transitions is the emergence of a chaotic 
state, therefore I shall discuss damage spreading transi- 
tions and behavior in these systems. 

Then I continue towards such nonequilibrium models 
which have no hermitian Hamiltonian and equilibrium 
Gibbs state. In this section I show cases when this is 
achieved by combining different competing dynamics (for 
example by connecting two reservoirs with different tem- 
peratures to the system) or by generating current from 
outside. Field theoretical investigations have revealed 
that model A systems arc robust against the intro- 
duction of various competing dynamics, which are lo- 
cal and do not conserve the order parameter (Grinstein 
et al., 1985). Furthermore it was shown, that this ro- 
bustness of the critical behavior persists if the compet- 
ing dynamics breaks the discrete symmetry of the sys- 
tem (Bassler and Schmittman, 1994) or it comes from 
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reversible mode coupling to a non-critical conserved field 
(Tauber and Racz, 1997). On the other hand if a com- 
peting dynamics is coupled to model B systems by an 
external drive (Schmittman and Zia, 1996) or by local, 
anisotropic order-parameter conserving process (Bassler 
and Racz, 1994, 1995; Schmittman, 1993; Schmittman 
and Zia, 1991) long-range interactions are generated in 
the steady state with angular dependence. The univer- 
sality class will be the same as that of the kinetic version 
of the equilibrium Ising model with dipolcir long-range 
interactions. 

As the number of neighboring interaction sites de- 
creases by lowering the spatial dimensionality of a sys- 
tem with short-ranged interactions the relevance of fluc- 
tuations increases. In equilibrium models finite-range in- 
teractions cannot maintain long-range order in d < 2. 
This observation is known as the Landau-Peierls argu- 
ment (Landau and Lifshitz, 1981). According to the 
Mermin- Wagner theorem (Mermin and Wagner, 1996), 
for systems with continuous symmetry long-range order 
do not exist even in d = 2. Hence in equilibrium mod- 
els phase transition universality classes exist for > 2 
only. One of the main open questions to be answered 
is whether there exist a class of noncquilibrium systems 
with restricted dynamical rules for which the Landau- 
Peierls or Mermin- Wagner theorem can be applied. 

A. Ising classes 

The equilibrium Ising model was introduced by (Ising, 
1925) as the simplest model for an uniaxial magnet but 
it is used in different settings for example binary by flu- 
ids or alloys as well. It is defined in terms of spin vari- 
ables Sj = ±1 attached to sites i of some lattice with the 
Hamiltonian 

H = -jY,SiSi.-BY,Si (46) 

where J is the coupling constant and B is the external 
field. In one and two dimensions it is solved exactly (On- 
sager, 1944), hence it plays a fundamental test-ground 
for understanding phase transitions. The Hamiltonian 
of this model exhibits a global, so called Z2 (up-down) 
symmetry of the state variables. While in one dimension 
a first order phase transition occurs at T = only (see 
section II) in two dimensions there is a continuous phase 
transition where the system exhibits conformal symmetry 
(Henkel, 1999) as well. The critical dimension is dc = 4. 
The following table (III) summarizes some of the known 
critical exponents of the Ising model. The quantum ver- 
sion of the Ising model, which in the simplest cases might 
take the form (in Id) 

H = -Jj2itaf + atat+^ + hat) , (47) 

i 

- where cr^'^ are Pauli matrices and t and h are couplings 

- for T > has been shown to exhibit the same critical 



exponent 


d = 2 


d = 3 


d = 4(MF) 


an 


O(log) 


0.1097(6) 





13 


1/8 


0.3265(7) 


1/2 


7 


7/4 


1.3272(3) 


1 




1 


0.6301(2) 


1/2 



TABLE III Static critical exponents of the Ising model 



behavior as the classical one (in the same dimension). 
For T = however quantum effects become important 
and the quantum Ising chain can be associated with the 
two dimensional classical Ising model such that the trans- 
verse field t plays the role of the temperature. In general 
a mapping can be constructed between classical d + 1 
dimensional statistical systems and d dimensional quan- 
tum systems without changing the universal properties, 
which has been widely utilized (Frdakin and Susskind, 
1978; Suzuki, 1971). The effects of disorder and bound- 
ary conditions are not discussed here (for recent reviews 
see (Alonso and Munoz, 2001; Igloi et al., 1993)). 



1. Correlated percolation clusters at Tc 

If we generate clusters in such a way that we join near- 
est neighbor spins of the same sign we can observe per- 
colation at Tc in 2d. While the order parameter percola- 
tion exponents /3p and 7p of this percolation (defined in 
section IV.B.l) were found to be different from the expo- 
nents of the magnetization (/3, 7) the correlation length 
exponent is the same: u = Vp. For 2d models with Z2 
symmetry the universal percolation exponents are (For- 
tunato, 2002; Fortunate and Satz, 2001): 

(3p = 0.049(4), 7p = 1.908(16) . (48) 

These exponents are clearly different from those of the 
ordinary percolation classes (Table XV) or from Ising 
class magnetic exponents. 

On the other hand by Fortuin-Kasteleyn cluster con- 
struction (Fortuin and Kasteleyn, 1972) the percolation 
exponents of the Ising model at Tg coincide with those of 
the magnetization of the model. 



2. Dynamical Ising classes 

Kinetic Ising models such as the spin-flip Glauber Ising 
model (Glauber, 1963) and the spin-exchange Kawasaki 
Ising model (Kawasaki, 1966) were originally intended 
to study relaxational processes near equilibrium states. 
In order to assure the arrival to an equilibrium state the 
detailed balance condition for transition rates {wi^j) and 
probability distributions {P{s)) are required to satisfy 

Wi^jP{s{i))=Wj^,P{s{3)) . (49) 
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Knowing that Peq{s) oc exp(— if(s)/(fcsT)) this entails 
the 



eM-^H{s)/{kBT)) 



(50) 



condition which can be satisfied in many different ways. 
Assuming spin-flips (which do not conserve the magne- 
tization (model A)) Glauber formulated the most gen- 
eral form in a magnetic field (h) 



Wi^ =Wi{l — tanh /iSj) w Wj(l — /is,) 



(51) 
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2.165(10) 


2.325(10) 


2.032(4) 
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0.737(1) 


0.667(8) 


1.362(19) 
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1/4 




0.225(10) 




0.41(2) 


1/2 





TABLE IV Critical dynamical exponents in the Ising model. 
Columns denoted by (A) and (B) refer, model A and model B 
dynamics. Data are from (Grassberger, 1995a; Jaster et ai, 
1999; Stauffer, 1996; Zheng, 1998, 2000, 2001; Zwerger, 1981) 



Wi = + 6si-iSi+i) - -Si{si-i + Si+i 



)) (52) 



where 7 = tanh 2 J/fcT, F and 6 are further parameters. 
The d = 1 Ising model with Glauber kinetics is exactly 
solvable. In this case the critical tcmpcratiirc is at T = 
and the transition is of first order. We recall that pr = 
e~kT plays the role of ^-^^ in Id and in the vicinity of 
T = critical exponents can be defined as powers of pt, 
thus e.g. that of the coherence length, i^±, via ^ oc Pt~^^ 
(see Section II). In the presence of a magnetic field B, 
the magnetization is known exactly. At T = 



m(T = 0,B) = sgn{B). 



(53) 



Moreover, for ^ 1 and B/kT <C 1 the the exact solu- 
tion reduces to 

m~2/i^; h = B/kBT. (54) 

In scaling form one writes: 

mr^r^gih^^) (55) 

where A is the static magnetic critical exponent. Com- 
parison of eqs. (54) and (55) results in f3s = and 
A = i'± . These values are well known for the Id Ising 
model. It is clear that the transition is discontinuous at 
B = 0, also when changing B from positive to negative 
values, see eq.(53). The order of limits are meant as: 
S ^ and then T ^ 0. The S = 0, T = 1 case is usu- 
ally referred as the Glauber-Ising model. The dynamical 
exponents are (Glauber, 1963; Majumdar et al., 1996): 

Zld Glauber = 2 , Glauber = 1/4 (56) 

Applying spin-exchange Kawasaki dynamics, which 
conserves the magnetization (model B) 



Wi 



1 

27 



1 



72 



{si-iSi + Si+isi + 2) 



(57) 



where 72 = tanh(2J/fcBr), the dynamical exponent is 
different. According to linear response theory (Zwcrgcr, 
1981) in one dimension, at the critical point (Tc = 0) it 
is: 



Note however, that in case of fast quenches to T = 
coarsening with scaling exponent 1/3 is reported (Cornell 
et al., 1991). Hence another dynamic Ising universality 
class appears with the same static but different dynami- 
cal exponents. 

Interestingly while the two dimensional equilibrium 
Ising model is solved the exact values of the dynam- 
ical exponents are not known. Table IV summarizes 
the known dynamical exponents of the Ising model in 
d= 1, 2, 3, 4. The d = 4 results are mean-field values. In 
Section IV. D. 2 I shall discuss an another fully nonequi- 
librium critical point of the d = 1 Ising model with com- 
peting dynamics (the NEKIM) , where the dynamical ex- 
ponents break the scaling relation (17) and therefore the 
magnetization is a non-Markovian process. For d > 3 
there is no non-Markovian effect (hence 9 is not indepen- 
dent) but for d = 2,3 the situation is still not completely 
clear (Zheng, 1998). 

In one dimension the domain walls (kinks) between up 
and down regions can be considered as particles. The 
spin-flip dynamics can be mapped onto particle move- 
ment 

Tii^TTi ~ -o^o* (59) 

or to the creation or annihilation of neighboring particles 



TTT^Tit 



(60) 



Therefore the T = Glauber dynamics is equivalent to 
the diffusion limited annihilation (ARW) mentioned al- 
ready in Section. II. By mapping the spin-exchange dy- 
namics in the same way more complicated particle dy- 
namics emerges, example: 



(61) 



■lid Kaw 



(58) 



one particle may give birth of two others or three par- 
ticle may coagulate to one. Therefore these models are 
equivalent to branching and annihilating random walks 
to be discussed in Section IV.D.l. 



3. Competing dynamics added to spin-flip 

Competing dynamics in general break the detailed bal- 
ance symmetry (49) and make the kinetic Ising model to 
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relax to a nonequilibrium steady state (if it exists) . Gen- 
erally these models become unsolvable, for an overview 
see (Racz, 1996). It was argued by (Grinstein et ai, 
1985) that stochastic spin-flip models with two states 
per site and updating rules of a short-range nature with 
Z2 symmetry should belong to the (kinetic) Ising model 
universality class. Their argument rests on the stability 
of the dynamic Ising fixed point in d = 4 — e dimen- 
sions with respect to perturbations preserving both the 
spin inversion and the lattice symmetries. This hypoth- 
esis has received extensive confirmation by Monte Carlo 
simulations (Castro and al., 1998)-(01iveira et ai, 1993), 
(H. W. J. Blote and Zia, 1990; Santos and Teixeira, 1995; 
Tamayo et ai, 1995) as well as from analytic calculations 
(Marques, 1989, 1990; T. Tome and Santos, 1991). The 
models investigated include Ising models with a competi- 
tion of two (or three (Tamayo et ai, 1995)) Glauber-like 
rates at different temperatures (Blote et ai, 1990; Gon- 
zalez-Miranda et ai, 1987; Marques, 1989; T. Tome and 
Santos, 1991), or a combination of spin-flip and spin- 
exchange dynamics (Garrido et ai, 1989), majority vote 
models (Castro and al., 1998; Santos and Teixeira, 1995) 
and other types of transition rules with the restrictions 
mentioned above (Oliveira et ai, 1993). 

Note that in all of the above cases the ordered state is 
fluctuating and non- absorbing. By relaxation processes 
this allows fluctuations in the bulk of a domain. It is 
also possible (and in Id it is the only choice) to gen- 
erate nonequilibrium two-state spin models (with short 
ranged interactions) where the ordered states are frozen 
(absorbing), hence by the relaxation to the steady state 
fluctuations occur at the boundaries only. In this case 
non-Ising universality appears which is called the voter 
model (VM) universality class (see Sect. IV. C). 

In two dimensions general Z2 symmetric update 

rules were investigated by (Achahbar et ai, 1996; Dor- 
nic et ai, 2001; Drouffe and Godreche, 1999; Oliveira 
et ai, 1993). The essence of this - following (Drouffe 
and Godreche, 1999) - is described below. Let us con- 
sider a two dimensional lattice of spins Si — ±1, evolving 
with the following dynamical rule. At each evolution 
step, the spin to be updated flips with the heat bath 
rule: the probability that the spin Si takes the value -1-1 
is P{si = 1) = p{hi), where the local field hi is the sum 
over neighboring sites sj and 



p{h) = ^{l + t&nh[l3{h)h]). 



(62) 



The functions p{h) and f3{h) are defined over integral 
values of h. For a square lattice, h takes the values 4, 2, 
0, —2, —4. We require that p{—h) 1 — p{h), in order 
to keep the up down symmetry, hence f3{—h) = I3{h) and 
this fixes p(0) = 1/2. The dynamics therefore depends 
on two parameters 
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FIG. 1 Phase diagram of 2d, Zi symmetric nonequilibrium 
spin models from (Drouffe and Godreche, 1999). Broken lines 
correspond to the noisy voter model (V), the Ising model (I) 
and the majority vote model (M). The low temperature phase 
is located in the upper right corner, above the transition line 
(full line). 



or equivalently on two effective temperatures 



Ti = 



1 



To 



1 



Defining the coordinate system 
2 



ti 



tanh ■ 



to 



tanh ■ 



with < ti,t2 < 1 this yields 



1 



Pi 



(1 + ti) 



P2 



1 



T2 



2t2 

l + tl 



(64) 



(65) 



(66) 



Pi=p{2), P2=p(4), 



(63) 



with 1/2 < pi,p2 < 1- One can call Ti and T2 as two tem- 
peratures, respectively associated with interfacial noise, 
and to bulk noise. Each point in the parameter plane 
(PijP2)j or alternatively in the temperature plane [ti, t2), 
corresponds to a particular model. The class of models 
thus defined comprises as special cases the Ising model, 
the voter and anti-voter models (Ligget, 1985), as well 
as the majority vote (Castro and al., 1998; Ligget, 1985) 
model (see Fig. 1). The P2 = 1 line corresponds to 
models with no bulk noise (T2 = 0), hence the dynamics 
is only driven by interfacial noise, defined above. The 
Pi — 1 line corresponds to models with no interfacial 
noise (Ti = 0), hence the dynamics is only driven by 
bulk noise. In both cases effects due to the curvature of 
the interfaces is always present. For these last models, 
the local spin aligns in the direction of the majority of its 
neighbors with probability one, if the local field is h — 2, 
i.e. if there is no consensus amongst the neighbors. If 
there is consensus amongst them, i.e. ii h = 4, the local 
spin aligns with its neighbors with a probability p2 < 1- 
Simulations (Oliveira et ai, 1993) revealed that the 
transition line between the low and high temperature re- 
gions is Ising type except for the endpoint {pi = l,p2 = 
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0.33(2) 


1.16(6) 


0.13(4) 


0.62(3) 



TABLE V Critical exponents of the d = 2 randomly driven 
lattice gas 



0.75), that is first-order and corresponds to the voter 
model class. The local persistence exponent was also 
found to be constant along the line 9i ^ 0.22 (Drouffe and 
Godreche, 1999) in agreement with that of the A-model 
(see Table IV) except for the VM point. The dynamics of 
this class of models may be described formally in terms 
of reaction diffusion processes for a set of coalescing, an- 
nihilating, and branching random walkers (Drouffe and 
Godreche, 1999). There are simulation results for other 
models exhibiting absorbing ordered state indicating VM 
critical behavior (Hinrichsen, 1997; Lipowski and Droz, 
2002a). 

It is important that nontrivial, nonequilibrium phase 
transition may occur even in one dimension if spin- 
exchange is added to spin- flip dynamics, the details will 
be discussed in Section IV. D. 2. 



4. Competing dynamics added to spin-exchange 

As mentioned in Section III model-B systems are 
more sensitive to competing dynamics. Local and 
anisotropic order parameter conserving processes gener- 
ate critical behavior that coincides with that of the ki- 
netic dipolar-interaction Ising model. In two dimensions 
both simulations and field theory (Praestraad and al., 
1994; Praestraad et al., 2000) predicts the critical ex- 
ponents The critical dimension is dc = 3. It is shown 
that the Langevin equation (and therefore the critical 
behavior) of the anisotropic diffusive system coincides 
with that of the randomly driven lattice gas system as 
well. Other systems in this universality class are the 
two-temperature model (Garrido et al., 1990), the ALGA 
model (Binder, 1981) and the infinitely fast driven lat- 
tice gas model (Achahbar et al., 2001). In the randomly 
driven lattice gas model particle current does not occur 
but an anisotropy can be found, therefore it was argued 
(Achahbar et al., 2001) that the particle current is not a 
relevant feature for this class. This argument gives the 
possibility to understand why some set of simulations 
of driven lattice systems (Valles and Marro, 1987) leads 
to different critical behavior than that of the canonical 
coarse-grained representative of this class, in which an 
explicit particle current jx is added to the continuous, 
model-B Ising model Hamiltonian: 



dt 



VC 



(67) 



(Marro and Dickman, 1999) (here 77 is a parameter, C, is 
a Gaussian noise). In this model one obtains mean- field 



exponents for 2 > d > 5 (with weak logarithmic correc- 
tions at d = 2) with (3=1/2 exactly. To resolve contra- 
dictions between simulation results of (Valles and Marro, 
1987) and (Leung, 1991) ref. (Zia et al, 2000) raised 
the possibility of the existence of another, extraordinary, 
"stringy" ordered phase in Ising type driven lattice gases 
that may be stable in square systems. 



5. Long-range interactions and correlations 

Universal behavior is due to the fact at criticality long- 
range correlations are generated that make details of 
short ranged interactions irrelevant. However one c;an 
also investigate the scaling behavior in systems with long- 
range interactions or with dynamically generated long- 
range correlations. If the Glauber Ising model (with non- 
conserving dynamics) changed to a nonequilibrium one 
in such a way that one couples a nonlocal dynamics 
(Droz et al., 1990) to it long-range isotropic interactions 
are generated and mean-field critical behavior emerges. 
For example if the nonlocal dynamics is a random Levy 
flight with spin exchange probability distribution 



P(r) cx 



1 



(68) 



effective long range interactions of the form Vef / oc r~'^~'^ 
are generated and the critical exponents change continu- 
ously as the function a and d (Bergersen and Racz, 1991). 
Similar conclusions for other nonequilibrium classes will 
be discussed later (Sect. IV. A. 6). 

The effect of power-law correlated initial conditions 



(</.(O)0(r)) 



-{d-a) 



in case of a quench to the or- 



dered phase in systems with non-conserved order 

parameter was investigated by (Bray et al., 1991). An 
important example is the (2-|-l)-dimensional Glauber- 
Ising model quenched to zero temperature. It was ob- 
served that long-range correlations are relevant if a ex- 
ceeds a critical value Uc. Furthermore, it was shown 
that the relevant regime is characterized by a continu- 
ously changing exponent in the autocorrelation function 
A{t) = [(j){r,t)(j){r,Q)] ~ t-(d-<^)/4^ whereas the usual 
short-range scaling exponents could be recovered below 
the threshold. These are in agreement with the simula- 
tions of the two-dimensional Ising model quenched from 
T = to T = 0. 



6. Damage spreading behavior 

The high temperature phase of the Ising model is 
chaotic. By lowering T a non-chaotic phase may emerge 
at Td with different types of transitions depending on the 
dynamics. The dynamics dependent DS critical behav- 
ior in different Ising models is in agreement with a con- 
jecture by (Grassberger, 1995c). Dynamical simulations 
with heat-bath algorithm in 2 and 3 dimensions (Grass- 
berger, 1995a; Gropengiesser, 1994; Wang and Suzuki, 
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1996) resulted in = Tc with a DS dynamical exponent 
coinciding with that of the Z-s of the replicas. In this case 
the DS transition picks up the Ising class universality of 
its replicas. With Glauber dynamics in 2d Tj, < Tc and 
DP class DS exponents were found (Grassberger, 1995b). 
With Kawasaki dynamics in 2d on the other hand the 
damage always spreads (Vojta, 1998). With Swendsen- 
Wmg dynamics in 2d > Tc and DP class DS behavior 
was observed (Hinrichsen et al, 1998). 

In one dimensional, nonequilibrium Ising models it is 
possible to design different dynamics showing either PC 
or DP class DS transition as the function of some con- 
trol parameter. This depends on the DS transition co- 
incides or not with the critical point. In the PC class 
DS case damage variables follow BARW2 dynamics (see 
Sects. IV. D.l and IV. D. 2) and Z2 symmetric absorbing 
states occur (Hinrichsen and Domany, 1997; Odor and 
Menyhard, 1998). 



B. Potts classes 

The generalization of the two-state equilibrium Ising 

model was introduced by (Potts, 1952), for an overview 
see (Wu, 1982). In the g-state Potts model the state 
variables can take q different values Si € (0, 1, 2, ...q) and 
the Hamiltonian is a sum of Kronecker delta function of 
states over nearest neighbors 
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FIG. 2 Schematic plot from (Wu, 1982) for qdd) (solid Hue). 
Open symbols correspond to continuous phase transition, 
filled symbol to a known first order transition. Below the 
dashed line the transition is continuous too. 
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TABLE VI Static exponents of the q states Potts model in 
two dimension 



H = -J S{si-Si.) 



(69) 



This hamiltonian exhibits a global symmetry described 
by the permutation group of q elements {Sq). The Ising 
model is recovered in the q = 2 case (discussed in Sec- 
tion III. A). The g-state Potts model exhibits a disordered 
high-temperature phase and an ordered low-temperature 
phase. The transition is first-order, mean-filed-like for 
q-s above the qdd) curve (shown in Figure 2) (and for 
q > 2 in high dimensions). The q = 1 limit can be shown 
(Fortuin and Kasteleyn, 1972) to be equivalent to the 
isotropic percolation (see Section IV.B.l) that is known 
to exhibit a continuous phase transition with dc = 6. 
The problem of finding the effective resistance between 
two node points of a network of linear resistors was solved 
by Kirchhoff in 1847. Ref. (Fortuin and Kasteleyn, 1972) 
showed that Kirchhoff's solution can be expressed as a 
q = limit of the Potts partition function. Further map- 
pings were discovered between the spin glass (Edwards 
and Anderson, 1975) and the q = 1/2 Potts model and 
the two dimensional g = 3, 4 cases to vertex models (see 
(Baxter, 1982)). From this figure we can see that for 
q > 2 Potts models continuous transitions occur in two 
dimensions only [q = 3,4). Fortunately these models are 
exactly solvable (see (Baxter, 1982)) and exhibit confor- 
mal symmetry as well as topological, Yang-Baxter invari- 
ance. The static exponents in two dimensions are known 
exactly (Table VI) 



1. Correlated percolation at Tc 

In 2d models with Z3 symmetry the critical point coin- 
cides with the percolation of site connected clusters and 
the following percolation exponents are reported (Fortu- 
nate and Satz, 2001): 



f3p = 0.075(14), 7p = 1.53(21) 



(70) 



In case of Fortuin-Kasteleyn cluster construction (For- 
tuin and Kasteleyn, 1972) the percolation exponents of 
the q state Potts model at Tc coincide with those of the 
magnetization of the model. 



2. Dynamical Potts classes 

The model-A dynamical exponents of the Potts 
classes are determined in two dimensions for g = 3, 4 by 
short-time Monte Carlo simulations (de Silva et al, 2002) 
(Table VII). The exponents were found to be the same 
for heat-bath and Metropolis algorithms. For the zero 
temperature local persistence exponent in one dimen- 
sion exact formulas have been determined. For sequential 
dynamics (Derrida et al., 1995) 
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(71) 
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TABLE VII Model-A dynamical exponents of the q states 
Potts model in two dimensions 



while for parallel dynamics Qi^p = 29i_s (Menon and Ray, 
2001). In a deterministic coarsening it is again different 
(see (Bray et al, 1994; Gopinathan, 2001). As we can 
see the dynamical universality class characterized by the 
^ oc t^/"^ characteristic length growth is split by different 
dynamics, as reflected by the persistence exponent. 

Similarly to the Ising model case in nonequilibrium 
systems the Potts model symmetry turned out to be a 
relevant factor for determining the universal behavior of 
transitions to fluctuating ordered states (Brunstein 
and Tome, 1998; Crisanti and Grassberger, 1994; Szabo 
and Czaran, 2001). On the other hand in case of nonequi- 
librium transitions to absorbing states a simulation 
study (Lipowski and Droz, 2002a) suggests first order 
transition for all g > 2 state Potts models in d > 1 di- 
mensions. The q = 2, d = 2 case corresponds to the 
Voter model class (Sect. IV.C) and in Id either PC class 
transition {q = 2) or N-BARW2 class transition {q = 3) 
(Sect.V.K) occurs. 

The DS transition oi a. q = 3, d = 2 Potts model 
with heat-bath dynamics was found to belong to the DP 
class (Sect.IV.A) because Tj, > (da Silva et al., 1997). 



3. Long-range interactions 

The effect of long-range interaction has also been in- 
vestigated in case of the one dimensional q = 3 Potts 
model (Glumac and Uzelac, 1998) with the Hamiltonian 



%<0 ' 



J 



5{Si,Sj) . 



(72) 



For cr < cTc ~ 0.65 a crossover from second order to first 
order (mean-field) transition was located by simulations. 
Similarly to the Ising model here we can expect to see 
this crossover if we generate the long-range interactions 
by the addition of a Levy type Kawasaki dynamics. 



C. XY model classes 

The classical XY model is defined by the Hamiltonian 
H = - jY,cos{Qi-ei.) . (73) 

with continuous 8^ € [0, 2tt] state variables. This model 
has a global U{1) symmetry. Alternatively the XY model 
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TABLE VIII Static exponents of the XY model in three di- 
mensions 



can be defined as a special N = 2 case of 0{N) symmetric 
models such that the spin vectors are two dimensional Sj 
with absolute value S? = 1 



H 



-J 



(74) 



<i,v> 



In this continuous model in two dimensions no local 
order parameter can take zero value according to the 
Mermin- Wagner theorem (Mermin and Wagner, 1996). 
The appearance of free vortexes (which are non-local) 
cause an unusual transition mechanism that implies that 
most of the thermodynamic quantities do not show 
power-law singularities. The singular behavior of the cor- 
relation length (^) and the susceptibility (x) is described 
by the forms for T > 



(75) 



where C is a non- universal positive constant. Conven- 
tional critical exponents cannot be used, but one can de- 
fine scaling dimensions. At Tc the two-point correlation 
function has the following long-distance behavior 



G{r) 



oc r 



-1/4 



(Inr) 



1/8 



(76) 



implying ija = 1/4, and in the entire low-temperature 
phase 



G{r) Qc r""-'^) 



(77) 



such that the exponent ria is a continuous function of 
the temperature, i.e. the model has a line of criti- 
cal points starting from Tc to T = 0. This is the so 
called Kosterlitz-Thouless critical behavior (Kosterlitz 
and Thouless, 1973) and corresponds to the conforinal 
field theory with c = 1 (Itzykson and Drouffe, 1989). 
This kind of transition can experimentally observed in 
many effectively two-dimensional systems with 0(2) sym- 
metry, such as thin films of superfluid helium and de- 
scribes roughening transitions of SOS models at crystal 
interfaces. In three dimensions the critical exponents 
of the 0(N) {N = 0, 1, 2, 3, 4) symmetric field theory have 
been determined by perturbative expansions up to sev- 
enth loop order (Guida and Zinn- Justin, 1998). The Ta- 
ble VIII summarizes these results for the XY case (for a 
more detailed overview see (Pelissetto and Vicari, 2000)). 

Exponents with model A dynamics in 2d have been 
determined by short-time simulations and logarithmic 
corrections to scaling were found (Ying et ai, 2001). Ta- 
ble IX summarizes the known dynamical exponents 
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TABLE IX Model-A dynamical exponents of XY model in 
two dimensions 



1. Long-range correlations 

Similarly to the Ising model (Bassler and Racz, 1995) 
studied the validity of the Mermin- Wagner theorem by 
transforming the two dimensional XY model to a non- 
equilibrium one using two-temperature, model-A dy- 
namics. They found that the Mermin- Wagner theorem 
does not apply for this case since effective long-range in- 
teractions are generated by the local nonequilibrium dy- 
namics. The universality class of the phase transition 
of the model coincides with that of the two temperature 
driven Ising model (see Table V). 

2. Self-propelled particles 

An other XY-like nonequilibrium model that exhibits 
ordered state for d < 2 dimensions is motivated by the de- 
scription of the 'flocking" behavior among living things, 
such as birds, slime molds and bacteria. In the simplest 
version of the self-propelled particle model (Vicsek et al, 
1995) each particle's velocity is set to a fixed magni- 
tude, 7;o. The interaction with the neighboring particles 
changes only the direction of motion: the particles tend 
to align their orientation to the local average velocity. In 
one dimension it is defined on the lattice as 

Xi,{t + 1) = Xt{t) + VQUi{t), 

Ui{t + l) = G({u{t))>j+^i, (78) 

where the particles are characterized by their coordi- 
nate Xi and dimensionlcss velocity w,,; and the function G 
that incorporates both the propulsion and friction forces 
which set the velocity in average to a prescribed value vq: 
G{u) > u for u < 1 and G{u) < u for u > 1. The distri- 
bution function P(x = ^) of the noise is uniform in the 
interval [—rj/2, ri/2]. Keeping vq constant, the adjustable 
control parameters of the model are the average density 
of the particles, p, and the noise amplitude r/. The order 
parameter is the average velocity (j) = {u) which vanishes 
as 

<l>{v,P)-U't(pry f°'-^<^^(^) , (79) 

I for J? > ric{p) 

at a critical rjc{p) value. 

This model is similar to the XY model of classical 
magnetic spins because the velocity of the particles, like 
the local spin of the XY model, has fixed length and 
continuous rotational U(l) symmetry. In the vq = and 
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low noise limit the model reduces exactly to a Monte- 
Carlo dynamics of the XY model. 

A field theory that included in a self-consistent way the 
non-equilibrium effects was proposed by (Tu and Toner, 
1995). They have shown that their model is different 
from the XY model for d < 4. The essential difference 
between the self-propelled particle model and the equi- 
librium XY model is that at different times, the "neigh- 
bors" of one particular " bird" will be different depending 
on the velocity field itself. Therefore, two originally dis- 
tant "birds" can interact with each other at some later 
time. They found a critical dimension dc = 4, below 
which linearized hydrodynamics breaks down, but owing 
to a Galilean invariance they could obtain exact scaling 
exponents ind = 2. For the dynamical exponent they got 
Z = 6/5. Numerical simulations (Czirok et al, 1997; Vic- 
sek et al., 1995) indeed found a long range ordered state 
with a continuous transition characterized by /3 = 0.42(3) 
in two dimensions. 

In one dimension the field theory and simulations 
(Czirok et al., 1999) provided evidence for a continuous 
phase transition with /? = 0.60(5), which is different from 
the mean-field value 1/2 (Stanley, 1971). 

D. 0(N) symmetric model classes 

As already mentioned in the previous section the 0{N) 
symmetric models are defined on spin vectors Sj of unit 
length S? = 1 with the Hamiltonian 

H = -J SiS; . (80) 
<»,»•> 

The most well known of them is the classical Heisenberg 

model that corresponds to iV = 3 being the simplest 
model of isotropic ferromagnets. The A'^ = 4 case cor- 
responds to the Higgs sector of the Standard Model at 
finite temperature. The A'^ = case is related to poly- 
mers and the N = 1 and N = 2 cases are the Ising and 
XY models respectively. The critical dimension is dc = 4 
and by the Mermin- Wagner theorem we cannot find finite 
temperature phase transition in the short range equilib- 
rium models for N > 2 below d = 3. The static critical 
exponents have been determined by e = 4 — d expansions 
up to five loop order (Gorishny et al., 1984), by exact 
RG methods (see (Berges et al, 2002) and the references 
therein), by simulations (Hasenbuch, 2001) and by se- 
ries expansions in 3D (see (Guida and Zinn- Justin. 1998) 
and the references there). In Table X I show the latest 
estimates from (Guida and Zinn- Justin, 1998) in three 
dimensions for N = 0,3.4. The N ^ oo limit is the 
exactly solvable spherical model (Berlin and Kac, 1952; 
Stanley, 1968). For a detailed discussion of the static 
critical behavior of the 0{N) models see (Pelissetto and 
Vicari, 2000). 

The dynamical exponents for model A are known ex- 
actly for the (A^ oo) spherical model (Janssen et al., 
1989; Majumdar et al, 1996) case. For other cases 
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N_ an I3_ 7 u± rjg 

_0_ 0.235(3) 0.3024(1) 1.597(2) 0.588(1) 0.028(2) 
_3_ -0.12(1) 0.366(2) 1.395(5) 0.707(3) 0.035(2) 

4|-0.22(2)| 0.383(4) | 1.45(1) 1 0.741(6) 1 0.035(4) 

TABLE X Static exponents of the 0(N) model in three di- 
mensions 
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TABLE XI Model A dynamical exponents of the 0(N) model 
in three dimensions 



e = 4 — d expansions up to two loop order exist (Majum- 
dar et al, 1996; Oerding et al, 1997). For a discussion 
about the combination of different dynamics see the gen- 
eral introduction Section III and (Tauber et al., 1998b). 



1. Correlated percolation at Tc 

In three dimensions for 0{2), 0(3) and 0(4) symmetric 
models the Fortuin-Kasteleyn cluster construction (For- 
tuin and Kastclcyn. 1972) results in percolation points 
and percolation exponents which coincide with the corre- 
sponding Tc-s and magnetization exponent values (Blan- 
chard and al., 2000). 



IV. GENUINE, BASIC NONEQUILIBRIUM CLASSES 

In this section 1 introduce such "genuine nonequilib- 
rium" universality classes that do not occur in dynam- 
ical generalizations of equilibrium systems. Naturally 
in these models there is no hcrmitian Hamiltonian and 
they are defined by transition rates not satisfying the de- 
tailed balance condition (49). They can be described by 
a master equation and the deduced stochastic action or 
Langevin equation if it exists. The most well known cases 
are reaction-diffusion systems with order-disorder transi- 
tions in which the ordered state may exhibit only small 
fluctuations, hence they trap a system falling in it (ab- 
sorbing state). They may occur in models of population 
(Albano, 1994), epidemics (Ligget, 1985; Mollison, 1977), 
catalysis (ZifF et al., 1986) or enzyme biology (Berry, 
2003) for example. There are also other nonequilibrium 
phase transitions for example in lattice gases with cur- 
rents (Evans, 2000; Evans et al., 1995, 1998; Kolomeisky 
et al, 1998) or in traffic models (Chowdhury et al., 2000), 
but in these systems the critical universality classes have 
not been explored yet. 

Phase transitions in such models may occur in low 
dimensions in contrast with equilibrium ones (Marro 



and Dickman, 1999). As it was already shown in Sec- 
tion III. A. 2 reaction-diffusion particle systems may be 
mapped onto spin-flip systems stochastic cellular au- 
tomata (Chopard and Droz, 1998) or interface growth 
models (see Sect. VI). The mapping however may lead 
to nonlocal systems that does not bear real physical rele- 
vance. The universality classes of the simple models pre- 
sented in this section constitute the fundamental building 
blocks of more complex systems. 

For a long time phase transitions with completely 
frozen absorbing states were investigated. A few uni- 
versality classes of this kind were known (Grassberger, 
1996),(Hinrichsen, 2000a), the most prominent and the 
first one that was discovered is that of the directed perco- 
lation (DP) (Kinzel, 1983). An early hypothesis (Grass- 
berger, 1982a; Grinstein et al, 1989; Janssen, 1981) was 
confirmed by all examples up to now. This DP hy- 
pothesis claims that in one component systems ex- 
hibiting continuous phase transitions to a single 
absorbing state (without extra symmetry and in- 
homogeneity or disorder) short ranged interac- 
tions can generate DP class transition only. De- 
spite the robustness of this class experimental observa- 
tion is still lacking (Grassberger, 1996; Hinrichsen, 2000c) 
probably owing to the sensitivity to disorder that cannot 
be avoided in real materials. 

A major problem of these models is that they are 
usually far from the critical dimension and critical fluc- 
tuations prohibit mean-field (MF) like behavior. Fur- 
ther complication is that bosonic field theoretical meth- 
ods cannot describe particle exclusion that may obvi- 
ously happen in d = 1. The success of the applica- 
tion of bosonic field theory in many cases is the con- 
sequence of the asymptotically low density of particles 
near the critical point. However in multi-component sys- 
tems, where the exchange between different types is non- 
trivial, bosonic field theoretical descriptions may fail. In 
case of the binary production (BP) models (Sect. V.F) 
bosonic RG predicts diverging density in the active phase 
contrary to the lattice model version of hard-core parti- 
cles (Carlon et al, 2001; Hinrichsen, 2001b,c; Odor, 2000, 
2001a; Park et al., 2001). Fermionic field theories on the 
other hand have the disadvantage that they are non-local, 
hence results exist for very simple reaction-diffusion sys- 
tems only (Brunei et al, 2000; Park et al, 2000; Wij- 
land, 2001). Other techniques like independent interval 
approximation (Krapivsky and Ben-Naim, 1997), empty 
interval method (ben Avraham et al, 1990), series ex- 
pansion (Essam et al, 1996) or density matrix renormal- 
ization (DMRG) are currently under development. 

The universal scaling-law behavior in these models is 
described by the critical exponents in the neighborhood 
of a steady state, hence the generalization of dynamical 
exponents introduced for "Out of equilibrium classes" 
(Sect. HI) (like Z, 9, A ...etc.) is used. Besides that 
there are genuinely non-equilibrium dynamical exponents 
as well to characterize spreading behavior, defined in 
Sect. I. D. For each class I discuss the damage spreading 
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TABLE XII Estimates for the critical exponents of directed 
percolation. Data are from: (Jensen, 1999a) (ld),(Voigt and 
Ziff, 1997) (2d), (Jensen, 1992)(3d), (Bronzan and Dash, 1974; 
Janssen, 1981)(4 - e). 



FIG. 3 Directed site percolation in d = 1 + 1 dimensions 



transitions, the effects of different boundary conditions, 
disorder and long-correlations generated by anomalous 
diffusion or by special initial states. 



A. Directed percolation classes 

The directed percolation (DP) introduced by (Broad- 
bent and Hammersley, 1957) is an anisotropic percolation 
with a preferred direction t. This means that this prob- 
lem should be d > 2 dimensional. If there is an object 
(bond, site etc.) at {xi,yj, ...,tk) it must have a nearest 
neighboring object at tk-i unless tk = (sec Figure 3). If 
we consider the preferred direction as the time we recog- 
nize a spreading process of an agent A that can not have 
spontaneous source: (f) -/^ A. This results in the possibil- 
ity of a completely frozen, so called absorbing state from 
which the system cannot escape if it has fallen into it. 
As a consequence these kinds of models may have phase 
transitions in d = 1 spatial dimension already. By in- 
creasing the branching probability p of the agent we can 
have a phase transition between the absorbing state and 
an active steady state with finite density of As. If the 
transition is continuous it is very likely that it belongs 
to a robust universality (DP) class. For a long time all 
examples of such absorbing phase transitions were found 
to belong to the DP class and a conjecture was advanced 
by (Grassberger, 1982a; Grinstein et al, 1989: Janssen, 
1981). This claims that in one component systems ex- 
hibiting continuous phase transitions to a single absorb- 
ing state (without extra symmetry, inhomogeneity or dis- 
order) short ranged interactions can generate DP class 
transition only. This hypothesis has been confirmed by 
all examples up to now, moreover DP class exponents 
were discovered in some systems with multiple absorb- 
ing states. For example in systems with infinitely many 
frozen absorbing states (Jensen, 1993a; Jensen and Dick- 
man, 1993a, b; Mendes et al, 1994) the static exponents 
were found to coincide with those of DP. Furthermore by 



models without any special symmetry of the absorbing 
states DP behavior was reported (Menyhard and Odor, 
1996; Odor and Menyhard, 1998; Park and Park, 1995) 
too. So although the necessary conditions for the DP 
behavior seem to be confirmed the determination of suf- 
ficient conditions is an open problem. There are many 
introductory works available now to DP (Grassberger, 
1996; Hinrichsen, 2000a; Kinzel, 1983; Marro and Dick- 
man, 1999) therefore I shall not go very deeply into the 
discussions of details of various representations. 

In the reaction-diffusion language the DP is built up 
from the following processes 



A0<^0A A^2A 2A^A 



(81) 



The mean-field equation for the coarse-grained particle 
density p{t) is 

dp 



dt 



This has the stationary stable solution 

I tor : (T < 7 



(82) 



(83) 



exhibiting a continuous transition at a = 7. A small 
variation of <t or 7 near the critical point implies a linear 
change of p, therefore the order parameter exponent in 
the mean-field approximation is P — 1. Near the critical 
point the 0{p) term is the dominant one, hence the den- 
sity approaches the stationary value exponentially. For 
a — J the remaining 0(/9^) term causes power-law de- 
cay: p oi t~^ indicating a — 1. To get information about 
other scaling exponents we have to take into account the 
diffusion term DV^ describing local density fluctuations. 
Using rescaling invariance two more independent expo- 
nents can be determined: = 1/2 and Z = 2 if d > 4. 
Therefore the upper critical dimension of directed perco- 
lation is dc = 4. 

Below the critical dimension the RG analysis of the 
Langevin equation 



dpjx, t) 
dt 



= DV^p{x,t) + {a--f)p{x,t)- 
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- {X + a)p^{x,t) + ^/p{x,t)ij{x,t) (84) 

is necessary (Janssen, 1981). Here ri{x, t) is the Gaussian 
noise field, defined by the correlations 



< V{x,t)r]{x^,t') > 



(85) 
rS'^ix- x')S{t-t^) (86) 



The noise term is proportional ^ p(x^ t) ensuring that in 
the absorbing state {p{x,t) = 0) it vanishes. The square- 
root behavior stems from the definition of p{x,t) as a 
coarse-grained density of active sites averaged over some 
mesoscopic box size. Note that DP universality occurs 
in many other processes like in odd offspring, branch- 
ing and annihilating random walks (BARWo) (see Sec- 
tion IV. A. 3) or in models described by field theory with 
higher order terms like p^{x,t) or W'^p{x,t), which are 
irrelevant under the RG transformation. This stochastic 
process can through standard techniques (Janssen, 1976) 
be transformed into a Lagrangian formulation with the 
action 



d'^xdt 



D 



(87) 

where is the density field and tp is the response field 
(appearing in response functions) and the action is in- 
variant under the following time-reversal symmetry 



(j3{x,t) ^ —llj{x,—t) , ll}(x,t) ~(j3{x,—t) 



This symmetry yields (Grassberger and de la Torre, 1979; 
Munoz et ai, 1997) the scaling relations 



/3 = 
4(5 + 277 = dz 



(89) 
(90) 



This field theory was found to be equivalent (Cardy and 
Sugar, 1980) to the Reggeon field theory (Abarbanel 
et ai, 1975; Brower et ai, 1978), which is a model of 
scattering elementary particles at high energies and low- 
momentum transfers. 

Perturbative e = A — d renormalization group analy- 
sis (Bronzan and Dash, 1974; Janssen, 1981) up to two- 
loop order resulted in estimates for the critical expo- 
nents shown in the Table XII. The best results ob- 
tained by approximative techniques for DP like improved 
mean- field (Ben-Naim and Krapivsky, 1994), coherent 
anomaly method (Odor, 1995), Monte Carlo simulations 
(Dickman and Jensen, 1991; Grassberger, 1989a, 1992a; 
Grassberger and de la Torre, 1979), series expansions 
(De'Bell and Essam, 1983; Essam et ai, 1988; Jensen, 
1996a, 1999a; Jensen and Dickman, 1993c, d; Jensen and 
Guttmann, 1995, 1996), DMRG (Carlon et ai, 1999; 
Hieida, 1998), and numerical integration of eq.84 (Dick- 
man, 1994) are also shown in Table XII. 

The local persistence probability may be defined 
as the probability {pi{t)) that a particular site never 
becomes active up to time t. Numerical simulations 
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1/t 



FIG. 4 Local slopes of the density decay in a bosonic BARWl 
model (Sect. IV. A. 3). Different curves correspond to A = 
0.12883, 0.12882, 0.12881 0.1288, 0.12879 (from bottom to 
top) (Odor and Menyhard, 2002). 



(Hinrichsen and Koduvely, 1998) for this in the 1+1 d 
Domany-Kinzel SCA (see Sect. IV. A. 2) found power-law 
with exponent 



Oi = 1.50(1) 



(91) 



The global persistence probability defined here as the 
probability {pg{t)) that the deviation of the global den- 
sity from its mean value does not change its sign up to 
time t. The simulations of (Hinrichsen and Koduvely, 
1998) in 1+1 dimensions claim 9g> 9i. This agrees with 
field theoretical e = 4 — d expansions (Oerding and van 
Wijland, 1998) that predict = 2 for d > 4 and for 
d < 4: 



9.=2-|+0(.^) 



(92) 



The crossover from isotropic to directed perco- 
lation was investigated by perturbative RG (Frey et al, 
1993) up to one loop order. They found that while for 
d > 6 the isotropic, for d > 5 the directed Gaussian 
fixed point is stable. For d < 5 the asymptotic behavior 
is governed by the DP fixed point. On the other hand 
in d = 2 exact calculations and simulations (Frojdh and 
den Nijs, 1997) found that the isotropic percolation is sta- 
ble with respect to DP if we control the crossover by a 
spontaneous particle birth parameter. It is still an open 
question what happens for 2 < d < 5. Crossovers to 
mean-field behavior generated by long-range interactions 
IV. A. 7 and to compact directed percolation IV. C will be 
discussed later. 

It was conjectured (Deloubriere and van Wijland, 
2002) that in Id "fermionic" (single site occupancy) and 
bosonic (multiple site occupancy) models may exhibit 
different critical behavior. An attempt for a fermionic 
field theoretical treatment of the DP in 1+1 d was 
shown in (Brunei et at, 2000; Wijland, 2001). This run 
into severe convergence problems and has not resulted 
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in precise quantitative estimates for the critical expo- 
nents. Although the bosonic field theory is expected to 
describe the fermionic case owing to the asymptotically 
low density at criticality it has never been proven rigor- 
ously. Since only bosonic field theory exists that ^ives 
rather inaccurate critical exponent estimates, in (Odor 
and Menyhard, 2002) simulations of a a BARWl process 
(101) with unrestricted site occupancy were performed to 
investigate the density {p{t)) decay of a DP process from 
random initial state. Figure 4 shows the local slopes of 
density decay defined as 



V , ' , ' 

« • ♦ • • — 



dlogpjt) 
dlogt 



(93) 



around the critical point for several annihilation rates 
(A). The critical point is estimated at Ac = 0.12882(1) 
(corresponding to straight line) with the extrapolated 
decay exponent a = 0.165(5), which agrees well with 
fermionic model simulation and series expansion results 
0.1595(1) (Dc^Bell and Essam, 1983). 

Note that in site restricted models already the 



AID ^ 0A 



2A 



(94) 



processes generate a DP class phase transition, while in 

the bosonic version the 2A-^ A process is also necessary 
to ensure an active steady state (without it the density 
blows up for (T > 7). 

Models exhibiting DP transitions have been reviewed 
in great detail in (Marro and Dickman, 1999) and in (Hin- 
richscn, 2000a) . In the next subsections I recall only three 
important examples. 



1. The Contact process 

The contact process is one of the earliest and sim- 
plest lattice model for DP with asynchronous update 

introduced by Harris (Grassberger and de la Torre, 1979; 
Harris, 1974b) to model epidemic spreading without im- 
munization. Its dynamics is defined by nearest-neighbor 
processes that occur spontaneously due to specific rates 
(rather than probabilities). In numerical simulations 
models of this type are usually realized by random se- 
quential updates. In one dimension this means that a 
pair of sites {si,Si+i} is chosen at random and an up- 
date is attempted according to specific transition rates 
w{si^t+dt, Si+i^t+dt\si^t, Si+i,t)- Each attempt to up- 
date a pair of sites increases the time t hy dt = 1/N, 
where A'^ is the total number of sites. One time step 
(sweep) therefore consists of N such attempts. The con- 
tact process is defined by the rates 

w{A,I\A,A) = w{I,A\A,A) = fi, (95) 
w{I,I\AJ) = w{I,I\I,A) = A, (96) 
w{A,A\A,I) = w{A,A\I,A) = 1, (97) 

where A > and n > arc two parameters (all other 
rates are zero). Equation (95) describes the creation of 



t+1 



FIG. 5 Update in the Domany-Kinzel model 



inactive (/) spots within active (A) islands. Equations 
(96) and (97) describe the shrinkage and growth of active 
islands. In order to fix the time scale, we chose the rate 
in Eq. (97) to be equal to one. The active phase is re- 
stricted to the region A < 1 where active islands are likely 
to grow. In one dimension series expansions and numeri- 
cal simulations determined the critical point and critical 
exponents precisely (De'Bell and Essam, 1983; Dickman 
and Jensen, 1991; Dickman and de Silva, 1998; Essam 
et ai, 1988; Jensen, 1996a; Jensen and Dickman, 1993c, d; 
Jensen and Guttmann, 1995, 1996). In two dimensions 
the order parameter moments and the cumulant ratios 
were determined by (Dickman and de Silva, 1998). 



2. DP-class stochastic cellular automata 

Cellular automata as the simplest systems exhibit- 
ing synchronous dynamics have extensively been stud- 
ied (Wolfram, 1983). By making the update rules prob- 
abilistic phase transitions as the function some control 
parameter may emerge. There are many stochastic cel- 
lular automata (SCA) that exhibit DP transition (Boc- 
cara and Roger, 1993) perhaps the first and simplest one 
is the (l-l-l)-dimensional Domany-Kinzel (DK) model 
(Domany and Kinzel, 1984). In this model the state 
at a given time t is specified by binary variables {si}, 
which can have the values A (active) and / (inactive). 
At odd times odd-indexed, whereas at even times the 
rest of the sites are updated according to specific con- 
ditional probabilities. This defines a cellular automaton 
with parallel updates (discrete time evolution) acting 
on two independent triangular sub-lattices (Fig. 5). The 
conditional probabilities in the Domany-Kinzel model 
P{si,t+i I Si-i,t , Si+i^t) are given by 

P{I\I,I) = 1, (98) 

PiA\A,A) = p2, (99) 
P{A\I,A) = P{A\A,I) = pi, (100) 

and P{I\si-i, Si+i) + P{A\si^i, Si+i) = 1, where < 
Pi < 1 and < P2 < 1 are two parameters. Equation 
(98) ensures that the configuration ...,/,/,/,... is the 
absorbing state. The process in Eq. (99) describes the 
creation of inactive spots within active islands with prob- 
ability 1 — p2- The random walk of boundaries between 
active and inactive domains is realized by the processes 
in Eq. (100). A DP transitions can be observed only if 
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where X G {0,1} is a two valued random variable sueli 
that Prob{X = I) = p (Odor and Szolnoki, 1996). The 
y = 3 case was introduced and investigated by (Bidaux 
et ai, 1989) in d = 1,2,3. The very first simulations in 
one dimension (Bidaux et al., 1989) suggested a counter- 
example to the DP conjecture. More precise spreading 
simulations of this model (Jensen, 1991), GMF + CAM 
calculations and simulations of the y < Q family in one 
and two dimensions have proven that this does not hap- 
pen for any case (Odor and Szolnoki, 1996). The transi- 
tions are either belong to DP class or first order. 



FIG. 6 Phase diagram of the Id Domany-Kinzel SCA (Hin- 
richsen, 1997). 



Pi > 5, when active islands are biased to grow (Wol- 
fram, 1983). The phase diagram of the Id DK model is 
shown in Fig. 6. It comprises an active and an inactive 
phase, separated by a phase transition line (solid line) 
belonging to DP class. The dashed lines corresponds to 
directed bond percolation (p2 = p{2 — pi)) and directed 
site percolation {pi = P2) models. At the special symme- 
try endpoint [pi = \, P2 = 1) compact domain growth 
occurs (CDP) and the transition becomes first order (see 
Section IV. C). The transition on the P2 = axis cor- 
responds to the transition of the stochastic version of 
Wolfram's rule-18 cellular automaton (Wolfram, 1983). 
This range-1 SCA generates A at time t only when the 
right or left neighbor was A at f — 1: 



t-1: 
t: 



All 
A 



IIA 
A 



with probability pi (Boccara and Roger, 1993). The crit- 
ical point was determined by precise simulations {p\ = 
0.80948(1)) (Odor and al., 1997). In the f ^ 00 limit the 
steady state is built up from II and lA blocks (Eloranta 
and Nummelin, 1992). This finding permits us to map 
this model onto an even simpler one, the rule-6/16 SCA 
with new variables: lA — > A and II — > I: 



t-1: 
t: 



I I 
I 



I A 
A 



A I 
A 



A A 
I 



By solving GMF approximations and applying Padc ap- 
proximations (Szabo and Odor, 1994) or CAM method 
(Odor, 1995) very precise order parameter exponent esti- 
mates were found: (3 = 0.2796(2). The DS phase struc- 
ture of the Id DK model was explored by (Hinrichscn 
et al., 1997) and DP class transitions were found. 

Another SCA example I mention is the family of 
range-4 SCA with an acceptance rule 



_ \ X if y<j:A^s{t,j)<6 



s{t + l,j) = 



3. Branching and annihilating random walks with 
odd number of offspring 

Branching and annihilating random walks (BARW) in- 
troduced by (Takayasu and Tretyakov, 1992) can be re- 
garded as generalizations of the DP process. They are 
defined by the following reaction-diffusion processes: 



A-^{m + 1)A kA 



A% & 9 A 



(101) 



The 2A — > A and 2A reactions dominating in the in- 
active phase are shown to be equivalent (Peliti, 1986) (see 
Sect.IV.C.l). Therefore the fc = 2 and to = 1 (BARWl) 
model differs from the DP process (81) that spontaneous 
annihilation of particles is not allowed. The action of the 
BARW process was set up by (Cardy and Tauber, 1996, 
1998) 



+u{l - V''")V'?5'] 



(102) 



otherwise 



The bosonic RG analysis of BARW systems (Cardy and 
Tauber, 1996) proved that for fc = 2 all the lower branch- 
ing reactions with to — 2, to — 4, ... are generated via fluc- 
tuations involving combinations of branching and annihi- 
lation processes. As a consequence for odd to the A ^ % 
reaction appears (via A 2 A 0). Therefore after the 
first coarse graining step in the BARWl (and in general 
in the odd m BARW (BARWo) cases) the action becomes 
the the same as that of the DP process. The fluctuations 
are relevant for d < 2 and the universal behavior is DP 
type, with dc = 2. This prediction was conflrmed by 
simulations (see for example (Jensen, 1993b)). 

For even m (B ARWe) , when the parity of the number 
of particles is conserved the spontaneous decay A — > is 
not generated, hence there is an absorbing state with a 
lonely wandering particle. This systems exhibits a non- 
DP class critical transition, which will be discussed in 
Sect.IV.D.l. 



4. DP with spatial boundary conditions 

For a review of critical behavior at surfaces of equi- 
librium models see (Igloi et al., 1993)). Ref. (Cardy, 
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1983a) suggested that surface critical phenomena may be 
described by introducing an additional surface exponent 
for the order parameter field which is generally indepen- 
dent of the other bulk exponents. In nonequilibrium sta- 
tistical physics one can introduce spatial, temporal (see 
Sect.IV.A.7) or mixed (see Sect.IV.A.5) boundary condi- 
tions. 

In DP an absorbing wall may be introduced by cut- 
ting all bonds (IBC) crossing a given ((i-l)-dimensional 
hyperplane in space. In case of reflecting boundary con- 
dition (RBC) where the wall acts like a mirror so that 
the sites within the wall are always a mirror image of 
those next to the wall. A third type of boundary condi- 
tion is the active boundary condition (ABC) where the 
sites within the wall are forced to be active. The density 
at the wall is found to scale as 



(103) 



with a surface critical exponent (3i > [3. Owing to the 
time reversal symmetry of DP (88) only one extra ex- 
ponent is needed to describe surface effects, hence the 
cluster survival exponent is the same 



(104) 



The mean lifetime of finite clusters at the wall is defined 
as 



where Ag 
relation 



{t) ~ lAI-""! (105) 
{p — Pc), and is related to /3i by the scaling 



n = I'll ~ 



(106) 



The average size of finite clusters grown from seeds on 
the wall is 



(107) 



Series expansions (Essam et al, 1996) and numerical sim- 
ulations (Lauritsen et al., 1997) in 1-1-1 dimensions indi- 
cate that the presence of the wall alters several expo- 
nents. However, the scaling properties of the correlation 
lengths (as given by and z/_l) are not altered. 

The field theory for DP in a semi-infinite geometry was 
first analyzed by (Janssen et al, 1988). They showed that 
the appropriate action for DP with a wall at a;_L = is 
given hy S = S'buik + ^'surface, where 



^surface 



dt As Vs 



(108) 



with the definitions (/)s = (/)(x||,xx = 0,f) and Vs = 

'!/'(x||,.Tx = 0, i). The surface term S'suifacc corresponds 
to the most relevant interaction consistent with the sym- 
metries of the problem and which also respects the ab- 
sorbing state condition. The appropriate surface expo- 
nents were computed to first order in e = 4 — cZ using 
renormalization group techniques: 

/3i = |- J + 0(e2)- (109) 



They also showed that the corresponding hyperscaling 
relation is 



relating (3\ to 



Vll + dl/± = Pi+(3 + 'yi 



7.4 + S + 0(e^). 



(110) 



(111) 



The schematic phase diagram for boundary DP is shown 
in Fig. 7, where A and Ag represent, respectively, the 
deviations of the bulk and surface from criticality. For 
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FIG. 7 Schematic mean field phase diagram for boundary 
DP (Frojdh et al, 2001). The transitions are labeled by 
0=ordinary, E=extraordinary, S=surface, and Sp=special. 

A > and for As sufficiently negative the boundary 
orders even while the bulk is disordered, there is a surface 
transition. For Ag < and A — > 0, the bulk orders in 
the presence of an already ordered boundary, there is 
an extraordinary transition of the boundary. Finally at 
A = Ag = 0, where all the critical lines meet, and where 
both the bulk and isolated surface are critical, we find a 
multi-critical point, i.e. the special transition. 

For Ag > and A ^ there is an ordinary transi- 
tion, since the bulk orders in a situation in which the 
boundary, if isolated, would be disordered. At the or- 
dinary transition, one finds just one extra indepen- 
dent exponent related to the boundary: this can be 
the surface density exponent /3i,dens- In Id the IBC and 
RBC cases belong to the same universality class (Lau- 
ritsen et al., 1998) that was identified as the ordinary 
transition. There are numerical data for the exponents 
of the extraordinary and special transitions (however see 
(Janssen et al., 1988) for an RG analysis). 

The best exponent estimates currently available were 
summarized in (Frojdh et al, 2001). Some of them are 
shown in Table XIII. In d = 1 the best results are from 
series expansions (Essam et al., 1996; Jensen, 1999b); in 
all other cases are from Monte-Carlo data. (Frojdh et al, 
1998; Howard et al, 2000; Lauritsen et al, 1998, 1997) 
The exponent ti was conjectured to equal unity, (Essam 
et al, 1996) although this has now been challenged by 
the estimate n = 1.00014(2) (Jensen, 1999b). 

It has been known for some time that the presence 
of an edge introduces new exponents, independent of 
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rf= 1 


d = 2 


Mean Field 


01 


0.733 71(2) 


1.07(5) 


3/2 


Si = ai 


0.423 17(2) 


0.82(4) 


3/2 


Tl 


1.000 14(2) 


0.26(2) 





71 


1.820 51(1) 


1.05(2) 


1/2 



TABLE XIII Critical exponents for DP in d 
for the ordinary transition at the boundary. 



1 and d 





those associated with the bulk or with a surface (see 
(Cardy, 1983a)). For an investigation showing numer- 
ical estimates in 2d and mean-field values sec (Frojdh 
et ai, 1998). Table IV. A. 4 summarizes results for the 
ordinary edge exponents. A closely related application is 
the study of spreading processes in narrow channels (Al- 
bano, 1997). 



Angle (a) 


7r/2 


37r/4 


TT 


57r/4 


P? (d = 2) 


1.6(1) 


1.23(7) 


1.07(5) 


0.98(5) 


(MF) 


2 


5/3 


3/2 


7/5 



TABLE XIV Numerical estimates for the ordinary 02 expo- 



nents for edge DP together with the mean field values, 
that /??(7r) =/3P 



Note 



5. DP with mixed (parabolic) boundary condition 
scaling 

Boundary conditions, which act both in space and time 
direction can also been investigated in dynamical sys- 
tems. These turn out to be related to hard-core repulsion 
effects of Id systems (Sect.V.B). (Kaiser and Turban, 
1994, 1995) investigated the 1-1-1 d DP process confined 
in a parabola-shaped geometry. Assuming an absorbing 
boundary of the form x = ±Ct'^ they proposed a general 
scaling theory. It is based on the observation that the 
coefficient of the parabola (C) scales as C ^ A^'^~^C 
under rescaling 



Aa 



A 



-p, 

(112) 

where A = |p — Pc| and Z is the dynamical exponent of 
DP. By referring to conformal mapping of the parabola 
to straight lines and deriving it in the mean-field approx- 
imation they claimed that this boundary is a relevant 
perturbation for u > 1/Z, irrelevant for cr < 1/Z and 
marginal for a = l/Z (sec Fig. 8). The marginal case 
results in C dependent non-universal power-law decay, 
while for the relevant case stretched exponential func- 
tions have been obtained. The above authors have given 
support to these claims by numerical simulations. 




FIG. 8 The space-time evolution of the critical, 1+1 dimen- 
sional directed site percolation process confined by parabola 
(Kaiser and Turban, 1994). (a) a < 1/Z, (b) a = 1/Z, (c) 
a > 1/Z. 



6. Levy flight anomalous diffusion in DP 

Levy flight anomalous diffusion generating long-range 
correlations was already mentioned in Sections III. A. 5 
and III.B.3 in case of equilibrium models. In non- 
equilibrium systems following the suggestion of (MoUi- 
son, 1977) ref. (Grassbcrgcr, 1986) introduced a varia- 
tion of the epidemic processes with infection probability 
distribution P{R), which decays with the distance R as 
a power-law like 



P{R) oc 



1 



(113) 



This can model long-range epidemics mediated by flies, 

wind, ... etc. He claimed that the critical exponents 
should depend continuously on a. This result was con- 
firmed in Id by estimating (3 based on CAM calcula- 
tions in (Marques and Ferreira, 1994). The study of 
anomalous diffusion was extended for GEP processes (see 
Sect. IV. B) and for annihilating random walks (see Sect. 
IV.C.4) too. The effective action of the Levy flight DP 
model is 



I 



df-x dt 



(114) 



where Djv denotes the normal. Da the anomalous dif- 
fusion constant and g is the interaction coupling con- 
stant. Field theoretical RG method up to first order in 
e = 2a — d expansion (Janssen et ai, 1999) gives: 
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FIG. 9 Estimates for the exponent (5 and the derived expo- 
nents and 1 in comparison with the fiold-thcorctic results 
(soUd hnos) and the DP exponents (dot-dashed hnes) (Hin- 
richsen and Howard, 1999). The quantities Ai and A2 rep- 
resent deviations from the scahng relations (116) and (117), 
respectively. 

v\\ = 1 + ^+0 (e^) , (115) 
Z=^ = a-e/7 + Oie^) . 

I/_L 

Moreover, it was shown that the hyperscaling relation 

r] + 26 = d/Z (5 = /3/i/||), (116) 

for the so-called critical initial slip exponent rj and the 
relation 

u\i - ujl{(t - d) - 2f3 = . (117) 

hold exactly for arbitrary values of a. Numerical simula- 
tions on H-l d bond percolation confirmed these results 
except in the neighborhood of a = 2 (Hinrichsen and 
Howard, 1999) (see Fig.9). 

7. Long-range correlated initial conditions in DP 

It is well known that initial conditions influence the 
temporal evolution of nonequilibrium systems. The 
"memory" of systems for the initial state usually de- 
pends on the dynamical rules. For example, stochastic 
processes with a finite temporal correlation length relax 
to their stationary state in an exponentially short time. 
An interesting situation emerges when a system under- 
goes a nonequilibrium phase transition in which the tem- 
poral correlation length diverges. This setup motivates 
the question whether it is possible construct initial states 
that affect the entire temporal evolution of such systems. 

Montc-Carlo simulations of critical models with ab- 
sorbing states usually employ two different types of initial 
conditions. On the one hand uncorrelated random initial 
conditions (Poisson distributions) arc used to study the 
relaxation of an initial state with a finite particle density 



towards the absorbing state. In this case the particle den- 
sity p{t) decreases on the infinite lattice asymptotically 
as 

p{t) ~ t-^'/'^ii . (118) 

On the other hand, in spreading simulations (Grass- 
berger and de la Torre, 1979), each run starts with a 
single particle localized active seed from where a 
cluster originates (this is a long-range correlated state). 
Although many of these clusters survive for only a short 
time, the number of particles n{t) averaged over many 
independent runs increases as 

{n{t)} ~ t+^ , (119) 

These two cases seem to represent extremal situations 
where the average particle number either decreases or 
increases. 

A crossover between these two extremal cases takes 
place in a critical spreading process that starts from a 
random initial condition of very low density. Here the 
particles are initially separated by empty intervals of a 
certain typical size, wherefore the average particle num- 
ber first increases according to Eq. (119). Later, when 
the growing clusters begin to interact with each other, the 
system crosses over to the algebraic decay of Eq. (118) - 
a phenomenon which is referred to as the "critical initial 
slip" of nonequilibrium systems (Janssen et al, 1989). 

In (Hinrichsen and Odor, 1998: Mcnyhard and Odor, 
2000) it was investigated whether it is possible to inter- 
polate continuously between the two extremal cases in 
cases of 1+1 dimensional DP and PC processes. It was 
shown that one can in fact generate certain initial states 
in a way that the particle density on the infinite lattice 
varies as 

p{t) ~ i« (120) 
with a continuously adjustable exponent k in the range 

-l3/iy\\ <K<+r]. (121) 

To this end artificial initial configurations with algebraic 
long-range correlations of the form 

C{r)^{s,s,+r)^r-^''-^\ (122) 

were constructed, where () denotes the average over many 
independent realizations, d the spatial dimension, and 
Si = 0, 1 inactive and active sites. The exponent cr is a 
free parameter and can be varied continuously between 
and 1. This initial condition can be taken into account 
by adding the term 

Sic = II j d''xi;ix,0)Mx) (123) 

to the action, where (pai^x) represents the initial particle 
distribution. The long-range correlations limit a ^ d 
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corresponds to a constant particle density and thus we 
expect Eq. (118) to hold {(poix) — const, is irrelevant 
under rescaling). On the other hand, the short-range 
limit (7 — » represents an initial state where active sites 
are separated by infinitely large intervals {4>o{x) = S'^{x)) 
so that the particle density should increase according to 
Eq. (119). In between we expect p{t) to vary algebraically 
according to Eq. (120) with an exponent k depending 
continuously on a. 

In case of the 1+1 d Domany-Kinzel SCA (see 
Sect. IV. A. 2) field-theoretical renormalization group cal- 
culation and simulations proved (Hinrichsen and Odor, 
1998) the exact functional dependence 

{7] for a < CTc 

(124) 
^{d — a - for ct > CTc 

with the critical threshold ac = P/u^. 



reversal symmetry of the DP process (88) is not broken 
therefore 

5 = 5^ (126) 
and derived a hyper-scaling law for the inactive phase 

r] = dz/2 (127) 
and for the absorbing phase 

r] + 6 = dz/2 . (128) 

They confirmed these relations by simulations and found 

that the dynamical exponents change continuously as the 
function of the disorder probability. An RG study by 
(Hooyberghs et ai, 2002) showed that in case of strong 
enough disorder the critical behavior is controlled by an 
infinite randomness fixed point (IRFP), the static expo- 
nents of which in Id are 



8. Quench disordered DP systems 



/3 = (3 - V5)/2 , i^x = 2 



(129) 



Perhaps the lack of experimental observation of the 
robust DP class lies in the fact that even weak disorder 
changes the critical behavior of such models. Therefore 
this section provides a view onto interesting generaliza- 
tions of DP processes that may be observable in physical 
systems. First (Noest, 1986, 1988) showed using Harris 
criterion (Harris, 1974a) that spatially quenched dis- 
order (frozen in space) changes the critical behavior of 
DP systems for d < 4. Ref. (Janssen, 1997a) studied the 
problem by field theory taking into account the disorder 
in the action by adding the term 



5 + 7 



dtijxj) 



(125) 



This additional term causes marginal perturbation 

and the stable fixed point is shifted to an unphysical re- 
gion, leading to runaway solutions of the flow equations in 
the physical region of interest. This means that spatially 
quenched disorder changes the critical behavior of DP. 
This conclusion is supported by the simulation results of 
(Moreira and Dickman, 1996) who reported logarithmic 
spreading behavior in two-dimensional contact process 
at criticality. In the sub-critical region they found Grif- 
fiths phase in which the time dependence is governed by 
non- universal power-laws, while in the active phase the 
relaxation of P{t) is algebraic. 

In 1+1 dimension (Noest, 1986, 1988) predicted 
generic scale invariance. Ref. (Webman et al., 1998) 
reported glassy phase with non-universal exponents in a 
1+1 d DP process with quenched disorder. Ref. (Cafiero 
et al., 1998) showed that DP with spatially quenched ran- 
domness in the large time limit can be mapped onto a 
non-Markovian spreading process with memory, in agree- 
ment with previous results. They showed that the time 



and ^^/^ oc Inr. For disorder strengths outside the at- 
tractive region of the IRFP disorder dependent critical 
exponents are detected. 

The temporally quenched disorder can be taken 
into the action by adding the term: 



(130) 



This is a relevant perturbation for the DP processes. 
Ref. (Jensen, 1996b) investigated the 1+1 d directed 
bond percolation (see Sect. IV. A. 2) with temporal dis- 
order via series expansions and Monte Carlo simulations. 
The temporal disorder was introduced by allowing time 
slices to become fully deterministic (pi = p2 = 1), 
with probability a. He found a dependent, continuously 
changing critical point and critical exponent values be- 
tween those of the the 1+1 d DP class and those of the 
deterministic percolation. This latter class is defined by 
the exponents: 

/3 = 0, ^ = 0, = 1, Z=l, I'll = 2, z/_L = 2 . 

(131) 

For small disorder parameter values violation of the Har- 
ris criterion is reported. 

If quenched disorder takes place in both space 
and time the corresponding term to action is 

5^5 + 7 j dtd'^x [il)(l>f . (132) 

and becomes an irrelevant perturbation to the Reggeon 
field theory. This has the same properties as the intrin- 
sic noise in the system and can be considered as being 
annealed. 
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B. Dynamical percolation (DyP) classes 

If we allow memory in the unary DP spreading pro- 
cess (Sect. IV. A) such that the infected sites may have a 
different re-infection probability (p) than the virgin ones 
(q) wc obtain different percolation behavior (Grassberger 
et ai, 1997). The model in which the re-infection proba- 
bility is zero is called the General Epidemic Model (GEP) 
(MoUison, 1977). In this case the epidemic stops in fi- 
nite systems but an infinite epidemic is possible in the 
form of a solitary wave of activity. When starting from 
a single seed this leads to annular growth patterns. The 
transition between survival and extinction is a critical 
phenomenon called dynamical percolation (Grassberger, 
1982b). Clusters generated at criticality are the ordinary 
percolation clusters of the lattice in question. Field theo- 
retical treatment was given by (Cardy, 1983b; Cardy and 
Grassberger, 1985; J.-Dalmaroni and Hinrichsen, 2003; 
Janssen, 1985; Mufioz et al., 1998). The action of the 
model is 



S= I d''xdt^^l)'^(t>-^(^t(t>-^'^(l>-r(t> + w(t> ds(j){s) 

(133) 

This is invariant under the non-local symmetry transfor- 
mation 
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FIG. 10 Isotropic site percolation in d = 2 dimensions 
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TABLE XV Critical exponents of the ordinary percolation 



(l){x,t) ^ -dttp{x,-t) , (134) 

which results in the hyperscaling relation (Mufioz et al., 
1997): 

rj + 2S+l = ^ (135) 

As in case of the DP, the relations (3 = f3' and S = a 
again hold. The upper critical dimension is rfc = 6. The 
dynamical critical exponents as well as spreading and 
avalanche exponents are summarized in (Munoz et al., 
1999). The dynamical exponents are Z = 1.1295 for 
d = 2, Z = 1.336 for d = 3 and Z = 2 for d = 6. 
Dynamical percolation was observed in forest fire models 
(Albano, 1994; Drossel and Schwabl, 1993) and in some 
Lotka-Volterra type lattice prey-predator models (Antal 
et al, 2001) as well. 

1. Isotropic percolation universality classes 

The ordinary percolation (Grimmett, 1999; Stauffer 
and Aharony, 1994) is a geometrical phenomenon that 
describes the occurrence of infinitely large connected 
clusters by completely random displacement of some 
variables (sites, bonds, etc) (with probability p) on lat- 
tices (sec Fig. 10). 

The dynamical percolation process is known to gen- 
erate such percolating clusters (see Sect. IV. B). At the 
transition point moments of the s cluster size distribu- 
tion ns{p) show singular behavior. The ordinary perco- 
lation corresponds to the q = I limit of the Potts model. 



That means its generating functions can be expressed in 
terms of the free energy of the 9 — > 1 Potts model. In the 
low-temperature dilute Ising model the occupation prob- 
ability (p) driven magnetization transition is an ordinary 
percolation transition as well. As a consequence the crit- 
ical exponents of the magnetization can be related to 
the cluster-size exponents. For example the susceptibil- 
ity obeys simple homogeneity form with p — Pc replacing 

XocIp-PcI"^ (136) 

Table XV summarizes the known critical exponents of 
the ordinary percolation. The exponents are from the 
overview (Bunde and Havlin, 1991). Field theoretical 
treatment (Bcnzoni and Cardy. 1984) provided an upper 
critical dimension dc = Q. The d = 1 case is special : Pc = 
1 and the order parameter jumps (/3 = 0). Furthermore 
here some exponents exhibit non-universal behavior by 
increasing the interaction length unless the we redefine 
the scaling variable (see (Stauffer and Aharony, 1994)). 

2. DyP with spatial boundary conditions 

There are very few numerical results exists for sur- 
face critical exponents of the dynamical percolation. 
The GEP process in 3d was investigated numerically by 
Grassberger (Grassberger, 1992b). The surface and edge 
exponents (for angle 7r/2) were determined in case of IBC. 
Different measurements (density and cluster simulations) 
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resulted in a single surface /3i = 0.848(6) and a single 
edge /3e = 1.36(1) exponent. 



3. Levy flight anomalous diffusion in DyP 

To model long-range epidemic spreading in a system 
with immunization the effect of Levy flight diffusion (113) 
was investigated by (.Jansscn et ai, 1999). The rcnor- 
malization group analysis of the GEP with anomalous 
diffusion resulted in the following e = 3a — d expansion 
results for the critical initial slip exponent: 



(137) 



for the order parameter (density of removed (immune) 
individuals) exponent: 



for the spatial correlation exponent 

and for the temporal correlation length exponent 



-ll = l + I^+0(e^ 



(138) 



(139) 



(140) 



C. Voter model (VM) classes 

Now we turn to models which can describe the spread- 
ing of opinion of voters arranged on regular lattices. They 
exhibit first order transition still dynamical scaling can 
be observed in them. The voter model (Durret, 1988; 
Liggct, 1985) is defined by the following spin-flip dynam- 
ics. A site is selected randomly which takes the "opinion" 
(or spin) of one of its nearest neighbors (with probabil- 
ity p). This rule ensures that the model has two ho- 
mogeneous absorbing states (all spin up or down) and 
is invariant under the Zi symmetry. General feature of 
these models is that dynamics takes place only at the 
boundaries. The action, which describes this behavior 
was proposed in (Dickman and Tretyakov, 1995; Peliti, 
1986) 



:^</.(l-</.)V''-V(c^t0-AV^ 



(141) 



is invariant under the symmetry transformation: 

^ 1 - (/), tp ^ -tp . (142) 

This results in the "hyperscaling" relation(Muhoz et ai, 
1997) 



that is valid for all first order transitions (/? = 0) with 
d < 2, hence dc = 2 is the upper critical dimension. It 
is also valid for all compact growth processes (where 
"compact" means that the density in surviving colonies 
remains finite as t oo). 

In one dimension at the upper terminal point of the 
DK SCA (Fig. 6 pi = \, P2 = 1) an extra Z2 symmetry 
exists between 1-s and 0-s, hence the scaling behavior is 
not DP class like but corresponds to the fixed point of 
the inactive phase of PC class models (Sect. IV. D). As 
a consequence compact domains of 0-s and 1-s grow such 
as the domain walls follow annihilating random walks 
(ARW) (see Scct.IV.C.l) and belong to the Id VM class. 
In Id the compact directed percolation (CDP) is also 
equivalent to the T = Glauber Ising model (see Sec- 
tions II,III.A). By applying non-zero temperature (corre- 
sponding to spin flips in domains) or symmetry breaking 
(like changing p2 or adding an external magnetic field) a 
first order transition takes place (/3 = 0) . 

In two (and higher) dimensions the p = 1 situa- 
tion corresponds to the pi = 3/4, p^ = 1 point in the 
phase diagram of symmetric models (see Fig.l). This 
model has a "duality" symmetry with coalescing random 
walks: going backward in time, the successive ancestors 
of a given spin follow the trail of a simple random walk; 
comparing the values of several spins shows that their as- 
sociated random walks necessarily merge upon encounter 
(Ligget, 1985). This correspondence permits us to solve 
many aspects of the kinetics. In particular, the calcula- 
tion of the density of interfaces pm{t) (i.e. the fraction 
of +— nearest neighbour (n.n.) pairs) starting from ran- 
dom initial conditions of magnetization m, is ultimately 
given by the probability that a random walk initially at 
unit distance from the origin, has not yet reached it at 
time t. Therefore, owing to the recurrence properties of 
random walks, the VM shows coarsening for d < 2 (i.e. 
Pm{t) when t 00). For the the 'marginal' d = 2 
case one finds the slow logarithmic decay (Frachebourg 
and Krapivsky, 1996; Krapivsky, 1992; Scheucher and 
Spohn, 1988): 



Pm{t) = (1 - m^) 
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(144) 



6 + r] = dz/2 



(143) 



with D is being the diffusion constant of the underlying 
random walk {D = 1/4 for the standard case of n.n., 
square lattice walks, when each spin is updated on aver- 
age once per unit of time). 

Simulating general, Z2 symmetric spin-flip rules in 2d 
(Dornic et al., 2001) conjectured that all critical Z2- 
symmetric rules without bulk noise form a co-dimension- 
1 'voter-like' manifold separating order from disorder, 
characterized by the logarithmic decay of both p and m. 
The critical exponents for this class are summarized in 
Table XVI. Furthermore ref. (Dornic et ai, 2001) found 
that this Z2 symmetry is not a necessary condition, the 
VM behavior can also be observed in systems without 
bulk fluctuations, where the total magnetization is con- 
served. Field theoretical understanding of these results 
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TABLE XVI Critical exponents of VM classes 



are still lacking. 



1. The 2A ^ (ARW) and the 2A ^ A models 

As it was mentioned in Sect. IV. C in one dimension 
the annihilating random walk and the voter model are 
equivalent. In higher dimensions this is not the case (see 
Sect. IV. E). The simplest reaction-diffusion model - in 
which identical particles follow random walk and anni- 
hilate on contact of a pair - is adequately described by 
mean- field- type equations in dc > 2 dimensions 



pit) (X t 



(145) 



but in lower dimensions fluctuations become relevant. 
Omitting boundary and initial condition terms, the field 
theoretical action is 

[ d'^xdt [ipidt(l) ~ DV^ct)) - A(l - (146) 



where, D denotes the diffusion coefficient and A is the 
annihilation rate. 

For d — dc = 2 the leading order decay of the ARW was 
derived exactly by Lee using field theoretical RG method 
(Lee, 1994): 
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ln(t)/t + 0(l/t) 



(147) 



For d = 1 (Lushnikov, 1987; Racz, 1985) predicted that 
the particle density decays as 



Pit) ^ A2{Dt)-'^' 



(148) 



This scaling law was confirmed by e expansion and the 
universal amplitude A2 was found to be 
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(149) 



The universal scaling behavior of the ARW was shown 
to be equivalent to that of the A + A A coagulation 
random walk process by (Peliti, 1986). The renormaliza- 
tion group approach provided universal decay amplitudes 
(different from those of the ARW) to all orders in ep- 
silon expansion. It was also shown (Domany and Kinzel, 
1984) that the motion of kinks in the compact version 
of directed percolation (CDP) (Essam, 1989) and the 
Glauber-Ising model (Glauber, 1963) at the T = tran- 
sition point are also described exactly by (148). These 
reactions have also intimate relationship to the EW in- 
terface growth model (see Sect. VLB). 




10 

t (MCS) 



FIG. 11 Parabola boundary confinement cluster simulations 
for CDP (Odor and Menyhard, 2000). Middle curves: number 
of active sites (C = 2, 1.5, 1.2, 1 top to bottom); Lower curves: 
survival probabifity (C = 2, 1.5, 1.2, 1 top to bottom); Upper 
curves: R'it) (C = 2, 1.5, 1.2, 1 top to bottom). 



2. Compact DP (CDP) with spatial boundary 
conditions 

By introducing a wall in CDP, the survival probability 
is altered and one obtains surface critical exponents just 
as for DP. With IBC, the cluster is free to approach and 
leave the wall, but not cross. For d = 1, this gives rise 
to (3[ = 2. On the other hand, for ABC, the cluster 
is stuck to the wall and therefore described by a single 
random walker for d = 1. By reflection in the wall, this 
may be viewed as symmetric compact DP which has the 
same /3' as normal compact DP, giving (3[ = 1 (Essam 
and Guttmann, 1995; Essam and TanlaKishani, 1994). 



3. CDP with parabolic boundary conditions 

Space-time boundaries are also of interest in CDP. 
Cluster simulations in 1-1-1 d and MF approximations 
(Dickman and ben Avraham, 2001; Odor and Menyhard, 
2000) for CDP confined by repulsive parabolic boundary 
condition of the form x = ztCf^ resulted in C depen- 
dent 6 and 77 exponents (see Fig.ll) similarly to the DP 
case (see Sect. IV. A. 5) in case of marginal condition: 
a — 1/2. In the mean-field approximations (Odor and 
Menyhard, 2000) similar results were obtained as for the 
DP (Kaiser and Turban, 1995). Analytical results can 
be obtained only in limiting cases. For narrow systems 
(small C) one obtains the following asymptotic behavior 
for the connectedness function to the origin: 



\2C7t) 



(150) 



Recently an analytical solution was derived for a re- 
lated problem (Dickman and ben Avraham, 2001). For a 
one-dimensional lattice random walk with an absorbing 
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boundary at the origin and a movable partial reflector 
(with probability r) S varies continuously between 1/2 
and 1 as r varies between and 1. 



4. Levy flight anomalous diffusion in ARW-s 

Long-range interactions generated by non-local diffu- 
sion in annihilating random walks results in the recovery 
of mean-field behavior. This has been studied by dif- 
ferent approaches. Particles performing simple random 
walks subject to the reactions A + B ^ (Sect.IV.C.l) 
and A + A ^ (Sect. V. A) in the presence of a quenched 
velocity field were investigated in (Zumofcn and Klafter, 
1994). The quenched velocity field enhances the diffusion 
in such a way that the effective action of the velocity field 
is reproduced if Levy fiights are substituted for the simple 
random walk motion. In the above mentioned reactions 
the particle density decay is algebraic with an exponent 
related to the step length distribution of the Levy flights 
defined in Eq. (113). These results have been confirmed 
by several renormalization group calculations (Deem and 
Park, 1998a,b; Oerding, 1996). 

The A + A ^ (f) process with anomalous diffusion 
was investigated by field theory (Hinrichsen and Howard, 
1999). The action of this model is 

Slip, ^] = j d'^x dt ^^tPidt - DnV^ - DaV)^ 



where no is the initial (homogeneous) density at t 
The density decays for a < 2 as : 



0. 



n{t) 



t-'^/'' for d < cr , 
hit ior d= dc = (7 
for d > cr . 



(152) 



The simulation results of the corresponding 1+1 d lat- 
tice model (Hinrichsen and Howard, 1999) can be seen in 
Fig. 12. It was also shown in (Hinrichsen and Howard, 
1999) that Levy flight annihilation and coagulation 
processes {A + A ^ A) are in the same universality 
class. 



• direct measurement 
O extrapolated exponents 
— field-theoretical result 



e 




FIG. 12 The anomalous annihilation process: the graph from 
(Hinrichsen and Howard, 1999) shows direct estimates and 
extrapolations for the decay exponent a, as a function of a. 
The solid line represents the exact result (neglecting log cor- 
rections at cr = 1). 



(see Sect.V.F) in one (Park et ai, 2001) and two dimen- 
sions (Odor et al., 2002) were also found to be insensi- 
tive to the presence of the parity conservation. Multi- 
component BARW2 models in Id (see Sect.V.K) gener- 
ate different, robust classes again (Cardy and Tauber, 
1996; Hooyberghs et ai, 2001; Odor, 2001b). By this 
time it is known that the BARW2 dynamics in single- 
component, single-absorbing state systems (without 
inhomogeneities, long-range interactions and other sym- 
metries) provides sufficient condition for the PC class 
(Cardy and Tauber, 1996). In single-component, 
multi-absorbing state systems the Z2 symmetry en- 
sures necessary but not sufficient condition, the BARW2 
dynamics (Sect.IV.D.l) of domain walls is also a neces- 
sary condition. Some studies have shown (Hwang et al, 
1998; Menyhard and Odor, 1996; Park and Park, 1995) 
that an external field that destroys the Z2 symmetry of 
absorbing states (but preserves the the BARW2 dynam- 
ics) yields a DP instead of a PC class transition in the 
system. Some other names for this class are also used, 
like directed Ising (DI) class, or BARW class. 



D. Parity conserving (PC) classes 

In an attempt to generalize DP and CDP like systems 
we arrive to new models, in which a conservation law is 
relevant. A new universality class appears among 1-1-1 di- 
mensional, single component, reaction-diffusion models. 
Although it is usually named parity conserving class (PC) 
examples have proved that the parity conservation itself 
is not a sufficient condition for the PC class behavior. 
For example in (Inui et ai, 1995) a Id stochastic cellular 
automaton with a global parity conservation was shown 
to exhibit DP class transition. Binary spreading process 



1. Branching and annihilating random walks with 
even number of offspring 

The generic models of the PC class, for which field the- 
oretical treatment exists are branching and annihilating 
random walks - introduced in Sect. IV. A. 3 with k ^ 2 
and even number (m) of offspring (BARWe) (ben Avra- 
ham et al., 1994; Janssen, 1997b; Jensen, 1993c, 1994; 
Lipowski, 1996; Zhong and Avraham, 1995). These con- 
serve the particle number mod 2, hence there are two 
distinct sectors in these models, an odd and an even par- 
ity one. In the even sector particles finally die out (S 7^ 0, 
ry = 0), while in the odd one at least one particle always 
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remain alive (5 = 0, 77 ^ 0). BARWe dynamics may also 
appear in multi-component models exhibiting Z2 sym- 
metric absorbing states in terms of the kinks between 
ordered domains. Such systems are the NEKIM (Sect. 
IV.D.2) and GDK models (Sect. IV.D.3) for example. 
It was conjectured (Menyhard and Odor, 2000) that in 
all models with Z2 symmetric absorbing states an under- 
lying BARWe process is a necessary condition for phase 
transitions with PC criticality. Sometimes it is not so 
easy to find the underlying BARWe process and can bee 
seen on the coarse grained level only (see example the 
GDK model, in which the kinks are spatially extended 
objects). This might have lead some studies to the con- 
jecture that the Z2 symmetry is a sufficient condition 
for the PC class (Hwang and Park, 1999). However the 
example of CDP (see Sect. IV. C) shows that this cannot 
be true. The field theory of BARWe models was investi- 
gated by (Cardy and Tauber, 1996, 1998). For this case 
the action 



S = 



d'^xdt 



+ o-(l-V^)V(/) 



A(l - V^2)0' + 



(153) 



is invariant under the simultaneous transformation of 
fields 
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TABLE XVII Critical exponents of BARWe. 
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FIG. 13 Local slopes (93) of the density decay in a bosonic 
BARW2 model. Different curves correspond to ct = 0.466, 
0.468 0.469, 0.47 (from bottom to top) (Odor and Menyhard, 
2002). 



tp —-0, (f) —(j) (154) 

Owing to the non-recurrence of random walks in c? > 2 
the system is in the active phase for a > and mean- 
field transition occurs with P = 1. However the survival 
probability of a particle cluster is finite for any ct > 
implying /3' — 0. Hence contrary to the DP class (i ^ (3' 
for d > 2. At c? = 2 random walks are barely recurrent 
and logarithmic corrections can be found. In this case 
the generalized hyperscaling law (Mendes et ai, 1994) is 
valid among the exponents 

2(^1 + ^^S' + 2r]' = dz. (155) 

In d = 1 however (3 — 13' holds owing to an exact duality 
mapping (Mussawisade et ai, 1998) and the hyperscaling 
is the same that of the DP (eq. 90). 

The RG analysis of BARWe for d < 2 run into difficul- 
ties. These stem from, the presence of another critical 
dimension d'^ = 4/3 (above which the branching reaction 
is relevant at cr = 0, and irrelevant for d < d'^) hence the 
d = I dimension cannot be accessed by controlled expan- 
sions from dc — 2. The truncated one- loop expansions 
(Cardy and Tauber, 1996) for d = 1 resulted in 

/3 = 4/7 , 3/7 , i^ii = 6/7 , Z = 2 (156) 

which are quite far from the numerical values determined 
by Jensen's simulations (Jensen, 1994) (Table XVII). 
Here the cluster exponents S and i] corresponding to the 
sector with even number of initial particles are shown. 



In case of odd number of initial particles they exchange 
values. 

It was conjectured (Deloubriere and van Wijland, 
2002) that in Id fermionic (single occupancy) and bosonic 
(multiple occupancy) models may have different critical 
behavior. Since only bosonic field theory exists, which 
gives rather inaccurate critical exponent estimates in 
(Odor and Menyhard, 2002) bosonic simulations were 
performed to investigate the density decay of BARW2 
from random initial state. Figure 13 shows the local 
slopes of density decay (oe// (93)) around the critical 
point for several branching rates (cr). The critical point 
is estimated at CTc = 0.04685(5), with the corresponding 
decay exponent a — 0.290(3). This value agrees with 
that of the PC class. 

If there is no explicit diffusion of particles besides the 
AA A ^ 3A processes (called DBAP by (Sud- 

bury, 1990)) an implicit diffusion can still be generated. 
By spatially asymmetric branching: A00 AAA and 
00yl AAA a diffusion may go on by two lattice steps: 
yl00 AAA 00A. As a consequence the decay pro- 
cess slows down, a single particle can not join a domain, 
hence domain sizes exhibit a parity conservation. This 
results in additional new sectors (besides the existing 
two BARW2 sectors) that depend on the initial condi- 
tions. For example in case of random initial distribution 
6 ~ 0.13(1) was measured by simulations (Hinrichsen and 
Odor, 1999a). Similar sector decomposition has been 
observed in diffusion of k-mer models (see for example 
(Barma and Dhar, 1994; Barma et ai, 1993)). 
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2. The NEKIM model 

An other important representative of PC class appear 
among nonequilibrium Ising models, in which the steady 
state is generated by kinetic processes in connection with 
heat baths at different temperatures (DeMasi et ai, 1985, 
1986; Droz et ai, 1989; Gonzalez-Miranda et at, 1987; 
Wang and Lebowitz, 1988). The research of them have 
shown that phase transitions are possible even in Id un- 
der nonequilibrium conditions (for a review see (Racz, 
1996)). In short ranged interaction models any non-zero 
temperature spin-flip dynamics cause disordered steady 
state. Menyhard proposed a class of general nonequilib- 
rium kinetic Ising models (NEKIM) with combined spin 
flip dynamics at T = and Kawasaki spin exchange dy- 
namics at T = cxD in which, for a range of parameters of 
the model, a PC-type transition takes place (Menyhard, 
1994). 

A general form (Glauber, 1963) of the Glauber spin-flip 
transition rate in one-dimension for spin = ±1 sitting 
at site i is: 



T ~ ( 1 

Wi = 2 + <5s,j_i.Si+i) ( 1 - -s,j(si_i + Sj+i 



(157) 



Here 7 = tanh {2J/kT), J denotes the coupling constant 
in the ferromagnetic Ising Hamiltonian, F and 6 are fur- 
ther parameters, which can in general, also depend on 
temperature. The Glauber model is a special case corre- 
sponding to ^ = 0, r = 1. There arc three independent 
rates: 

^TTT - ^(l + '5)(l-7), WiTi = ^(l + ^)(l+7) 



(158) 



In the NEKIM model T = is taken, thus 7 = 1, ^ttt = 
and r, S are the control parameters to be varied. 
The Kawasaki spin-exchange rate of neighboring spins 



Wii+i{Si,Si+i) = ^(l-SiSi+i)[l-^(Si_lSi-f Si+iSi+2)]- 

(159) 

At T = CX3 (7 = 0) the above exchange is simply an 
unconditional nearest neighbor exchange: 



1 



(160) 



where p^^ is the probability of spin exchange. The tran- 
sition probabilities in eqs.(157) and (160) are responsible 
for the basic elementary processes of kinks (K). Kinks 
separating two ferromagnetically ordered domains can 
carry out random walks with probability 



Prw oc 2i«|fx = r(l - S) , 



(161) 



while two kinks getting into neighboring positions will 
annihilate with probability 




FIG. 14 Phase diagram of the two-parameter model. The 
transversal dotted line indicates the critical point that was 
investigated in more detail (Menyhard and Odor, 2000). 



{will responsible for creation of kink pairs inside of 
ordered domains at T ^ 0). In case of the spin exchanges, 
which act only at domain boundaries, the process of main 
importance here is that a kink can produce two offspring 
at the next time step with probability 



PK-^ZK OC pex- 



(163) 



Pan OC W^i = +'5) 



(162) 



The abovementioned three processes compete, and it de- 
pends on the values of the parameters F, 5 and Pex what 
the result of this competition will be. It is important to 
realize that the process K 3K can develop into prop- 
agating offspring production only if Prw > Pan, i-e. the 
new kinks are able to travel on the average some lattice 
points away from their place of birth and can thus avoid 
immediate annihilation. It is seen from the above def- 
initions that 5 < is necessary for this to happen. In 
the opposite case the only effect of the K 3K process 
on the usual Ising kinetics is to soften domain walls. In 
the NEKIM model investigations the normalization con- 
dition Prw + Pan + Pk^3k = 1 waS Set. 

The phase diagram determined by simulations and 
GMF calculations (Menyhard and Odor, 1995, 2000) is 
shown in Fig. 14. The line of phase transitions separates 
two kinds of steady states reachable by the system for 
large times: in the Ising phase, supposing that an even 
number of kinks are present in the initial states, the sys- 
tem orders in one of the possible ferromagnetic states of 
all spins up or all spins down, while the active phase is 
disordered from the point of view of the underlying spins. 
The cause of disorder is the steadily growing number of 
kinks with time. While the low-level, iV = 1, 2 GMF so- 
lutions for the SCA version of NEKIM exhibit first order 
transitions, for > 2 this becomes continuous. GMF 
approximations (up to A^ = 6) with CAM extrapolation 
found /3 ~ 1 (Menyhard and Odor, 1995). Recent high 
precision Monte Carlo simulations (Menyhard and Odor, 
2000) resulted in critical exponents (3 = 0.95(2) (see Fig. 
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FIG. 15 /3e// = "^dTogT of kinks (circles) near the critical 
point (e = |5 — and linear extrapolation (dashed line) to 
the asymptotic value (/3 = 0.95(2)). Simulations were per- 
formed on a Id NEKIM ring of size L = 24000 (Menyhard 
and Odor, 2000). 



15) and S = 0.280(5) at the dotted line of the phase dia- 
gram. 

Ref. (Mussawisade et ai, 1998) have shown that an 
exact duality mapping exist in the phase diagram of the 
NEKIM: 



Prw 



/ 

Pan ~ Pan-, 

Pan 

+ 2p 

ex •) 

' Prw Pan • 



(164) 



The regions mapped onto each other have the same phys- 
ical properties. In particular, the hne p^x = maps onto 
the hne prw — Pan and the fast-diffusion hmit to the hmit 
Pex oo. There is a self-dual hne at 



6 = 



1 -Pe 



(165) 



By various static and dynamical simulations spin and 
kink density critical exponents have been determined in 
(Menyhard and Odor, 1996) and as the consequence of 
the generalized hyperscaling law for the structure factor 



S{0,t) = L[< Af2 > - < M >2] oc r 
and kink density 



) >cx r 



the exponent relation 



22/ 



(166) 



(167) 



(168) 



is established. Spins-clusters at the PC point grow by 
compact domains as in the Glauber point albeit with dif- 
ferent exponents (Menyhard and Odor, 1998). The spin- 
cluster critical exponents in magnetic field are summa- 
rized in Table (II). The global persistence {9g) and time 
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PC 


.00(1) 


.444(2) 
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TABLE XVIII Simulation data for static and dynamic criti- 
cal spin exponents for NEKIM. 



h 


0.0 


0.01 


0.05 


0.08 


0.1 


DP 


p 


1.0 


0.281 


0.270 


0.258 


0.285 


0.2767(4) 


1 




0.674 


0.428 


0.622 


0.551 


0.5438(13) 



TABLE XIX CAM estimates for the kink density and its 
fluctuation exponents. 



autocorrelation exponents (A) were determined both at 
the Glauber an at the PC critical points (Menyhard and 
Odor, 1997) and are shown in Table (XVIII). While at 
the Glauber point the scaling relation (17) is satisfied by 
these exponents it is not the case at PC criticality, there- 
fore the magnetization is non-Markovian process here. 

By applying an external magnetic field h that breaks 
the Z2 symmetry the transition type of the model changes 
to DP type (see Table XIX). (Menyhard and Odor, 1996) 

In (Menyhard and Odor, 2002) the symmetry of the 
spin updates was broken in such a way that two different 
types of domain walls emerged. The following changes to 
the Glauber spin-slip rates (158) with P = 1, (5 = were 
introduced: 



0, 



P+ < 1/2 , 



(169) 
(170) 



while the spin-exchange part remains the same. In the 
terminology of domain walls as particles the following 
reaction-diffusion picture arises. Owing to the sym- 
metry breaking there are two kinds of domain walls 
jt= A and tj= B, which can only occur alternately 
( ...B..A..B..A..B...A...) owing to the spin background. 
Upon meeting AB — s- happens, while in the oppo- 
site sequence, BA, the two domain-walls are repulsive 
due to (169). The spin exchange leads to A ABA 
and B <-> BAB type of kink reactions, which together 
with the diffusion of As and B-s leads to a kind of two- 
component, coupled branching and annihilating random 
walk (see Sect.V.K). There are two control parameters 
in this model: Pex that regulates the kink production- 
annihilation and p+ that is responsible for the local 
symmetry breaking (170). Simulations show that for 
Pex ~> 0, p+ < 0.5 an absorbing phase emerges with 
N-BARW2 class exponents (owing to the pairwise order 
of kinks hard-core effects cannot play a role), while the 
transition on the Pex > line belongs to the 1-1-1 d DP 
class. Since the AB reaction breaks the parity con- 
servation of species (but preserves the global parity con- 
servation) the necessary conditions for N-BARW2 class 
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can be cased. On the other hand the occurrence of the 
DP transition introduces a zero branching rate condition 
for N-BARW2 universal behavior. This study and the 
results for the generalized contact process (Sect. V.K.I) 
emphasize that the conditions for the N-BARW2 class 
should be further investigated. 

A generahzation is the probabilistic cellular automaton 
version of NEKIM, which consists in keeping the spin-flip 
rates given in eqs.(158) and prescribing synchronous up- 
dating. In this case the K 3K branching is generated 
without the need of additional, explicit spin-exchange 
process and for certain values of parameter-pairs (F, 6) 
with (5 < PC-type transition takes place. The phase 
boundary of NEKIM-CA in the (r, -5) plane is simi- 
lar to that in Fig. 14 except for the highest value of 
r = 1, (5c = cannot be reached, the limiting value is 
5c = -.065. 

An other possible variant of NEKIM was introduced 
in (Mcnyhard and Odor, 2000) in which the Kawasaki 
rate eq.(159) is considered at some finite temperature, 
instead of T = oo, but keeping T = in the Glauber- 
part of the rule. By lowering the temperature the spin 
exchange process acts against the kink production and a 
PC class transition occurs. In this case the active phase 
part of the phase diagram shrinks. For more details see 
((Odor et al, 1999)). 

The DS transition of this model coincides with the 
critical point and the scaling behavior of spin and kink 
damages is the same as that of the corresponding NEKIM 
variables (Odor and Menyhard, 1998). 



3. Parity conserving stochastic cellular automata 

The first models in which non-DP class transition to 
absorbing state was firmly established were Id stochastic 
cellular automata defined by (Grassberger et al., 1984). 
In these models the 00 and 11 pairs follow BARW2 dy- 
namics and their density is the order parameter that 
vanishes at some critical point. The PC class critical 
exponents estimated by simulations for these models by 
(Grassberger, 1989b). While in the "A" model the crit- 
ical point coincides with the DS transition point and both 
of them are PC type, in the "B" model the DS tran- 
sition occurs in the active phase - where the symmetry 
of replicas is broken - and therefore the DS exponents 
belong to the DP class (Odor and Menyhard, 1998). 

Another stochastic CA, which may exhibit PC class 
transition and which is studied later from different direc- 
tions is also introduced here. It points out that the under- 
lying BARW2 dynamics of domain walls in Z2 symmetric 
systems can sometimes be seen on coarse grained level 
only. A generalization of the Domany-Kinzel stochastic 
cellular automaton (Domany and Kinzcl, 1984) (sec Sec- 
tion IV. A. 2) was introduced by (Hinrichsen, 1997). This 
model has n+1 states per site: one active state Ac and 
n different inactive states /i , J2 ,...,/„ . The conditional 
updating probabilities are given by (A;, ^ = 1, . . . , n; k ^ 



I) 





P{Ik\h,h) 


= 1, 


(171) 


P{Ac\Ac, Ac) = \- 


nP{Ik\Ac, Ac) 


= q, 


(172) 


P{Ac\h,Ac:) 


= PiAc\Ac,Ik) 


= Pk , 


(173) 


P{h\Ik,Ac) 


^P{Ik\Ac,Ik) 








P{Ac\h,Ii) 


= 1, 


(174) 



and the symmetric case pi , . . . , p„ = p was explored. 
Equations (171)-(173) arc straightforward generaliza- 
tions of Eqs. (98)-(100). The only different process is 
the creation of active sites between two inactive domains 
of different colors in Eq. (174). For simplicity the prob- 
ability of this process was chosen to be equal to one. 

For n = 1 the model defined above reduces to the orig- 
inal Domany-Kinzel model. For n = 2 it has two Z2 
symmetrical absorbing states. The phase diagram of this 
model is very similar to that of the DK model (Fig. 6) 
except the transition line is PC type. If we call the re- 
gions separating inactive domains Ji and I2 as domain 
walls (denoted by K), they follow BARW2 process 

K^3K 2K^% KID^m (175) 

but for 1 > g > the size of active regions, hence 
the domain walls stays finite. Therefore a the obser- 
vation of the BARW2 process is not so obvious and 
can be done on coarse grained level only (except at 
the endpoint at g = 0, where active sites really look 
like kinks of the NEKIM model). Series expansions for 
the transition point and for the order parameter crit- 
ical exponent resulted in /? = 1.00(5) (Jensen, 1997) 
that is slightly higher than the most precise simula- 
tion results (Menyhard and Odor, 2000) but agrees with 
the GMF+CAM estimates (Mcnyhard and Odor, 1995). 
Similarly to the NEKIM model the application of ex- 
ternal symmetry breaking field changes the PC class 
transition into a DP class one (Hinrichsen, 1997). The 
other symmetrical endpoint (g = 1, p = ^) in the 
phase diagram again shows different scaling behavior 
(here three types of compact domains grow in competi- 
tion, the boundaries perform annihilating random walks 
with exclusion (see. Section V.B). 

Models with n > 3 symmetric absorbing states in Id 
do not show phase transitions (they are always active). 
In terms of domain walls as particles they are related 
to > 1 component N-BARW2 processes, which ex- 
hibit phase transition for zero branching rate only (see 
Sect.V.K) (Cardy and Tauber, 1996; Hooyberghs et al., 
2001). 

GDK type models - exhibiting n symmetric absorb- 
ing states - can be generalized to higher dimensions. In 
two dimensions Hinrichsen's spreading simulations for 
the n = 2 case yielded mean-field like behavior with 
6 = 1, r] = and z = 1 leading to the conjecture 
that 1 < dc < 2. A similar model exhibiting Potts- 
like Zn symmetric absorbing states in = 2 yielded 
similar spreading exponents but a first order transition 
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FIG. 16 Phase diagram of the MM from (Brown et al, 1997). 
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FIG. 17 Domain-wall dynamics in the interacting monomer- 
monomer model. Left part: branching; right part: annihila- 
tion 



(/3 = a = 0) forn = 2 (Lipowski and Droz, 2002a). In 
three dimensions the same model seems to exhibit mean- 
field like transition with (3=1. The verification of these 
findings would require further research. 

4. PC class surface catalytic models 

In this subsection I discuss some one-dimensional, sur- 
face catalytic-type reaction-diffusion models exhibiting 
PC class transition. Strictly speaking they are multi- 
component models, but I show that the symmetries 
among species enables us to interpret the domain-wall 
dynamics as a simple BARW2 process. 

The two-species monomer- monomer (MM) model was 
first introduced by (Zhuo et al., 1993). Two monomers, 
called A and adsorb at the vacant sites of a one- 
dimensional lattice with probabilities p and g, respec- 
tively, where p + q = l. The adsorption of a monomer at 
a vacant site is affected by monomers present on neigh- 
boring sites. If either neighboring site is occupied by 
the same species as that trying to adsorb, the adsorption 
probability is reduced by a factor r < 1, mimicking the 
effect of a nearest-neighbor repulsive interaction. Unlike 
monomers on adjacent sites react immediately and leave 
the lattice, leading to a process limited by adsorption 
only. The basic reactions are 

0^A, $^B, AB^$ (176) 

The phase diagram, displayed in Fig. 16 with p plotted vs. 
r, shows a reactive steady state containing vacancies bor- 
dered by two equivalent saturated phases (labeled A and 
B). The transitions from the reactive phase to either of 
the saturated phases are continuous, while the transition 
between the saturated phases is first-order discontinuous. 
The two saturated phases meet the reactive phase at a 
bi-critical point at a critical value of r = Vc. In the case 
of r = 1, the reactive region no longer exists and the only 
transition is a first-order discontinuous line between the 
saturated phases. Considering the density of vacancies 
between unlike species as the order parameter (that can 



also be called a species "C") the model is the so called 
"three species monomer-monomer model" . Simulations 
and cluster mean-field approximations were applied to in- 
vestigate the phase transitions of these models (Bassler 
and Browne, 1996, 1997, 1998; Brown et al., 1997). As 
Fig. 17 shows if we call the extended objects filled with 
vacancies between different species as domain-walls (C) 
we can observe C 3C and 2C BARW2 processes 
in terms of them. These C parity conserving processes 
arise as the combination of the elementary reaction steps 
(176). The reactions always take place at domain bound- 
aries, hence the Z2 symmetric A and B saturated phases 
are absorbing. 

The interacting monomer-dimer model (IMD) (Kim 
and Park, 1994) is a generalization of the simple 
monomer-dimer model (Ziff et al, 1986), in which parti- 
cles of the same species have nearest-neighbor repulsive 
interactions. The IMD is parameterized by specifying 
that a monomer (A) can adsorb at a nearest-neighbor 
site of an already-adsorbed monomer (restricted vacancy) 
at a rate ta^a with < va < 1, where Jca is an ad- 
sorption rate of a monomer at a free vacant site with 
no adjacent monomer-occupied sites. Similarly, a dimer 
{B2) can adsorb at a pair of restricted vacancies {B in 
nearest-neighbor sites) at a rate rsks with < rs < 1, 
where fcs is an adsorption rate of a dimer at a pair of free 
vacancies. There arc no nearest-neighbor restrictions in 
adsorbing particles of different species and the AB 
desorption reaction happens with probability 1. The case 
rA = fB = ^ corresponds to the ordinary noninteracting 
monomer-dimer model which exhibits a first-order phase 
transition between two saturated phases in one dimen- 
sion. In the other limiting case rA = rB= 0, there exists 
no fully saturated phase of monomers or dimers. How- 
ever, this does not mean that this model no longer has 
any absorbing states. In fact, there are two equivalent 
{Z2 symmetric) absorbing states in this model. These 
states comprise of only the monomers at the odd- or 
even-numbered lattice sites. A pair of adjacent vacan- 
cies is required for a dimer to adsorb, so a state with 
alternating sites occupied by monomers can be identified 
with an absorbing state. The PC class phase transition of 
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the rA = rB = infinite repulsive case has been investi- 
gated in (Hwang et ai, 1998; Kim and Park, 1994; Kwon 
and Park, 1995; Park et ai, 1995; Park and Park, 1995). 
As one can see the basic reactions are similar to those of 
the MM model (eq. (176)) but the order parameter here 
is the density of dimers (K) that may appear between 
ordered domains of alternating sequences: 'OAO..AOA.' 
and 'AOA..OAO', where monomers are on even or odd 
sites only. The recognition of an underlying BARW2 
process (175) is not so easy in this case, still considering 
regions between odd and even filled ordered domains one 
can identify domain wall random-walk, annihilation and 
branching processes through the reactions with dimers as 
one can see on the examples below. The introduction of 
Z2 symmetry-breaking field, that makes the system pre- 
fer one absorbing state to the other was shown to change 
that transition type from PC to DP (Park and Park, 



1995). 
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FIG. 18 log(pfe(t)) versus log(t) in NEKIM simulations for 
a = 0, 0.1, 0.2..., 1 initial conditions (from bottom to top 
curves) (Menyhard and Odor, 2000). 
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5. NEKIM with long-range correlated initial 
conditions 

The effect of initially long-range correlations has al- 
ready been discussed in 2d Ising models (Sect.III.A.5) 
and in Id bond percolating systems (Sect. IV. A. 7). 
In both cases continuously changing decay exponents 
have been found. In case of the NEKIM model (see 
Sect. IV.D.2) simulations (Menyhard and Odor, 2000; 
Odor et al., 1999) showed that the density of kinks Pk{t) 
changes as 

Pk{t) oc (177) 

by starting the system with two-point correlated kink 
distributions of the form (122). The k((t) changes lin- 
early between the two extremes f3/i^\\ ~ ±0.285 shown 
in Fig. 18. This behavior is similar to that of the DP 
model case (Sect. IV. A. 7), but here one can observe a 
symmetry: 

(7 <-> 1 — (T, K — K (178) 

which is related to the duality symmetry of the NEKIM 
model (164). 

6. GDK with spatial boundary conditions 

The surface critical behavior of the PC class has been 
explored through the study of the GDK model (Sect. 
IV.D.3) (Frojdh et ai, 2001; Howard et ai, 2000; Laurit- 
sen et ai, 1998). The basic idea is that on the surface one 



FIG. 19 Schematic mean field boundary phase diagram for 
BARW from (Frojdh et ai, 2001). See text for an explanation 
of the labeling. 

may include not only the usual BARW2 reactions (175) 
but potentially also a parity symmetry breaking A ^ 
reaction. Depending on whether or not the A — > re- 
action is actually present, we may then expect different 
boundary universality classes. Since the time reversal 
symmetry (88) is broken for BARW2 processes two inde- 
pendent exponents (/3i,sccd, /3i,dcns) characterize the sur- 
face critical behavior. 

The surface phase diagram for the mean field theory of 
BARW (vahd for d > dc = 2) is shown in Fig. 19. Here 
fm, <^ms the rates for the branching processes A — > 
{m + 1)A in the bulk and at the surface, respectively, and 
/is is the rate for the surface spontaneous annihilation 
reaction A — > 0. Otherwise, the labeling is the same as 
that for the DP phase diagram (see Figure 7). The fig > 
corresponds to the parity symmetry breaking RBC. 

For fig = (IBC) parity conserving case the surface 
action is of the form 

Ss = J d'^-^X^^ dt ^ (72,, (1 - i^fi^s<t's) , (179) 

where tps = V'(xj[,a::_L = 0, i) and (f)s — (/>(x||,xj_ = 0, i). 



34 





Id (IBC) 


Id (RBC) 


2d (O) 


2d (Sp) 


/^l jsccd 


2.06(2) 


1.37(2) 








/3l jdcns 


1.34(2) 


2.04(2) 


3/2 


1 


Tl 


1.16(4) 


1.85(4) 


1 


1 


71 


2.08(4) 


2.77(4) 


1/2 


1/2 



^ I ^ TABLE XX Critical boundary exponents of the PC class in 

Us M-s d = 1, 2 for ordinary and special cases. 



FIG. 20 Schematic surface phase diagrams for BARW in d = 
1 for (a) am < cr^,criticai, and (b) Um = cfm.criticai (Frojdh 
et al., 2001). See text for an explanation of the labeling. 



In = 1 the boundary and bulk transitions inaccessible 

to controlled perturbative expansions, but scaling analy- 
sis shows that surface branching is irrelevant leading to 
the Sp* and Sp special transitions. For the /is > (RBC) 
parity symmetry breaking case the surface action is 

//>T ^ 

+M.(V« - l)4>e\- (180) 

and the RG procedure shows that the stable fixed point 
corresponds to the ordinary transition. Therefore in Id 
the phase diagram looks very differently from the mean- 
field case (Fig. 20). One can differentiate two cases cor- 
responding to (a) the annihilation fixed point of the bulk 
and (b) the PC critical point of the bulk. As on can see 
in both cases the ordinary transition (0,0*) corresponds 
to /is > 0, RBC and the special transitions {Sp,S-p*) to 
/tg = 0, IBC. The ABC condition obviously behaves as 
if there existed a surface reaction equivalent to ^ A, 
and thus it belongs to the normal transition universal- 
ity class. By scaling considerations the following scaling 
relations can be derived: 

n = !^|| -/3l,dens , (181) 

-I- dvji_ = ^i,seed + /3dens + 7l • (182) 

Rcf. (Howard et al., 2000) showed that on the self-dual 
line of the Id BARWc model (see Sect. IV.D.2 and (Mus- 
sawisade et al., 1998)) the scaling relations between ex- 
ponents of ordinary and special transitions 

/^°seed = /3l,dens (183) 

and 

/^Mccd^/^Mcn.- (184) 

hold. Relying on universality they claim that they should 
be valid elsewhere close to the transition line. Numerical 
simulations support this hypothesis as shown in Table 
XX. 



E. Branching with kA — > annihilation 

In Sect. IV.C.l the 2A annihilating random 
walk (ARW) has already been introduced. By adding 
branching processes to it wc defined the BARWo and 
BARWe models exhibiting continuous phase transitions 
(see Sects. IV.A.3 and IV.D.l). Now we generalize this 
construction to m-branching and kA — > annihilation 
type of models (BkARW) formulating it by the field the- 
oretical action: 

S = j d^xdt[jp{dt - DV'^)(I) - A(l - t{j'')(j)'' 

+a{l - V'")V"/'] • (185) 

These systems for fc > 2 exhibit an upper critical dimen- 
sion: dc = 2/{k-l) (Cardy and Tauber, 1998; Lee, 1994) 
with the mean-field exponents 

a = P = l/{k-l), Z = 2, uj_ = l/2. (186) 

At dc (which falls below physical dimensions for k > 3) 
the decay has logarithmic corrections: 

p{t)=Ak{ln{t)/tf/^'^-'^ . (187) 

So for the AAA process in one dimension this gives 
the decay behavior (Lee, 1994): 

p{t) = {^-^)HHt)/t)i+0{ti) , (188) 

which is the dominant behavior of the Id bosonic PCPD 
model at the transition point (see Sect.V.F). Note that 
for fc = 3, m = 1.2 the field theory of (Cardy and Tauber, 
1998) predicts DP class transitions in d = 1 owing to the 
BARWo terms generated by the renormalization. 

F. General nA —> {n + k)A, mA —^{m — l)A processes 

For a long time reaction-diffusion models with only 
single parent branching have been investigated, although 
there was and early forgotten numerical result by (Grass- 
berger, 1982c) claiming a non-DP type of continuous 
phase transition in a model where particle production can 
occur by the reaction of two parents. Later it turned out 
than in such binary production lattice models, where soli- 
tary particles follow random walk (hence they behave like 
a coupled system) different universal behavior emerges 
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indeed (sec Sect. V.F). In this subsection I shall discuss 

the mean-field classes in models 

nA-^{n + k)A, mAA{m- l)A, (189) 

with n > 1, m > 1, fc > 0, Z > and m — l > 0. In low di- 
mensions the site restricted and bosonic versions of these 
models exhibit different behavior. The field theory of 
the n = m = I = 2, bosonic model was investigated in 
(Howard and Tauber, 1997) and concluded non-DP type 
of criticality with dc = 2 (see Sect. V.F). For other cases 
no rigorous field theoretical treatment exist. Numerical 
simulations in one and two dimensions for various n = 3 
and n = 4 models resulted in somewhat contradictory re- 
sults (Kockelkoren and Chate, 2003a; Odor, 2003a; Park 
et ai, 2002). There is a disagreement in the value of the 
upper critical dimension, but in any case dc seems to be 
very low {d^ = 1 — 2) hence the number of non-mean-field 
classes in such models is limited. On contrary there is a 
series of mean-field classes depending on n and m. 

1. The n = m symmetric case 

In such models there is a continuous phase transition 
at finite production probability: Uc = characterized 

by the order parameter exponents (Odor, 2003a; Park 
et al, 2002): 

/3^^ = 1, aMF = 1/n . (190) 

These classes generalize the mean-field class of the DP 
(Sect. IV. A) and binary production models (Sect. V.F). 
For referencing it in Table XXXII of the Summary I call 
it as: PARWs (symmetric production and m-annihilating 
random walk class). 

2. The n > m symmetric case 

In this case the mean-field solution provides a first or- 
der transition (see (Odor, 2003a)), hence it does not im- 
ply anything with respect to possible classes for models 
below the critical dimension (rf < dc). Note however, that 
by higher order cluster mean- field approximations, when 
the diffusion plays a role the transition may turn into a 
continuous one (see for example (Menyhard and Odor, 
1995; Odor et al., 1993; Odor and Szolnoki, 1996)). 

3. The n < m case 

In this case the critical point is at zero production 
probability CTc = 0, where the density decays as = 

l/(m — 1) as in case of the n = 1 branching and m = I an- 
nihilating models (BkARW classes sec Sect. IV. E), but 
the steady state exponent is different: (3^^^ = l/{m — n), 
defining an other series of mean-field classes (PARWa) 
(Odor, 2003a). Again cluster mean-field approximations 



may predict the appearance of other a > transitions 
with different critical behavior. 



V. UNIVERSALITY CLASSES OF MULTI-COMPONENT 
SYSTEMS 

First I recall some well known results (see refs. in 
(Privman, 1996)) for multi-component reaction-diffusion 
systems without particle creation. From the viewpoint 
of phase transitions these describe the behavior in the 
inactive phase or in case of some N-component BARW 
models right at the critical point. Then I show the ef- 
fect of particle exclusions in Id. Later universal behavior 
of more complex, coupled multi-component systems are 
discussed. This field is quite new and some of the results 
are still under debate. 



A. The ^-1-5^0 classes 

The simplest two-component reaction-diffusion model 
involves two types of particles undergoing diffusive ran- 
dom walks and reacting upon contact to form an inert 
particle. The action of this model is: 

S = J d''xdi[^jA{di - DAy^)'pA + ^B{dt - DbV^)4>b 

- \{l-ijAi^B)<t>A<t>B]~ p{^){^A(0) + 4>Bm (191) 

where Da and Db denotes the diffusion constants of 
species A and B. In d < dc = 4 dimensions and for 

homogeneous, initially equal density of A and B parti- 
cles (po) the density decays asymptotically as (Bramson 
and Lebowitz, 1988; Burlatskii and Ovchinnikov, 1978; 
Kang and Rcdner, 1985; Lee and Cardy, 1995; Ovchin- 
nikov and Zel'dovich, 1978; Toussaint and Wilczek, 1983) 

PA{t) = pB{t) oc Csf^t-"!^ , (192) 

where A = p(0) - C"/9'*/2(0) + C is a universal C 
is a non- universal constant. This slow decay behavior is 
due to the fact that in the course of the reaction, local 
fluctuations in the initial distribution of reactants lead to 
the formation of clusters of like particles that do not react 
and they will be asymptotically segregated for < 4. The 
asymptotically dominant process is the diffusive decay of 
the fluctuations of the initial conditions. Since this is a 
short ranged process the system has a long-time memory 
- appearing in the amplitude dependence - for the initial 
density p(0). For d < 2 controlled RG calculation is 
not possible, but the result (192) gives the leading order 
term in e = 2 — d expansion. For Da ^ Db case a RG 
study (Lee and Cardy, 1995) found new amplitude but 
the same exponents. 

The persistence behavior in Id with homogeneous, 
equal initial density of particles (po = Pa{^) + Pb{^)) was 
studied by (O'Donoghue and Bray, 2001). The probabil- 
ity p{t), that an annihilation process has not occurred 
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at a given site ("type I persistence") has the asymptotic 
form 

p{t) ~ const. + . (193) 

For a density of particles p » 1, 9i is identical to 
that governing the persistence properties of the one- 
dimensional diffusion equation, where 9i w 0.1207. In 
the case of an initially low density, pq << 1, 9i 1/4 
was found asymptotically. The probability that a site 
remains unvisited by any random walker ("type II per- 
sistence" ) decays in a stretched exponential way 

pit) ~ exp(-const.xpo/^i^/*) (194) 
provided po << 1. 



B. AA 0, BB with hard-core repulsion 

The next simplest two-component model in which par- 
ticle blocking may be effective in low dimensions was 
investigated first in the context of a stochastic cellular 
automaton model. At the symmetric point of the GDK 
model (see Sect. IV. D. 3) compact domains of II and 
12 grow separated hy A = Ac — II and B — Ac ~ 12 
kinks that cannot penetrate each other. In particle lan- 
guage this system is a reaction-diffusion model of two 
types A + A^ <l>, B + B ^ <!) with exclusion AB BA 
and special pairwise initial conditions (because the 
domains are bounded by kinks of the same type): 

A. . .A. . .B.B. .B B. .A. .A. . 

In the case of homogeneous, pairwise initial conditions 
simulations (Odor and Menyhard, 2000) showed a den- 
sity decay of kinks p oc characterized by a power-law 
with an exponent somewhat larger than a = 0.5. The 
a = 1/2 would have been expected in case of two copies 
of ARW systems which do not exclude each other. Fur- 
thermore the deviation of a from 1/2 showed an initial 
density dependence. Ref. (Odor and Menyhard, 2000) 
provided a possible explanation based on permutation 
symmetry between types, according to which hard-core 
interactions cause marginal perturbation resulting in 
non-universal scaling. The situation is similar to that 
of the compact directed percolation that is confined by 
paraboHc boundary conditions (see Sect. IV. C. 3) if we as- 
sume that AB and BA pairs exert parabolic space-time 
confinement on coarsening domains. Non-universal scal- 
ing can also be observed at surface critical phenomena. 
Similarly here AB, BA pairs produce 'multi-surfaces' in 
the bulk. However simulations and independent interval 
approximations in a similar model predict logarithmic 
corrections to the single component decay with the form 
p ~ t-^/Vln(t) (Majumdar et al, 2001). Note that both 
kind of behavior may occur in case of marginal perturba- 
tions. Cluster simulations (Odor and Menyhard, 2000) 



also showed initial density {pii{Q)) dependent survival 
probability of /2-s in the sea of /1-s: 

Pi2(t)(xt-^^P''^'^^K (195) 

However this reaction-diffusion model with homoge- 
neous, random initial distribution of As and B-s ex- 
hibits a much slower density decay. An exact duality 
mapping helps to understand the coarsening behavior. 
Consider the leftmost particle, which may be either A or 
B, and arbitrarily relabel it as a particle of species X. 
For the second particle, we relabel it as Y if it is the 
same species as the initial particle; otherwise wc relabel 
the second particle as X. We continue to relabel each 
subsequent particle according to this prescription until 
all particles are relabeled from {A,B} to {X,Y}. For 
example, the string 

AABABBBA--- 

translates to 

XYYYYXYY ■■■ . 

The diffusion of the original A and B particles at equal 
rates translates into diffusion of the X and Y parti- 
cles. Furthermore, the parallel single-species reactions, 
A + A ^ % and B + B 0, translate directly to two- 
species annihilation X + Y ^ (d (see Sect. V.A) in 
the dual system. The interesting point is that in the 
X -h F ^ model blockades do not exist, because XY 
pairs annihilate, and there is no blockade between XX 
and YY pairs. Therefore the density decay should be 
proportional to t~'^^^. Simulations confirmed this for the 
A + A^ <D, B + B ^ 9 (Odor, 2001c) model, neverthe- 
less corrections to scaling were also observed. The pair- 
wise initial condition transforms in the dual system to 
domains of ..XYXY.. separated by YY and XX pairs, 
which do not allow X and Y particles to escape each 
other. 



C. Multi-species Ai + Aj <l) classes 

By generalizing the diffusion-limited reactions AB 
(Sect. V.A) for g > 2 species 

Ai +Aj^iD (196) 

in d > 2 dimensions the asymptotic density decay for 
such mutual annihilation processes with equal rates and 
initial densities is the same as for single-species pair an- 
nihilation A A —^ 0. 

In = 1, however, particles of different types can not 
pass each other and a segregation occurs for all q < oo. 
The total density decays according to a g dependent 
power law, p oc with 

a={q-l)/2q (197) 
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exactly (Dcloubricrc et al, 2002). These findings were 
also supported through Monte Carlo simulations. Spe- 
cial initial conditions such as ...ABCDABCD... prevent 
the segregation and lead to decay of the 2A ^ $ model 
(Sect.IV.C.l). 



D. Unidirectionally coupled ARW classes 

Besides the symmetrically coupled ARW systems dis- 
cussed above now I turn to unidirectionally coupled ARW 
(Sect.IV.C.l) models: 



A + A 
B + B 
C + C 



A^ A + B 
B-f B + C 
C^C + D 



(198) 



introduced and analyzed with RG technique and simula- 
tions by (Goldschmidt, 1998). This kind of coupling was 
chosen because A B would constitute a spontaneous 
death process of A particles leading to exponential den- 
sity decay. On the other hand quadratic coupling of the 
form A + A ^ B + B leads to asymptotically decoupled 
systems (Howard and Tauber, 1997). The mean- field 
theory is described by the rate equation for the density 
Pi{x, t) at level i: 



dpi{x,t) 
dt 



(199) 

which relates the long time behavior of level i to level 



1/2 

ni[t) oc n^l^ 



(200) 



By inserting into this the exact solution for ARW 
(Sect.IV.C.l) one gets: 




for d < 2 , 

for d = dc = 2 , 

for d > 2 , 



(201) 



The action of a two-component system with fields a, a, b, 
b for equal annihilation rates (A) takes the form 

S = J d'^x J dt a{dt - DV'^)a - A(l - a^)a^ + 

+b{dt - DV'^)b - A(l - 6^)6^ + a{l - b)aa\ (202) 

The RG solution is plagued by IR-divergent diagrams 
similarly to UCDP (see Sect. VI. H) that can be inter- 
preted as eventual non-universal crossover to the decou- 
pled regime. Simulation results exhibiting finite par- 
ticle numbers and coupling strengths - really show the 
breakdown of scaling (201), but the asymptotic behavior 
could not be determined. Therefore the results (201) are 
valid for an intermediate time region. 



E. DP coupled to frozen field classes 

One of the first generalizations of absorbing phase tran- 
sition models were such systems, which exhibit many 
absorbing states, hence the conditions of DP hypothe- 
sis (Grassberger, 1982a; Janssen, 1981) (Sect. IV) are 
not satisfied. Several variants of models with infinitely 
many absorbing states containing frozen particle config- 
urations have been introduced. The common behavior 
of these models that non-diffusive (slave) particles are 
coupled to a DP-like (order parameter) process. In case 
of homogeneous, imcorrelated initial conditions DP class 
exponents were found, on the other hand in cluster simu- 
lations - that involve correlated initial state of the order 
parameter particles - initial density dependent scaling 
exponents (77 and S) arise. These cluster exponents take 
the DP class values only if the initial density of the slave 
particles agrees with the "natural density" that occurs 
in the steady state. The first such models introduced by 
(Jensen, 1993a; Jensen and Dickman, 1993a) were the so 
called pair contact process (PCP) (see Sect. V.E.I) and 
dimer the reaction model. These systems seem to be sin- 
gle component ones by defining rules for the pairs, but 
the isolated, frozen particles behave as a second compo- 
nent. In the threshold transfer process (TTP) (Mendes 
et al., 1994) the two components defined as: the '2'-s fol- 
lowing DP process and the 'I'-s which decay or reappear 
as: ^ 1 0. 

For the PCP model defined by the simple processes 
(208) a set of coupled Langevin equations were set up 
(Mufioz et al., 1996, 1998) for the fields ni(x,i) and 
n2(x, f) (order parameter): 

= [r2 + DiV^ - uini - win{\ni -\- ^Jniqi 



dt 



[n + DiV^ - itma - wmi] + (203) 



where Di,ri,Ui, and Wi are constants and 771 (x,f) and 
772 (x,t) are Gaussian white noises. Owing to the mul- 
tiple absorbing states and the lack of the time reversal 
symmetry (88) a generalized hyperscaling law (155) has 
been derived by (Mendes et al., 1994). As discussed in 
(Munoz et al., 1998) this set of equations can be simpli- 
fied by dropping the Di, ui, and noise terms in the ni 
equation, and then solving that equation for m in terms 
of n2- Substituting in the n2 equation, one obtains 



(9712 (X,t) 
dt 



= £'2V^n2(x, t) + m2n2(x, t) — U2nl{x, t) + 

+ W2(ri/wi - ni(x,0))n2(x,i)e Jo 

+ V«2(x,t)772(x,t) , (204) 



where ni(x, 0) is the initial condition of the rii field, and 
TO2 = r2 — W2ri/wi. The "natural density" (Jensen and 
Dickman, 1993a), n""*, then corresponds to the uniform 
density, ni{t = 0) = ri/wi, for which the coefficient 
of the exponential term vanishes, and we get back the 
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Langevin equation of DP (84) . This derivation provides 
a simple explanation for the numerical observation of DP 
exponents in case of natural initial conditions. However 
it does not take into account the long-time memory and 
the fluctuations of passive particles (with power-law time 
and p dependences (Odor et ai, 1998)). Therefore, some 
of the terms omitted in this derivation (as for instance 
the term proportional to in the equation for rii) can- 
not be safely ehminated (Marques et al., 2001) and this 
simplified theory does not generate critical fluctuations 
for its background fleld. To overcome these problems a 
more rigorous field theoretical analysis involving path in- 
tegral representation of a two-component variant of this 
class (with M' activity and 'i?' slave): 



A^B 



(205) 



has been done (Mufioz et ai, 2002) and provided some 
evidence that the homogeneous state critical properties 
of the activity field are DP like irrespectively of the crit- 
icality of the slave field. For this model (Muhoz et ai, 
2002) could also prove the numerical result (Odor et at, 
1998) that the slave fleld decays as 



(206) 



However this treatment has still not provided theoreti- 
cal proof for the initial density dependent spreading ex- 
ponents observed in simulations (Jensen, 1993a; Jensen 
and Dickman, 1993a; Mendes et ai, 1994; Odor et ai, 
1998) and by the numerical integration of the Langevin 
equation (Lopez and Munoz, 1997). Furthermore the sit- 
uation is much more complicated when approaching crit- 
icality from the inactive phase. In particular, the scal- 
ing behavior of ua in this case seems to be unrelated to 
ub (this is similar to the diffusive slave field case (Odor 
et ai, 2002) Sect. V.F.3). In this case it is more diffi- 
cult to analyze the field theory and dynamical percola- 
tion type of terms are generated that can be observed 
in 2d by simulations and by mean-field analysis (Muhoz 
et at, 1996, 1998; Wijland, 2002). Very recently it was 
claimed, based on the field theoretical analysis of the 
GEP model (Sect. IV. B) - that exhibits similar long-time 
memory terms - that the cluster variables should follow 
stretched exponential decay behavior (J.-Dalmaroni and 
Hinrichsen, 2003). 

In two dimensions the critical point of spreading (ps) 
moves (as the function of initial conditions) and do not 
necessarily coincide with the bulk critical point {pc)- The 
spreading behavior depends on the coefficient of the ex- 
ponential, non-Markovian term of (204). For positive 
coefficient the Ps falls in the inactive phase of the bulk 
and the spreading follows dynamical percolation (see 
Sect. IV. B). For negative coefficient the Ps falls in the 
active phase of the bulk and spreading exponents are 
non-universal (like in Id) but satisfy the hyper-scaling 
(155). 

Simulations and GMF analysis (Marques et ai, 2001; 
Mufioz et ai, 2002) in the inactive phase in Id showed. 
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FIG. 21 Initial concentration dependence of the exponent rj 
for PGP model (Odor et ai, 1998). Linear regression gives a 
slope 0.320(7) between r) - -qop and pi(0) - pT^ 



that the steady state density of slave particles approaches 
the natural value by a power-law 



I nat 
\Pl 



Pi I « \P-Pc\ 



Pi 



(207) 



with /3i ~ 0.9 for PCPD and /3i = 1 for TTP transfer 
models. This difference might be the consequence that 
in TTP models the slave field fluctuates and relaxes to r 
quickly, while in the PCP case it is frozen. Due to this 
fact one might expect initial condition dependence here 
again. 



1. The pair contact process model 

Up to now I discussed spreading processes with unary 
particle production. Now I introduce a family of systems 
with binary particle production (i.e for a new particle 
production two particles are needed to collide). The PCP 
model is deflned on the lattice by the following processes: 



2A 3A, 



2A 



(208) 



such that reactions take place at nearest-neighbor (NN) 
sites and we allow single particle occupancy at most. The 
order parameter is the density of NN pairs p2- The PCP 
exhibits an active phase for p < pc] for p > Pc the system 
eventually falls into an absorbing configuration devoid 
of NN pairs, (but that contains a density pi of isolated 
particles). The best estimate for the critical parame- 
ter in one dimension is Pc = 0.077090(5) (Dickman 
and da Silva, 1998). Static and dynamic exponents cor- 
responding to initially uncorrelated homogeneous state 
agree well with those of 1+1 d DP (XII). Spreading expo- 
nents that involve averaging over all runs, hence involving 
the survival probability are non- universal (see Fig. 21) 
(Odor et ai, 1998). The anomalous critical spreading of 
PCP can be traced to a long memory in the dynamics 
of the order parameter p2 , arising from a coupling to an 
auxiliary field (the local particle density p) , that remains 
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frozen in regions where p2 = 0. In (Odor et al., 1998) 
a slight variation of the spreading critical point (as the 
function of /Oi(0)) was observed (similarly to the 2d case), 
but more detailed simulations (Dickman, 1999) suggest 
that this can be explained by strong corrections to scal- 
ing. Simulations provided numerical evidence that the pi 
exhibits anomalous scaling as eq.(207) to p""** = 0.242(1) 
with a DP exponent (Odor et al., 1998) for p < pc and 
with /3i = 0.9(1) for p > p,, (Marques et al, 2001). The 
DS transition point and the DS exponents of this model 
coincide with the critical point and the critical exponents 
of the PGP (Odor et al. 1998). 

The effect of an external particle source that cre- 
ates isolated particles, hence do not couple to the or- 
der parameter was investigated by simulations and by 
GMF-I-CAM approximations (Dickman et al, 2001b). 
While the critical point pc shows a singular dependence 
on the source intensity, the critical exponents appear to 
be unaffected by the presence of the source, except pos- 
sibly for a small change in /3. 

The properties of the two dimensional PCP in case 
of homogeneous, uncorrelated initial conditions was in- 
vestigated by simulations (Kamphorst et al, 1999). In 
this case all six NN of a pair was considered for reactions 
(208). The critical point is located at Pc = 0.2005(2). By 
determining a, (3/v±_ and Z exponents and order parame- 
ter moment ratios by simulations, the universal behavior 
of the 2-1-1 dimensional DP class (XII) was confirmed. 
The spreading exponents are expected to behave as de- 
scribed in Sect. V.E. 



F. DP with coupled diffusive field classes 

The next question one can pose following Section V.E 
is whether a diffusive field coupled to a DP process is 
relevant. Prominent representatives of such models are 
binary particle production systems with explicit diffu- 
sion of solitary particles. The critical behavior of such 
systems are still under investigations. The annihilation- 
fission (AF) process is defined as 

2A -^0, 2A^{n + 2)A, A0 <^ . (209) 

The corresponding action for bosonic particles was de- 
rived from master equation by (Howard and Tauber, 
1997) 

S = j d'^xdt[ij{dt - D\/^)4> - A(l - tp^)(j)^ 

+a{l - ^'')^Ij^4>^] . (210) 

Usually bosonic field theories describe well the critical 
behavior of "fermionic" systems (i.e. with maximum one 

particle per site occupation). This is due to the fact that 
at absorbing phase transitions the occupation number 
vanishes. In this case however the active phase of bosonic 
and fermionic models differ significantly: in the bosonic 
model the particle density diverges, while in the fermionic 



model there is a steady state with finite density. As a 
consequence the bosonic field theory cannot describe the 
active phase and the critical behavior of the fermionic 
particle system. 

As one can see this theory lacks interaction terms linear 
in the field variable (p (masslcss) in contrast with the DP 
action (87). Although the field theory of bosonic AF 
has turned out to be non-renornmalizible, (Howard and 
Tauber, 1997) concluded that critical behavior cannot be 
in DP class. In fact the upper critical behavior is = 2 
that is different from that of DP and PC classes (Odor 
et al, 2002). In the Langevin formulation 

^^^^1^ = DV'^pix, t) + (ncr - 2\)p^{x, t) + p{x, t)r]{x, t) 

(211) 

the noise is complex: 

< r^{x,t) > = (212) 
< r]{x,t)r]{x' ,t') > = [n{n + 3)a-2X]6'^{x-x'){t-t') , 

that is again a new feature (it is real in case of DP and 
purely imaginary in case of CDP and PC classes). 

Eq.(211) without noise gives the MF behavior of the 
bosonic model: for na > 2 A the density diverges, while 
for na < 2X it decays with a power-law with ~ 1. 

The MF description in the inactive phase of the bosonic 
model was found to be valid for the 2d fermionic AF 
system too (Odor et al, 2002). Here the pair density 
decays as p2{t) ex (Odor et al, 2002) in agreement 
with the MF approximation. Contrary to this for A < 
Ac p and p2 seem to be related by a logarithmic ratio 
p{T)/ p2(t) oc ln(t). This behavior could not be described 
by the mean-field approximations. 

The field theory suggests that the critical behavior of 
the Id bosonic model in the inactive phase is domi- 
nated by the 3A ^ process, that has an upper critical 
dimension dc = 1. Therefore in one dimension the parti- 
cle density decays with a power-law exponent: a = 1/2 
with logarithmic corrections (see Sect. IV. E). This be- 
havior has been confirmed by simulations in case of the Id 
annihilation- fission model (Odor and Menyhard, 2002). 

The field theoretical description of the fermionic AF 
process run into even more serioiis difficulties (Tauber, 
2000) than that of the bosonic model and predicted an 
upper an critical dimension dc = 1 that contradicts sim- 
ulation results (Odor et al, 2002). For the fermionic AF 
system mean-field approximations (Carlon et al., 2001; 
Odor et al, 2002) give a continuous transition with ex- 
ponents 

/3=1, 13' = 0, Z = 2, u\\=2, a = 1/2, ?? = 0. 

(213) 

These MF exponents are distinct from those of other well 

known classes (DP,PC,VM ...). They were confirmed in 
a 2d fermionic AF model, with logarithmic corrections, 
indicating dc = 2 (Odor et al, 2002). An explanation 
for the new type of critical behavior based on symme- 
try arguments are still missing but numerical simulations 
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suggest (Hinrichsen, 2001c; Odor, 2000) that the behav- 
ior of this system can be described (at least for strong 
diffusion) by coupled sub-systems: single particles per- 
forming annihilating random walk coupled to pairs (B) 
following DP process: B 2B, B %. The model has 
two non-symmetric absorbing states: one is completely 
empty, in the other a single particle walks randomly. Ow- 
ing to this fluctuating absorbing state this model does 
not oppose the conditions of the DP hypothesis. It was 
conjectured by (Henkel and Hinrichsen, 2001) that this 
kind of phase transition appears in models where (i) soli- 
tary particles diffuse, (ii) particle creation requires two 
particles and (iii) particle removal requires at least two 
particles to meet. The exploration of other conditions 
that affect these classes are still under investigation. 



1. The PCPD model 

A PCPD like model was introduced in an early work 
by (Grassberger, 1982c). His preliminary simulations in 
Id showed non-DP type transition, but the model has 
been forgotten for a long time. The PCPD model in- 
troduced by (Carlon et ai, 2001) is controlled by two 
parameters, namely the probability of pair annihilation 
p and the probability of particle diffusion D. The dy- 
namical rules are 



AAd), %AA AAA 
AA^m 
A0 ^ M 



with rate (l-p)(l-D)/2 

with rate p(l — D) 

with rate D . (214) 



The mean-field approximation gives a continuous transi- 
tion at p = 1/3. For p < Pc{D) the particle and pair 
densities exhibit singular behavior: 

p(oo) oc {pc - p)'^ P2(oo) cx {pc - pY^ (215) 

while aX p — Pc{D) they decay as: 

p(t)cxt-", p2{t)oct-"' , (216) 

with the exponents: 

a = 1/2, a2 - 1, /? - 1, /32 = 2 . (217) 

According to pair mean-field approximations the phase 
diagram can be separated into two regions (see Fig. 22). 
While for D > 1/7 the pair approximation gives the same 
Pc{D) and exponents as the simple MF, for D < 1/7-s 
the transition line breaks and the exponents are different 



1, 



a2 



1, f3 = l, f32 



1 



(218) 



In the entire inactive phase the decay is characterized by 
the exponents: 



1, a2=2. 



(219) 
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FIG. 22 Schematic phase diagram of the Id PCPD model. 
Circles correspond to simulation and DMRG results, solid line 
to mean-field, dashed line to pair-approximation. 
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The DMRG (Carlon et ai, 2001) method and simulations 
of the Id PCPD model (Hinrichsen, 2001b) resulted in 



TABLE XXI Summary of resuhs for Id, n = 1 AF model. 
The non-universal critical parameter pc of the parallel model 
is shown here. 



agreeing Pc{D) values but for the critical exponents no 
clear picture was found. They could not clarify if the 
two distinct universality suggested by the pair mean-field 
approximations was really observable in the Id PCPD 
model. It is still a debated topic whether one new class, 
two new classes or continuously changing exponents oc- 
cur in Id. Since the model has two absorbing states (be- 
sides the vacuum state there is another one with a sin- 
gle wandering particle) and some exponents were found 
to be close to those of the PC class {Z = 1.6 — 1.87, 
/3/i/j^ = 0.47-0.51) (Carlon et at, 2001) suspected that 
the transition (at least for low-D values) is PC type. 
However the lack of Z2 symmetry, parity conservation 
and further numerical data (Hinrichsen, 2001b; Odor, 
2000) exclude this possibility. Note, that the MF expo- 
nents are also different from those of the PC class. Sim- 
ulations and CAM calculations for the one dimensional 
n = 1 AF model (Odor, 2000, 2003b) corroborated the 
two new universality class prospect (see Fig. 23 and Ta- 
ble XXI) . The order parameter exponent (/3) seems to be 
very far from both of the DP and PC class values (Odor, 
2000, 2003b). 

The two distinct class behavior may be explained on 
the basis of competing diffusion strengths of particles and 
pairs (i.e. for large D-s the explicit diffusion of lonely 
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nontrivial exponent estimates, the latter one by Hinrich- 
sen provided a hypothesis that the ultimate long time 
behavior should be characterized by DP behavior. 

If we replace the annihilation process 2A ^ by co- 
agulation 2A ^ ^ in (209) we get the annihilation- 
coagulation model. GMF approximations and simula- 
tions of this model resulted in similar phase diagram than 
that of the PCPD model albeit without any sign of two 
distinct regions. In agreement with this CAM approxi- 
mations and simulations for the Id model found the same 
kind of continuous transition independently from D, with 
exponents in agreement with those of the PCPD in the 
low-D region (Odor, 2001a; Park and Kim, 2002). Again 
the spatial symmetry of particle production was found 
to be irrelevant. An exact solution was found by (Henkel 
and Hinrichsen, 2001) for the special case in Id, when 
the diffusion rate is equal to the coagulation rate, cor- 
responding to the inactive phase according to which the 
particle decay is like of ARW: p oc 



particles is stronger). Similar behavior was observed in 
case of Id models with coupled (conserved) diffusive field 
(see Sect. V.H). However a full agreement has not been 
achieved in the literature with respect the precise val- 
ues of the critical exponents. The low- 13 a is supported 
by (Park and Kim, 2002) who considered the case with 
coagulation and annihilation rates three times the diffu- 
sion rate. On the other hand the high-_D a of Table XXI 
coincides with that of (Kockelkoren and Chate, 2003a), 
who claim a single value for < 13 < 1. By assuming 
logarithmic corrections it was shown (Odor, 2003b) that 
a single universality class can be supported indeed with 
exponents 

a = 0.21(1), /? = 0.40(1), 
Z = 1.75(15), = 0.38(1) , (220) 

however there is no strong evidence for such corrections. 
Although the upper critical dimension is expected to be 
at dc = 2 (Odor et ai, 2002) one may not exclude the 
possibility of a second critical dimension (d^ = 1) or 
topological effects in Id that may cause logarithmic cor- 
rections to scaling. The spreading exponent r] seems to 
change continuously by varying D. Whether this is true 
asymptotically or the effect of some huge correction to 
scaling is still not clear. The simulations of (Odor, 2000) 
confirmed that it is irrelevant whether the particles pro- 
duction is spatially symmetric: Aid A AAA or spatially 
asymmetric: AAfb AAA, fbAA AAA. Recent sim- 
ulations and higher level GMF approximations suggest 
(Odor, 2003b; Odor et ai, 2002) that the pecuUarities 
of the pair approximation are not real; for N > 2 clus- 
ter approximations the low-D region scaling disappears. 
Recently two studies (Dickman and de Menenzes, 2002; 
Hinrichsen, 2002) reported non-universality in the dy- 
namical behavior of the PCPD. While the former one by 
Dickman and Menezes explored different sectors (a reac- 
tive and a diffusive one) in the time evolution and gave 



2. Cyclically coupled spreading with pair annihilation 

In this section I show an explicit two-component re- 
alization of the PCPD class. A cyclically coupled two- 
component reaction-diffusion system was introduced by 
(Hinrichsen, 2001c) 

A->2A, A ^ 0, A^ B, 2B ^ A, B% ^ %B 

(221) 

which mimics the PCPD model (Sect. V.F.I) by mapping 
pairs to As and single particles to B-s. This model is 
a coupled DP+ARW system. Its 1 -I- Id critical space- 
time evolution pattern looks very similar to that of the 
PCPD model. The appearance of the evolution in space- 
time seems to be a particular feature of this class. It is 
built up from compact domains with a cloud of lonely 
particles wandering and interacting with them. Further- 
more this model also has two non-symmetric absorbing 
states: a completely empty one and another with a sin- 
gle wandering B. By fixing the annihilation and diffusion 
rates of B-s {r = D — 1) the model exhibits continuous 
phase transition by varying the production rate of ^-s 
and the A —^ B transmutation rate. The simulations in 
Id showed that pA (x Pb for large times and resulted in 
the following critical exponent estimates 

a = 0.21(2), /3 = 0.38(6), /?' = 0.27(3), 

Z= 1.75(5), 1^11= 1.8(1), (222) 

satisfying the generalized hyperscaling relation (155). 
These exponents are similar to those of PCPD model 
in the high diffusion region (see Table XXI) that is rea- 
sonable since Z3 = 1 is fixed here. 

3. The parity conserving annihilation-fission model 

As we have seen parity conservation plays an important 
role in unary production systems. In case of BARW pro- 
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cesses it changes the universahty of the transition from 
DP (Sect.IV.A.3) to PC (Sect.IV.D.l) class. The ques- 
tion arises if one can see similar behavior in case of bi- 
nary production systems. Recently (Park et al., 2001) 
investigated a parity conserving representative (n = 2) 
of the Id AF model (209). By performing simulations 
for I0W--D-S they found critical exponents that are in the 
range of values determined for the corresponding PCPD 
class. They claim that the conservation law docs not af- 
fect the critical behavior and that the binary nature of 
the offspring production is sufficient condition for this 
class (see however Sect. V.F.2, where there is no such 
condition). 

The two dimensional version of the parity conserving 
AF model was investigated by GMF and simulation tech- 
niques (Odor et al, 2002). ' While the = 1,2 GMF 
approximations showed similar behavior as in case of 
the PCPD model (Sect. V.F.I) (including the two class 
prediction for N = 2) the A^ = 3, 4 approximations do 
not show D dependence of the critical behavior: /3 = 1, 
/32 = 2 was obtained for D > 0. Large scale simulations 
of the particle density confirmed the mean-field scaling 
behavior with logarithmic corrections. This result can be 
interpreted as a numerical evidence supporting that the 
upper critical dimension in this model is: dc = 2. The 
pair density decays in a similar way but with an addi- 
tional logarithmic factor to the order parameter. This 
kind of strongly coupled behavior at and above critical- 
ity was observed in case of the PCP model too (see Sect. 
V.E.I). At the D = endpoint of the transition line 
2+1 d class DP criticality (see Sect. IV. A) was found 
for p2 and for p — p{pc)- In the inactive phase for p{t) 
we can observe the two-dimensional ARW class scaling 
behavior (sec Sect.IV.C.l), while the pair density decays 
as p2 oc t~'^. Again like in d = 1 the parity conservation 
seems to be irrelevant. 



G. BARWe with coupled non-difTusive field class 

Similarly to the PCP (Sect.V.E.l) the effect of in- 
finitely many frozen absorbing states has been investi- 
gated in case of a BARWe model. A parity conserving 
version of the Id PCP model (Sect. V.E.I) was intro- 
duced by (Marques and Mendes, 1999), in which pairs fol- 
low a BARW2 process, while lonely particles are frozen. 
Simulations showed that while the critical behavior of 
pairs in case of homogeneous, random initial distribu- 
tion belongs to the PC class (Sect. IV.D.l), the spread- 
ing exponents satisfy hyperscaling (155) and change con- 
tinuously by varying the initial particle density. These 
results are similar to those found in the PCP model. 
Again long-memory effects arc responsible for the non- 
universal behavior in case of seed-like initial conditions. 
The slowly decaying memory was confirmed by studying 
a one-dimensional, interacting monomer-monomer model 
(Park and Park, 2001) by simulations. 



H. DP with diffusive, conserved slave field classes 

Motivated by model C results (Sect. Ill), it is an 
obvious question what happens with the phase transi- 
tion to an absorbing state of a reaction-diffusion sys- 
tem if a conserved secondary density is coupled to a 
non-conserved order parameter. One can deduce from 
the BARWl spreading process (Sect.IV.A.3) a two- 
component, reaction-diffusion model (DCF) (Kree et al., 
1989; Wijland et al, 1998) that exhibits total particle 
density conservation as follows: 



A + B-^2B, 



(223) 



By varying the initial particle density {p = pa(0)-I-/9b(0)) 
continuous phase transition occurs. General field theo- 
retical investigation was done by (Wijland et al., 1998) 
(for the equal diffusion case: Da = Db hy (Kree et al., 
1989)). The mean- field exponents that are valid above 
rfc = 4 are shown in Table XXII. The rescaled action of 
this model is 

S[ip, ip, ^, V^] = / d'^xdt lifidt -A)ip + lp{dt + X{a - A)) V 

-\-pTpA^ + gtp1p{tp — 1p) + Ull)tl){(fi + Tp) 

-pB{0)5{t) ^] (224) 

where tp and (f) are auxiliary fields, defined such that their 
average values coincide with the average density of B par- 
ticles and the total density of particles respectively. The 
coupling constants are related to the original parameters 
of the master equation by 



M = 1 - Db/Da g ■- 

vi = V2 = -U3 = k/DA u ■■ 

pb{o) = pb{Q)I^p 



ky^/DA \<T = fc(p^f - p)/Da 
-k^/DA X = Db/Da 

(225) 



If one omits from the action Eq. (224) the initial time 
term proportional to pb(0), then the remainder is, for 
/U = (i. e. Da = Db), invariant under the time reversal 
symmetry 



■iHx, t) 

tp{x,t) 
(p{x,t) 
Tp{x, t) 



-ip{x,-t) 
-ip{x,-t) 
Tp{x, -t) 
ip{x, -t) 



(226) 



Its epsilon expansion solution (Kree et al., 1989) and sim- 
ulation results (de Freitas et al., 2000; Fulco et al., 2001) 
are summarized in Table XXII. Interestingly RG pre- 
dicts Z = 2 and v± = 2/d \n all orders of perturbation 
theory. 

The breaking of this symmetry for p ^ 0, that is, when 
the diffusion constants Da and Db are different causes 

different critical behavior for this system. For Da < Dg 
RG (Wijland et al, 1998) predicts new classes with Z = 
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d 


P 


Z 




1 

4 - e 


0.44(1) 

1 - e/8 


2 
2 


2.1(1) 

2/d 



TABLE XXII Summary of results for DCF classes for Da 
Db 



d 


P 


Z 




1 

4-e 


0.33(2) 
1 


2 


2 
2/d 



TABLE XXIII Summary of results for DCF classes for Da < 
Db 



2, P = 1, i/± = 2/d, but simulations in Id (Fulco et ai, 
2001) show different behavior (see Table XXIII). For 
non-poissonian initial particle density distributions the 
critical initial slip exponent 77 varies continuously with 
the width of the distribution of the conserved density. 
The Da = extreme case is discussed in Section V.I. 

For Da > Db no stable fixed point solution was found 
by RG hence (Wijland et ai, 1998) conjectured first or- 
der transition for which signatures were found in 2d by 
simulations (Oerding et ai, 2000). However e expansion 
may break down in case of the occurrence of another 
critical dimension d'^ < dc = 4: for which simulations in 
Id (Fulco et ai, 2001) provided numerical support (see 
Table XXIV). 



I. DP with frozen, conserved slave field classes 

If the conserved field coupled to the BARWl pro- 
cess (eq. 223) is non-diffusive non-DP universality class 
(NDCF) behavior is reported again (Pastor-Satorras and 
Vespagnani, 2000; Rossi et a/., 2000). The corresponding 
action can be derived from (224) in the Da = limit: 

S= J d'^xdtllfiidt + r - DV^)(p + i^{dt - XV^))i^ 
+giptp{tp - ip) + uip'4){ip + 7^) +vi {iptjj) ^ 

+V2iptp{'4>ip — tlxfi) + vaifTpip'tjj (227) 

By neglecting irrelevant terms (227) is invariant under 
the shift transformation 



ijj ^ + A, r ^ r — V2A 



(228) 
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4-e 


0.67(1) 
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TABLE XXIV Summary of results for DCF classes for Da > 
Db 
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Z 


^^11 


a 


V 


1 


0.14(1) 


0.28(1) 




1.5(1) 


2.5 






2 


0.50(5) 


0.64(1) 


1.59(3) 


1.55(4) 


1.29(8) 


2.22(3) 


0.29(5) 


3 


0.90(3) 


0.84(2) 


1.23(4) 


1.75(5) 


1.12(8) 


2.0(4) 


0.16(5) 


4 


1 


1 


1 


2 


1 


2 






TABLE XXV Summary of results for NDCF classes 



where A is any constant. The field theoretical analysis of 

this action has run into difficulties (Pastor-Satorras and 
Vespagnani, 2000). The main examples for the NDCF 
classes are the conserved threshold transfer process and 
the conserved reaction-diffusion model (Pastor-Satorras 
and Vespagnani, 2000; Rossi et ai, 2000). Furthermore 
the models described by the NDCF classes embrace a 
large group of stochastic sandpile models (Jensen, 1998) 
in particular fixed-energy Manna models (Dhar, 1999; 
Dickman et ai, 1998; Manna, 1991; Munoz et ai, 2001). 
The upper critical dimension = 4 was confirmed by 
simulations (Liibeck and Hucht, 2002). 

It was also shown (Alava and Munoz, 2001) that these 
classes describe the depinning transition of quenched 
Edwards- Wilkinson (see Sect. VI. C) or linear inter- 
face models (LIM) (Barabasi and Stanley, 1995; Halpin- 
Healy and Zhang, 1995) owing to the fact that quenched 
disorder can be mapped onto long-range temporal cor- 
relations in the activity field (Marsili, 1994). However 
this mapping could not be done on the level of Langevin 
equations of the representatives of NDCF and LIM mod- 
els and in Id this equivalence may break down (Alava 
and Munoz, 2001; Dickman et al., 2001a; Kockelkoren 
and Chatc, 2003b). The critical exponents determined 
by simulations (Dickman et ai, 2002; Liibeck, 2001, 2002; 
Pastor-Satorras and Vespagnani, 2000; Rossi et al., 2000) 
and GMF-hCAM method in Id (Dickman, 2002) are sum- 
marized in Table XXV. Similarly to the PCP these 
models exhibit infinitely many absorbing states, there- 
fore non-universal spreading exponents are expected (in 
Table XXV the exponent 77 corresponding to natural ini- 
tial conditions is shown). 



J. Coupled N-component DP classes 

From the basic reaction-diffusion systems one can gen- 
erate N-component ones coupled by interactions sym- 
metrically or asymmetrically. In (Janssen, 1997b, 2001) 
Jansscn introduced and analyzed by field theoretical RG 
method (up to two loop order) the quadratically coupled, 
N-species generalization of the DP process of the form: 



Act 

Aa + Ap 



2A„ 



kAa + 1 An, 



(229) 



where k,l may take the values (0,1). He has shown 
that the multi-critical behavior is always described by 
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FIG. 24 Flow of the interspecies couplings in the two- 
component, DP model under renormalization. D means de- 
coupled, S symmetric, U unidirectional fixed points (Janssen, 
2001). 



the Reggeon field theory (DP class), but this is unstable 
and leads to unidirectionally coupled DP systems (see 
Fig.24). 

He has also shown that by this model the linearly, uni- 
directionally coupled DP (UCDP) (see Sect.VI.H) case 
can be described. The universality class behavior of 
UCDP is discussed in Sect.VI.H. 

In one dimension, if BARWo type of processes are 
coupled (which alone exhibit DP class transition (sec 
Sect. IV. A. 3) hard-core interactions can modify the 
phase transition universality (see Sect. V.L.I). 



K. Coupled N-component BARW2 classes 

Bosonic, A''-component BARW systems with two off- 
springs (N-BARW2), of the form 



Aa + 2AfS 



(230) 

(231) 
(232) 



were introduced and investigated by (Cardy and Taubcr, 
1998) via field theoretical RG method. These models ex- 
hibit parity conservation of each species and permutation 
symmetry on types (generalization for 0{N) symme- 
try is violated by the annihilation term). The A —^ 3 A 
processes turns out to be irrelevant, because like pairs an- 
nihilate immediately. Models with (231) branching terms 
exhibit continuous phase transitions at zero branching 
rate. The universality class is expected to be indepen- 
dent from A'' and coincides with that of the N ^ oo 
(N-BARW2) model that could be solved exactly. The 
critical dimension is dc = 2 and for d < 2 the exponents 
are 



while at d = dc = 2 logarithmic correction to the density 
decay is expected. Simulations on a d = 2 (fermionic) 
lattice model confirmed these results (Odor, 2001c). 

In one dimension it turns out that hard-core inter- 
actions can be relevant and different universal behav- 
ior emerges for fermionic models (see Sect. V.L). The 
bosonic description of fermionic models in one dimension 
works (at least for static exponents) in case of pairwise 
initial conditions (see Sect. V.B) when different types 
of particles do not make up blockades for each other. 
Such situation happens when these particles are gener- 
ated as domain walls of N+1 component systems exhibit- 
ing Sj\r+i symmetric absorbing states (see Sects. V.K.I, 
IV.D.2). 



1. Generalized contact processes with n > 2 
absorbing states in Id 

The generalized contact process has already been in- 
troduced in Sect. IV. D. 3 with the main purpose to show 
an example for PC class universality class transition in 
case of Z2 symmetric absorbing states. The more general 
case with n > 2 permutation symmetric absorbing states 
was investigated using DMRG method by (Hooyberghs 
et ai, 2001) and turned out to exhibit N-BARW2 tran- 
sition. In the one dimensional model, where each lattice 
site can be occupied by at most one particle (A) or can 
be in any of n inactive states (0i, 02 . . • 0^) the reactions 
are: 



AA^AiDk,hA 

A^Dk.dikA^Mk 
A(Dk,i/}kA^ AA 
%kh ^ A%u 0feA 



13=1, Z = 2, a = d/2, vn=2/d, v^ = l/d 



(233) 



with rate A/n (234) 

with rate /Xk (235) 

with rate 1 (236) 

(fc 7^ I) with rate 1 (237) 

The original contact process (Sect. IV.A.l), corresponds 
to the n = 1 case, in which the reaction (237) is 
obviously absent. The reaction (237) in the case n > 2 
ensures that configurations as (0i0i . . . ^SSj^j . . . ^j^j), 
with i ^ j are not absorbing. Such configurations do 
evolve in time until the different domains coarsen and 
one of the n absorbing states (0i0i . . . 0i), (0202 ■ • • 02), 
• • • (0ri0n • ■ • 0n) IS reached. For generalized contact 
processes with n = 2, simulations (Hinrichsen, 1997) 
and a DMRG study (Hooyberghs et ai, 2001) proved 
that the transition falls in the PC class if /ii = /i2, while 
if the symmetry between the two absorbing states was 
broken (/xi ^ ^2) a DP transition was recovered. 

The DMRG study for n = 3 and n = 4 showed that, 
the model is in the active phase in the whole parameter 
space and the critical point is shifted to the limit of in- 
finite reaction rates. In this limit the dynamics of the 
model can be mapped onto the zero temperature n-state 
Potts model (see also the simulation results of (Lipowski 
and Droz, 2002a)). It was conjectured by (Hooyberghs 
et ai, 2001) the model is in the same N-BARW2 uni- 
versality class for all n > 3. By calling a domain wall 



between %i and 0j as Xij one can follow the dynamics 
of such variables. In the limit A ^ oo, X^j coincides 
with the particle A and in the limit /i ^ oo Xij coin- 
cides with the bond variable For finite values of 
these parameters one still can apply this reasoning at a 
coarse-grained level. In this case X^ is not a sharp do- 
main wall, but an object with a fluctuating thickness. 
For n = 2 it was shown in Section IV. D. 3 that such do- 
main wall variables follow BARW2 dynamics (175). For 
n > 2 one can show that besides the N-BARW2 reactions 
(230), (231) (involving two types of particles maximum) 
reaction types involving three different domains {i ^ j, 
i^k and j ^ k): 

Xij XikXkj XikXkj — » Xij (238) 

occur with increasing importance as n ^ oo. These re- 
actions break the parity conservation of the N-BARW2 
process, therefore the numerical findings in (Hooyberghs 
et ai, 2001) for n = 3, 4 indicate that they are proba- 
bly irrelevant or the conditions for N-BARW2 univer- 
sal behavior could be relaxed. Owing to the fact that the 
Xij variables are domain walls they appear in pairwise 
manner, hence hard-core exclusion effects are ineffec- 
tive for the critical behavior i Id. To see the effect of 
pairwise initial conditions for dynamical exponents see 
Section V.B. 

For n = 3 upon breaking the global symmetry 
to a lower one, one gets a transition either in the di- 
rected percolation (Sect. IV. A), or in the parity conserv- 
ing class (Sect. IV. D), depending on the choice of param- 
eters (Hooyberghs et al., 2001). Simulations indicate 
(Lipowski and Droz, 2002b) that for this model local 
symmetry breaking may also generate PC class transi- 
tion). 

L. Hard-core 2-BARW2 classes in Id 

Besides the effects of coupling interactions in low di- 
mensions blockades generated by hard-core particles may 
also play an important role. The effect of particle exclu- 
sion (i.e. AB ^ BA) in 2-BARW2 models (Sect.V.K) 

was investigated in (Kwon et al., 2000; Odor, 2001c). 
For d = 2 the bosonic field theoretical predictions (Cardy 
and Tauber, 1998) were confirmed (Odor, 2001c) (mean- 
field class transition with logarithmic; corrections). In 
one dimension however two types of phase transitions 
were identified at zero branching rate {a = 0) depending 
on the arrangement of offspring relative to the parent in 
process (231). Namely if the parent separates the two 
offsprings (2-BARW2s): 

A A BAB (239) 

the steady state density is higher than in the case when 
they are created on the same site (2-BARW2a): 

A ^ ABB (240) 
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exponent 


N-BARW2 


N-BARW2S 


N-BARW2a 




2 


2.0(1)10.915(2) 


8.0(4)13.66(2) 


z 


2 


4.0(2)11.82(2)* 


4.0(2)11.82(2)* 


a 


1/2 


0.25(1)10.55(1)* 


0.25(1)10.55(1)* 


(3 


1 


0.50(1)11.00(1) 


2.0(1)11.00(1) 



TABLE XXVI Summary of critical exponents in one dimen- 
sion for N-BARW2 like models. The N-BARW2 data are 
quoted from (Cardy and Tauber, 1996). Data divided by 
" I" correspond to random vs. pairwise initial condition cases 
(Hooyberghs et al, 2001; Mcnyhard and Odor, 2002; Odor 
and Menyhard, 2000). Exponents denoted by * exhibit slight 
initial density dependence. 



at a given branching rate, because in the former case 
they are unable to annihilate each other. This results in 
different order parameter exponents for the symmetric 
(2-BARW2s) and for the asymmetric (2-BARW2a) cases 

I3s = 1/2, /3„ = 2. (241) 

This result is in contrast with a widespread belief that 
the bosonic field theory (where AB ^ BA is allowed) 
can describe these systems (because in that case the crit- 
ical behavior is different (Sect. 233)). This observation 
led (Kwon et ai, 2000) to the conjecture that in one- 
dimensional, reaction-diffusion systems a series of new 
universality classes should appear if particle exclusion is 
present. Note however, that since the transition is at 
(T = in both cases the on-critical exponents do not de- 
pend on how particles are created and they can be identi- 
fied with those described in Sect. (V.B). In (Odor, 2001c) 
a set of critical exponents satisfying scaling relations have 
been determined for this two new classes shown in Table 
XXVI. 



1. Hard-core 2-BARWo models in Id 

Hard-core interactions in the two-component, one- 
offspring production model (2-BARWl) were investi- 
gated in (Odor, 2001d). Without interaction between 
different species one would expect DP class transition. 
By introducing the AB <A BA blocking to the two- 
component model: 

A-^AA B-^BB (242) 
AA^-^9 BB^-^% (243) 

a DP class transition at cr = 0.81107 was located. Note 
that the effect exerted by different species on each other 
is irrelevant now unlike for the case of coupled ARW 
(Sect. V.B). On the other hand if we couple the two sub- 
systems by production: 

A^AB A^BA (244) 
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B^AB B^BA (245) 
AA^-^$ BB^-^% (246) 

a continuous phase transition emerges at cr = rate - 
therefore the on-critical exponents are the same as those 
described in Sect. (V.B) - and the order parameter expo- 
nent was found to be j3 = 1/2. Therefore this transition 
belongs to the same class as the 2-BARW2s model (see 
Sect.V.L). The parity conservation law, which is relevant 
in case of one component BARW systems (PC versus DP 
class) turns out to be irrelevant here. This finding re- 
duces the expectations suggested by (Kwon et at, 2000) 
for a whole new series of universality classes in Id sys- 
tems with exclusions. In fact the blockades introduced 
by exclusions generate robust classes. In (Odor, 2001d) 
a hypothesis was set up that in coupled branching 
and annihilating random walk systems of N-types 
of excluding particles for continuous transitions 
at cr = two universality classes exist, those of 
2-BARW2S and 2-BARW2a models, depending on 
whether the reactants can immediately annihilate (i.e. 
when similar particles are not separated by other type(s) 
of particle(s)) or not. Recent investigations in similar 
models (Lipowski and Droz, 2001; Park and Park, 2001) 
are in agreement with this hypothesis. 

2. Coupled binary spreading processes 

Two-component versions of the PCPD model (Sect. 
V.F.I) with particle exclusion in Id were introduced and 
investigated by simulations in (Odor, 2002) with the aim 
to test if the hypothesis for N-component BARW sys- 
tems set up in (Odor, 2G01d) (Sect.V.L. 1) can be apphed 
for such models. The following models with the same 
diffusion and annihilation terms {AA 0, BB 0) 
as in (Sect. V.B) and different production processes were 
investigated. 

1) Production and annihilation random walk model (2- 
PARW): 

AA ^ AAB, AA ^ BAA, (247) 
BB ^ BBA, BB ^ ABB . (248) 

2) Symmetric production and annihilation random 
walk model (2-PARWS): 

AA AAA, (249) 
BB BBB . (250) 

These two models exhibit active steady states for a > 
with a continuous phase transition at CTc = 0. There- 
fore the exponents at the critical point are those of the 
2-component ARW model with exclusion (Sect. V.B). 
Together with the exponent P — 2 result for both cases 
this indicates that they belong to the N-BARW2a class. 
This also means that the hypothesis set up for N-BARW2 






time (MCS) 



FIG. 25 Space-time evolution of the 2-PARWAS model at 
the critical point (Odor, 2002). Black dots correspond to A 
particles, others to B-s. 

systems (Odor, 2001d) (Sect.V.L. 1) can be extended. 

3) Asymmetric production and annihilation random 
walk model (2-PARWA): 

AB ^ ABB, AB ^ AAB, (251) 
BA ^ BAA, BA ^ BBA . (252) 

This model does not have an active steady state. The 
AA and BB pairs annihilate themselves on contact, 
while if an A and B particle meet an AB ABB A 
process eliminates blockades and the densities decay 
with the p oc law for a > 0. For a = the 

blockades persist, and in case of random initial state 
p oc decay (see Sect. V.B) can be observed here. 

4) Asymmetric production and annihilation ran- 
dom walk model with spatially symmetric creation (2- 
PARWAS): 

AB ^ ABA, AB ^ BAB, (253) 
BA ^ BAB, BA ^ ABA . (254) 

In this case AB blockades proliferate by production 
events. As a consequence of this an active steady state 
appears for a > 0.3253(1) with a continuous phase tran- 
sition. The space-time evolution from random initial 
state shows (Fig. 25) that compact domains of alternating 
..ABAB.. sequences separated by lonely wandering parti- 
cles are formed. This pattern is very similar to what was 
seen in case of one-component binary spreading processes 
(Hinrichsen, 2GGlc): compact domains within a cloud of 
lonely random walkers, except that now domains are built 
up from alternating sequences only. This means that the 
..AAAA... and ...BBBB ... domains decay by this anni- 
hilation rate and the particle blocking is responsible for 
the formation of compact clusters. In the language of 
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coupled DP + ARW model (Hinrichsen, 2001c) the pairs 
following DP process are the AB pairs now, which cannot 
decay spontaneously but through an annihilation process: 
AB + BA 0. They interact with two types of particles 
executing annihilating random walk with exclusion. The 
simulations resulted in the critical exponent estimates: 
/3 = 0.37(2), a = 0.19(1) and Z = 1.81(2) which agree 
fairly well with those of the PCPD model in the high 
diffusion rate region (Odor, 2000). 



i it t t f i i t i i 1 t t I ^Pi"s 
.'4*. . . ,''4\ kinks 



/\ surface 
< + / ^ 



VI. INTERFACE GROWTH CLASSES 



FIG. 26 Mapping between spins, kinks and surfaces 



Interface growth classes are strongly related to the ba- 
sic universality classes discussed so far and can be ob- 
served in experiments more easily. For example one 
of the few experimental realizations of the robust DP 
class (IV.A) is related to a dcpinning transition in in- 
homogeneous porous media (Buldyrev et al, 1992) (see 
Sect. VI. F). The interface models can either be defined 
by continuum equations or by lattice models of solid-on- 
solid (SOS) or restricted solid-on-solid (RSOS) types. In 
the latter case the height variables hi of adjacent sites 
are restricted 



\hi - hi+i\ < 1 



(255) 



The morphology of a growing interface is usually charac- 
terized by its width 



1/2 



(256) 



In the absence of any characteristic length, growth pro- 
cesses are expected to show power-law behavior of the 
correlation functions in space and height and the Family- 
Vicsek scaling (Family and Vicsek, 1985) form: 

iy(i,i) = i«/^/(L/e||(t)), (257) 
describes the surface, with the scaling function f{u) 



fin) 



if u < 1 
const, if M 1 



(258) 



Here a is the roughness exponent and characterizes the 
stationary regime in which correlation length ^| | (t) has 
reached a value larger than the system size L. The ratio 
(3 = a/Z is called as the growth exponent and character- 
izes the short time behavior of the surface. Similarly to 
equilibrium critical phenomena, these exponents do not 
depend on the microscopic details of the system under 
investigation. Using these exponents it is possible to di- 
vide growth processes into universality classes (Barabasi 
and Stanley, 1995; Krug, 1997). The (257) scaling form 
of W'^ is invariant under A the rescaling 

Ax, A^t, h{x, t) A-"h{x, t) (259) 



Recently anomalous roughening has been observed in 
many growth models and experiments. In these cases the 
measurable aioc roughness exponent is different from a 
and may satisfy different scaling law (Brii and al., 1998; 
Dasgupta et al., 1996; H. Yang and Lu, 1994; Jeffries 
et al, 1996; Krug, 1997; Lopez and Rodriguez, 1996, 
1997; Lopez and Schmittbuhl, 1998; Morel and al., 1998; 
Olivcira, 1992; Sarma and al., 1996; Schroeder and al., 
1993). 

Surfaces in 1 dimensional systems can be mapped 

onto a time step of a d dimensional particle reaction- 
diffusion or spin models. For example the Id Kawasaki 
spin model corresponding to the K 3A' proc;c;ss with 
random walk of kinks the mapping onto the Id surface is 
shown in Fig. 26. 



This means that to a unique {sj} spin configuration at 
time t corresponds a spatial profile {hiit)] by accumu- 
lating the spin values 



^'(*) =E*J 



(260) 



In Id the surface can also be considered as a random 
walker with fluctuation 



(261) 



hence the roughness exponent is related to the dynamical 
exponent Z^j, as 



a = l/Zu 



(262) 



The a = 1/2 corresponds to uncorrelated (or finite corre- 
lation length) random walks. If a > 1/2 the surface ex- 
hibits correlations, while if a < 1/2 the displacements in 
the profile are anti-correlated. Since the surfaces may ex- 
hibit drifts the fluctuations around the mean is measured 
defining the local roughness (Hurst) exponent. Using this 
surface mapping (Sales et al., 1997) have characterized 
the different classes of Wolfram's Id cellular automata 
(Wolfram, 1983). 

One can show that by coarse graining the Id Kawasaki 
dynamics 



Wi 



1 
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[1 - CTiCTi+i + A((7i+i - CTj)] 



(263) 
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a mapping can be done onto the KPZ equation (273), 
and the surface dynamics for A 7^ (corresponding to 
anisotropic case) is in the KPZ class (Sect. VI. D) while for 
A = it is in the Edward Wilkinson class (Sect. VLB). 
While these classes are related to simple random walk 
with Zw = 1/a = 2 the question arises what surfaces are 
related to other kind of random walks (example Levy- 
flights or correlated random walks etc ...). Recently it 
was shown that globally constrained random walks (i.e. 
when a walker needs to visit each site an even number of 
times) can be mapped onto surfaces with dimer-type of 
dynamics (Noh et ai, 2001) with Z^^ = 3 = 1/a. 

By studying correspondence between lattice models 
with absorbing states and models of pinned interfaces in 
random media (Dickman and Munoz, 2000) established 
the scaling relation 

P=l-l3/uii (264) 

that was confirmed numerically for d = 1,2,3,4 contact 
processes (Sect. IV. A. 1). The local roughness exponent 
was found to be smaller than the global value indicating 
anomalous surface growth in DP class models. 

In interface models different types of transitions may 
take place. Roughening transitions may occur be- 
tween smooth phase characterized by finite width W 
(in an infinite system) and rough phase when the width 
and diverges in an infinite system (but saturates in fi- 
nite ones) by varying some control parameters (e). Near 
the transition point the spatial (^||) and growth direction 
correlations (^_l) diverge as 

C|| oc e^'ii (265) 
C_L OC e"^ (266) 

(note that in RD systems ^|| denote temporal correlation 
length). In the smooth phase the heights hi{t) are corre- 
lated below While in equilibrium models roughening 
transitions exist in d > 1 dimensions only in nonequilib- 
rium models this may occur in d = 1 as well. 

An other surface transition is the depinning transi- 
tion, when as the consequence of changing some control 
parameter (usually an external force F) the surface starts 
propagating with speed v and evolves in a rough state. 
Close to the transition v is expected to scale as 

v(xiF- F^f (267) 

with the 9 velocity exponent and the correlation expo- 
nents diverge. Known depinning transitions (in random 
media) are related to absorbing phase transitions with 
conserved quantities (see Sects. VI. C, VI. E). 

In the rough phase a so-called faceting phase tran- 
sition may also take place when up-down symmetrical 
facets appear. In this case the surface scaling behavior 
changes (see Section VI.I) . 



A. The random deposition class 

The random deposition is the simplest surface growth 

process that involves uncorrelated adsorption of particles 
on top of each other. Therefore columns grow indepen- 
dently, linearly without bounds. The roughness exponent 
a (and correspondingly Z) is not defined here. The width 
of the surface grows as cx i^/"^ hence /3 = 1/2 in all 
dimensions. An example for such behavior is shown in a 
dimer growth model in Sect. VI.I.l. 



B. Edwards-Will<inson (EW) classes 

As it was mentioned in Sect. VI growth models of 
this class can easily mapped onto spins with sym- 
metric Kawasaki dynamics or to particles with ARW 
(Sect.IV.C.l). If we postulate the following translation 
and reflection symmetries 

x^x-l-Ax t^t+At h^h+Ah x ^ -x h^-h 

(268) 

we are led to the Edwards- Wilkinson (EW) equation (Ed- 
wards and Wilkinson, 1982) 

at/i(x, t)=v + aV'^hix, t) + C(x, t), (269) 

which is the simplest stochastic differential equation that 
describes a surface growth with these symmetries. Here 

V denotes the mean growth velocity, a the surface tension 
and C, the zero-average Gaussian noise field with variance 

< C(x,i)C(x',0 >= 2D6'^-\yi-x'){t-t') (270) 

This equation is linear and exactly solvable. The critical 
exponents of EW classes are 

^ = (^-^), Z = 2 (271) 



C. Quenched EW classes 

In random media linear interface growth is described 
by the so called quenched Edwards- Wilkinson equation 

dthi-K, t) = aV^h{x, t)+F + r/(x, /i(x, t)), (272) 

where F is a constant, external driving term and 
r/(x, ft,(x, t)) is the quenched noise. The corresponding 
linear interface models (LIM) exhibit a depinning transi- 
tion at Fc- The universal behavior of these models were 
investigated in (Kim et at, 2001; Narayan and Fisher, 
1993; Nattermann et ai, 1992) and it was shown to be 
equivalent with NDCF classes (Sect. V.I). Analytical 
studies (Narayan and Fisher, 1993) predict a = (4 — rf)/3 
andZ = 2-(2/9)(4-(i). 



d 


a 




Z 


1 


1/2 


1/3 


3/2 


2 


0.38 


0.24 


1.58 


3 


0.30 


0.18 


1.66 



TABLE XXVII Scaling exponents of KPZ from (Barabasi 
and Stanley, 1995) 



D. Kardar-Parisi-Zhang (KPZ) classes 

If we drop the h — > —h symmetry from (268) we can 
add a term to the (269) equation that is the most rele- 
vant term in renormahzation sense breaking the up-down 
symmetry: 

dth{x,t) =w + (7V^/t(x,t) + A(V/i(x,t))2+C(x,f) (273) 

that is the Kardar-Parisi-Zhang (KPZ) equation (Kardar 
et al., 1986). Here again v denotes the mean growth 
velocity, a the surface tension and C, the zero-average 
Gaussian noise field with variance 

< C(x, i)C(x', t') >= 2£)5'^-i(x - x')(t - t') (274) 

This is non- linear, but exhibits a tilting symmetry as the 
result of the Galilean invariance of (273): 

/i ^ /i' + ex, X ^ x' - \et, t^t' (275) 

where e is an infinitesimal angle. As a consequence the 
scaling relation 

a + Z = 2 (276) 

holds in any dimensions. In one dimension the critical 
exponents are known exactly whereas for > 1 dimen- 
sions numerical estimates exist (Barabasi and Stanley, 
1995) (see Table XXVII). The upper critical dimension 
of this model is debated. Mode coupling theories and var- 
ious phenomenological field theoretical schemes (Halpin- 
Healy, 1990; Lassig, 1995; Lassig and Kinzelbach, 1997) 
settle to dc = 4. Contrary to analytical approaches nu- 
merical solution of the KPZ equation (Moser et al., 1991), 
simulations (Ala-Nissila et al., 1993; Kim and Koster- 
litz, 1989; Marinari et al, 2000; Tang et al, 1992; Wolf 
and Kertesz, 1987), and the results of real-space renor- 
mahzation group calculations (Castellano et al., 1998a, 
1999, 1998b) provide no evidence for a finite dc- Fur- 
thermore, the only numerical study (Tu, 1994) of the 
mode-coupling equations gives no indication for the ex- 
istence of a finite dc either. Recently simulations of the 
restricted solid-on-solid growth models were used to build 
the width distributions of d = 2 — 5 dimensional KPZ in- 
terfaces. The universal scaling function associated with 
the steady-state width distribution was found to change 
smoothly as d is increased, thus strongly suggesting that 
d = 4 is not an upper critical dimension for the KPZ 
equation. The dimensional trends observed in the scal- 
ing functions indicate that the upper critical dimension 
is at infinity (Marinari et al., 2002). 
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1. Multiplicative noise systems 

In the hope of classifying nonequilibrium phase transi- 
tion classes according to their noise terms (G. Grinstein 
and Tu, 1996; Tu et al., 1997) introduced and studied 
systems via the Langevin equation 

dtn{x, t) = D V^n(a;, t)—r n{x, t)—u n{x, t)'^+n{x, t) r]{x, t) , 

(277) 

exhibiting real multiplicative noise (MN) : 

(7?(x,t)) = , {ri{x,t)ri{x' ,t')) = 2v 5'^{x - x')5{t - t') , 

(278) 

that is proportional to the field (n). Since for DP class 
the noise is square root proportional to the field; for 
ARW and PC systems the noise is imaginary one may 
expect distinct universal behavior for such systems. Ref. 
(Howard and Tauber, 1997) on the other hand argued on 
field theoretical basis that such 'naive' Langevin equa- 
tions fail to accurately describe systems controlled by 
particle pair reaction processes, where the noise is in 
fact 'imaginary' and one should derive a proper Langevin 
equation by starting from master equation of particles. 
Therefore they investigated the simplest RD systems 
where both real and imaginary noise is present and com- 
pete: (a) 2A -^%,2A^ 2B, 2B 2A, and 25 ^ (b) 
2A^ % and 2A ^ {n + 2)A (see Sect. V.F) by setting 
up the action first - derived from the master equation 
of particles. In neither case were they able to recover 
the MN critical behavior reported in (G. Grinstein and 
Tu, 1996; Tu et al, 1997). Therefore they suspected that 
there might not be real RD system possessing the MN 
behavior. 

On the other hand (G. Grinstein and Tu, 1996; Tu 

et al., 1997) have established connection of MN systems 
via the Cole-Hopf transformation: n{x,t) = e'*^^'*) to 
the KPZ theory. They have shown in Id that the phase 
diagram and the critical exponents Z , vi^ and /? of the 
two systems agree within numerical accuracy. They have 
found diverging susceptibility (with continuously chang- 
ing exponent as the function of r) for the entire range of 
r. 



E. Quenched KPZ classes 

In random media nonlinear interface growth is de- 
scribed by the so called quenched KPZ equation 

(Barabasi and Stanley, 1995), 

dth{yi, t) = tTV2/i(x, t) + A(V/i(x, t)f +F + r]{x, /i(x, t)) 

(279) 

where _F is a constant, external driving term and 
r]{x,h{yi,t)) is the quenched noise (do not fluctuate in 
time). Its universal behavior was investigated in (Bul- 
drcv et al., 1993; Leschorn, 1996) and predicted a ~ 0.63 
in one dimension, a ~ 0.48 in two dimensions and 
a ~ 0.38 in three dimensions. It was shown numerically 
that in Id this class is described by 14-1 d directed perco- 
lation depinning (Tang and Leschhorn, 1992). In higher 
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Langovin equation 


a 


/3 


z 


dth = -KV*h + C 


4-d 


4-d 


4 


dth = i/V^/i + Cd 


-d 


-d 
4 


2 


dth = -KV'h + Cd 


'2-d 
2 


2-d 
X 


4 


dth = -KV^h + AiV2(V/i)2 + C 


4-d 
3 


4-d 
8+d 


8+d 
3 






d < 1 


dth = -KV^h + XiV^iVhf + Cd 


2-d 
3 


2-d 
10+d 


10+d 
3 






d > 1 




2-d 
2 


2-d 
8 


4 



TABLE XXVIII Summary of continuum growth classes dis- 
cussed in this section following ref. (Barabasi and Stanley, 
1995). 



dimensions however it is related to percolating directed 
surfaces (Barabasi et al., 1996). 

F. Other continuum growth classes 

For continuum growth models exhibiting the symme- 
tries 

X ^ X + Ax t+At h h+Ah x ^ -x (280) 

the possible general Langevin equations with relevant 
terms were classified as follows (Barabasi and Stanley, 
1995). 

• The deterministic part describes conservative or 
nonconservative process (i.e. the integral over the 

entire system is zero or not). Conservative terms 
are V^/i, V/i and \/^{\/hf. The only relevant 
nonconservative terms is the (V/i)^. 

• The system is linear or not. 

• The noise term is conservative (i.e. the result of 
some surface diffusion) with correlator 

(Cd(x, i)Cd(x', t/)) = {-2DaV^ + D^V^)5(x - x')<5(t - 1') 

(281) 

or nonconservative like in eq.(270) (as the result of 
adsorption, desorption mechanisms). 

Analyzing the surface growth properties of such systems 
besides the EW and KPZ classes five other universality 
classes were identified (see Table XXVIII). 

G. Classes of mass adsorption-desorption, aggregation 
and chipping models 

Based on the interest in sclf-organizcd critical systems 
in which different physical quantities exhibit power law 
distributions in the steady state over a wide region of 
the parameter space (Bak et al, 1987) a family of lattice 
models, in which masses diffuse, aggregate on contact. 



and also chip off a single unit mass was introduced by 
(Majumdar et ai, 1998, 2000a). Self-organized critical- 
ity has been studied in a variety of model systems ranging 
from sandpiles to earthquakes. A particularly simple lat- 
tice model due to Takayasu, with mass diffusion, aggre- 
gation upon contact and adsorption of unit masses from 
outside at a constant rate, exhibits self-organized critical- 
ity (Takayasu et ai, 1988; Takayasu and Takayasu, 1997): 
the steady-state mass distribution has a nontrivial power 
law decay for large masses in all dimensions (Takayasu 
et al, 1988). These mass adsorption-desorption models 
in one dimension arc defined as follows. A site i is chosen 
randomly and then one of the following events can occur: 

1. Adsorption: With rate q, a single particle is ad- 
sorbed at site i; thus rrn ^ rrn + l. 

2. Desorption: With rate p, a single particle is des- 
orbed from site i; thus mi mi — 1 provided 

TTii > 1. 

3. Diffusion and Aggregation: With rate 1, the mass 
mi at site i moves to a nearest neighbor site [either 
(i — 1) or (i + 1)] chosen at random. If it moves to 
a site which already has some particles, then the 
total mass just adds up; thus mj — > and mi±i — > 
mi±i + mi. 

4. Chipping: With rate w a bit of mass at the site 
"chips" off, e. provided rrij > 1 a single particle 
leaves site i 

While the Takayasu model (p = 0, w = 0) does not have 

a phase transition in the steady state, by introducing a 
nonzero desorption rate p induces a critical line pdo) in 
the p — q plane. For fixed q, if one increases p from 0, one 
finds that for all p < Pc{q), the steady state mass distri- 
bution has the same large m behavior as in the Takayasu 
case, i.e., 

P(m) oc m"^ (282) 

where the exponent r is the Takayasu exponent and is 
independent of q. For p = Pc{q), we find the steady state 
mass distribution still decays algebraically for large m, 
but with a new critical exponent Tc which is bigger than 
the Takayasu exponent Tt- For p > Pc{q), we find that 

P(m) ~ exp(-m/m*) (283) 

for large m where m* is a characteristic mass that di- 
verges if one approaches Pc{q) from the p > pdq) side. 
The critical exponent Tc is the same everywhere on the 
critical line pdq)- This phase transition occurs in all 
spatial dimensions including d = 1. The t exponent was 
determined for the mean-field and one dimensional cases 
((Majumdar et al, 2000b; Takayasu et al, 1988)) 

Tf^=3/2, Tf-^ = 5/2, tI'^ = 4/3, t,1'^ = 1.833, 

(284) 

although the location of dc is not known. 
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This model can also be mapped onto an interface dy- 
namics, if we interpret the configuration of masses as an 
interface profile regarding local height variable. 

The phase transition of the model can be qualitatively 
interpreted as a nonequilibrium wetting transition of the 
interface. The corresponding surface growth exponents 
are (Majumdar et al., 2000b) 

I^MF ^ jyg^ ^ ^ 2, ^^'^ = 0.358, Z = 2 . (285) 

In the p = q = conserved model (CM) a nonequi- 
librium phase transition occurs by varying the chipping 
rate or the average mass per site p. There is a critical 
line Pciw) in the p-w plane that separates two types of 
asymptotic behaviors of P(rn). For fixed w, as p is varied 
across the critical value pc{w), the large m behavior of 
P(m) was found to be, 

[ e-'"/"'* p < p,{w), 

P{m) I p = Pc{w), 

m""^ + infinite aggregate p > Pc{w). 

(286) 

As one increases p beyond pc, this asymptotic algebraic 
part of the critical distribution remains unchanged but 
in addition an infinite aggregate forms. This means that 
all the additional mass (p — Pc)V (where V is the vol- 
ume of the system) condenses onto a single site and does 
not disturb the background critical distribution. This 
is analogous, in spirit, to the condensation of a macro- 
scopic number of bosons onto the single k = mode in 
an ideal Bose gas as the temperature goes below a cer- 
tain critical value. Ref. (Rajesh and Majumdar, 2001) 
proved analytically that the mean field phase boundary, 
pc{w) = y/w + 1 — 1, is exact and independent of the 
spatial dimension d. They also provided unambiguous 
numerical evidence that the exponent t = 5/2 is also 
independent of d. The corresponding growth exponents 
are: Z = 2, a = 2/3 (Majumdar et al, 2000a). Even 
though the single site distribution P{m) may be given 
exactly by the mean field solution, that does not prove 
that mean field theory or product measure is the exact 
stationary state in all dimensions. 

The left-right asymmetric version of the CM model 
was also studied in (Majumdar et al., 2000a). This has a 
qualitatively similar phase transition in the steady state 
as the CM model but exhibits a different class phase tran- 
sition owing to the mass current in this system. Simula- 
tions in one dimension predict: Z = 1.67, a = 0.67. 

H. Unidirectionally coupled DP classes 

As it was mentioned in Section V.J in case of coupled, 
multi-species DP processes field theoretical RG analysis 
(Jansscn, 1997b) predicts DP criticality with an unsta- 
ble, symmetrical fixed point, such that sub-systems with 
unidirectionally coupled DP behavior emerge. This was 
shown to be valid for linearly coupled N-component, con- 
tact processes too. Unidirectionally coupled DP systems 



of the form (UCDP) 



A 


^ 2A 


A- 


^A + B 
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^ 2B 


B- 


^B + C 


C 


^ 2C 


C- 


^C + D 



(287) 

were investigated by (Goldschmidt et al, 1999; Tauber 
et al., 1998a) with the motivation that such models 
can describe interface growth models, where adsorption- 
desorption are allowed at terraces and edges (see 
Sect.VI.H.l). The simplest set of Langevin equations for 
such systems was set up by (Alon et al., 1998): 

dt(t)k{x,t) = a(t>k{yi,t)-X4>l{x,t)+D\/^(l)k{x,t) + 
+ p<^fe_i(x,t) + Cfe(x,i) , (288) 

where Cfc are independent multiplicative noise fields for 
level k with correlations 

(a(x,t)) = 0, (289) 
(a(x,t)0(x',O) = 2r</.fc(x,t)4,i^''(x-x')^(t-t') , 

for fc > 0, while for the lowest level (fc = 0) 0-1 = is 
fixed. The parameter a controls the offspring production, 
p is the coupling and A is the coagulation rate. As one 
can see the A; = equation is just the Langevin equation 
of DP (84). The mean- field solution of these equations 
that is valid above rfc = 4 results in critical exponents for 
the level k: 

P'm'f = 2"' • (290) 

and = 1/2 and J^^^ = 1 independently of k. For 

d < dc field theoretical RG analysis of the action for 
k < K levels 

S = Ef=o' / d'^x dt I Vfc (rdt - DV^ - a) ct>k 

-fitpkcpk-i + i V'fc (^fc - V'fc) (t>k I , (291) 

was performed in (Goldschmidt et al., 1999; Tauber 
et al., 1998a). The RG treatment run into several 
difficulties. Infrared-divergent diagrams were encoun- 
tered (Goldschmidt, 1998) and the coupling constant p 
was shown to be a relevant quantity (that means it di- 
verges under RG transformations). Ref. (Goldschmidt 
et al, 1999) argued that this is the reason why scaling 
seems to break down in simulations for large times (in 
lattice realization p, is limited). The exponents of the 
one-loop calculations for the first few levels, correspond- 
ing to the interactive fixed line as well as results of lattice 
simulations are shown in Table XXIX. These scaling ex- 
ponents can be observed for intermediate times but it is 
not clear if in the asymptotically long time they drift to 
the decoupled values or not. 
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Pi 
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0.80(4) 
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0.157(4) 
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0.75(10) 


3/4 - 0(6) 


2/Zi 


1.26(1) 


1.10(2) 


1.03(2) 




2/Z2 


1.25(3) 


1.12(3) 


1.04(2) 


l + e/24 + 0(e2) 


2/Zs 
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1.15(30) 


1.21(15) 





TABLE XXIX Critical exponents of UCDP (Goldschmidt 
et al, 1999). 



The main representatives of these classes are certain 
monomer adsorption-dcsorption models (Sect.VI.H.l) 
and polynuclear growth models (PNG) with depinning 
transitions (Kertesz and Wolf, 1989: Lchner et al, 1990; 
Toom, 1994a, b). The latter type of systems are defined 
by parallel update dynamic rules and coupled DP pro- 
cesses emerge in a co-moving frame. 

1. Monomer adsorption-desorption at terraces 

Here I show an example how a coupled particle sys- 
tem is related to an interface growth model. Refs. (Alon 
et al, 1996, 1998) defined SOS and RSOS models that 
can be mapped onto UCDP (Sect. VI. H). In this mod- 
els adsorption an desorption processes may take place at 
terraces and edges. For each update a site i is chosen at 
random and an atom is adsorbed 

hi ^ hi + 1 with probability q (292) 

or desorbed at the edge of a plateau 

hi — > min(/ij, /ij+i) with probability (1 — q)/2 , 
hi min(/ii, /ij-i) with probability (1 — q)/2 . 

(293) 

Identifying empty sites at a given layer as A particles, 
the adsorption process can be interpreted as the decay 
of A particles {A—* 0), while the desorption process cor- 
responds to A particle production {A — > 2 A). These 
processes generate reactions on subsequent layers, hence 
they are coupled. The simulations in Id have shown 
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that this coupling is relevant in the upward direction 
only hence the model is equivalent to the UCDP pro- 
cess. Defining the order parameters on the k — th layer 
as 

1 

where hi is the height at site i, they are expected to scale 
as 

rik^iqc-qf"" ■ fc = 1,2,3,... (295) 

By varying the growth rate (q) these models exhibit a 
roughening transition at qc = 0.189 (for RSOS) and at 
qc = 0.233(1) (for SOS) from a non-moving, smooth 
phase to a moving, rough phase in one spatial dimen- 
sion. The /3'^ (and other exponents) take those of the Id 
UCDP class (see Table XXIX). 

The scaling behavior of the interface width is charac- 
terized by many different length scales. At criticality it 
increases as 

W'^{t)(XThit (296) 

with T ~ 0.102(3) for RSOS (Hinrichsen and Mukamel, 
2003). 

I. Unidirectionally coupled PC classes 

Similarly to the UCDP case (see Sect. VI. H) surface 
growth processes of dimers stimulated the introduction 
of unidirectionally coupled BARW2 (Sect.IV.D.l) models 
(Hinrichsen and Odor, 1999a,b): 

A^3A B^3B C^3C 

2A^0 2B^0 2C^0 (297) 

A^A + B B^B + C C^C + D,... 

generalizing the concept of UCDP. The mean field ap- 
proximation of the reaction scheme (297) looks as 

dfUA = (TUA - Xn\ , 

dtUB = (TUB — ^^n^ + l^nA , (298) 
dtuc = cmc - Xric + tJ^ns , • • • 

where nA,nB,nc correspond to the densities no,ni,n2 
in the growth models, a and A are the rates for offspring 
production and pair annihilation respectively. The coeffi- 
cient fj, is an effective coupling constant between different 
particle species. Since these equations are coupled in only 
one direction, they can be solved by iteration. Obviously, 
the mean-field critical point is ctc = 0. For small values 
of a the stationary particle densities in the active state 
are given by 



corresponding to the mean field critical exponents 

= /3^^^ = l/2, P'J^ ^l/A,.... (300) 

These exponents should be valid for d > dc ^ 2. Solv- 
ing the asymptotic temporal behavior one finds = 1, 
implying that (5f ^ = 2"'=. 

The effective action of unidirectionally coupled 
BARW2's should be given by 

S[lpo,tpl,tp2, ■ ■ . ,l/io,V'l,V'2, ■ • ■] = 

/oo 
d'^xdt Y,{Mdt - DV^)^k - A(l - ^1)^1 + 

k=0 

+(7(1 - ipl)ipk'il'k + m(i - ■4'k)i>k-i'4>k^i I (301) 

where V'-i = "0-1 = 0- Here the fields and ipu repre- 
sent the configurations of the system at level k. Since 
even the RG analysis of the one component BARW2 
model suffered serious problems (Cardy and Tauber, 
1996) the solution of the theory of (301) seems to be 
hopeless. Furthermore one expects similar IR diagram 
problems and diverging coupling strengths as in case of 
UCDP that might be responsible for violations of scaling 
in the long time limit. 

Simulations of a 3-component model in Id coupled by 
instantaneous particle production of the form (297) re- 
sulted in decay exponents for the order parameter defined 
as (294): 

5a = 0.280(5) , 5b = 0.190(7) , 5c = 0.120(10) , (302) 

For further critical exponents see Sect. VI. 1. 1. 

It would be interesting to investigate parity-conserving 
growth processes in higher dimensions. Since the up- 
per critical dimension d'^ is less than 2, one expects the 
roughening transition ~ if still existing - to be described 
by mean-field exponents. In higher dimensions, n-mers 
might appear in different shapes and orientations. 

The main representatives of these classes are certain 
dimer adsorption-desorption models (Sect.VI.H.l) and 
polynuclear growth models (PNG) with depinning tran- 
sitions (Hinrichsen and Odor, 1999a). The latter type 
of systems are defined by parallel update dynamic rules 
and coupled DP processes emerge in a co-moving frame. 

1. Dimer adsorption-desorption at terraces 

Similarly to the monomer case (Sect.VI.H.l) dimer 
adsorption-desorption models were defined (Hinrich- 
sen and Odor, 1999a,b; Noh et ai, 2000). With the 
restriction that desorption may only take place at the 
edges of a plateau the models can be mapped onto 
the unidirectionally coupled BARW2 (Sect.VI.I). The 
dynamical rules in d = 1 are defined in Fig. 27. The 
mapping onto unidirectionally coupled BARW2 can be 
seen on Fig. 28. 
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FIG. 27 Dimers are adsorbed with probability p and desorbed 
at the edges of terraces with probability 1 — p. Evaporation 
from the middle of plateaus is not allowed (Hinrichsen and 
Odor, 1999a). 



ccccccccccc 




FIG. 28 Extended particle interpretation. Dimers are ad- 
sorbed {2A 0) and desorbed {A — > 3A) at the bottom 
layer. Similar processes take place at higher levels (Hinrich- 
sen and Odor, 1999a). 



Such dimer models can be defined in arbitrary spatial 
dimensions. In (Hinrichsen and Odor, 1999a) four one 
dimensional variants were investigated: 

1) Variant A is a restricted solid-on-solid (RSOS) model 
evolving by random sequential updates. 

2) Variant B is a solid-on-solid (SOS) model evolving by 
random sequential updates. 

3) Variant C is a restricted solid-on-solid (RSOS) model 
evolving by parallel updates. 

4) Variant D is a solid-on-solid (SOS) model evolving by 
parallel updates. 

Variants A and B exhibit transitions in contrast with 
PNG models, in which only parallel update rules permit 
roughening transitions. By varying the adsorption rate 
p the phase diagram shown in Fig. 29 emerges for RSOS 
and SOS cases. If p is very small, only a few dimers are 
adsorbed at the surface, staying there for a short time 
before they evaporate back into the gas phase. Thus, the 
interface is anchored to the actual bottom layer and does 
not propagate. In this smooth phase the interface with 
growths logarithmically until it saturates to a finite value 
(even for L — s- cx)). 

As p increases, a growing number of dimers covers the 
surface and large islands of several layers stacked on top 
of each other are formed. Approaching a certain critical 
threshold pc the mean size of the islands diverges and the 
interface evolves into a rough state with the the finite-size 
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; roughening! 
transition 



faceting 
transition 



1/2 



Smooth 
RoUo phase 



W = a lii[t F(te^ii)] 



SOS 



Rough phase 

algebraic growth 
crosses over to : 
W ~ In (t) 



Faceted 
phase: 
W - In (t) 




FIG. 29 Phase diagram of Id dimer models. 



variant 


A 


B 


C 


D 


restriction 


yes 


no 


yes 


no 


updates 


random 


random 


parallel 


parallel 


Pc 


0.3167(2) 


0.292(1) 


0.3407(1) 


0.302(1) 


a 


0.172(5) 


0.23(1) 


0.162(4) 


0.19(1) 


Z 


1.75(5) 


1.75(5) 


1.74(3) 


1.77(5) 


So 


0.28(2) 


0.29(2) 


0.275(10) 


0.29(2) 


5i 


0.22(2) 


0.21(2) 


0.205(15) 


0.21(2) 
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0.14(2) 


0.14(3) 


0.13(2) 


0.14(2) 


a 


1.2(1) 


undefined 


1.25(5) 


undefined 


P 


0.34(1) 


0.50(1) 


0.330(5) 


0.49(1) 



TABLE XXX Numerical estimates for the four variants of 
the dimer model at the roughening transition p = pc (upper 
part) and at the transition p = 0.5 (lower part). 



scaling form 



W^{L,t) ~ a In tG{t/L^) 



(303) 



The order parameter defined on the fc-th layer as (294) 
exhibits unidirectionally coupled BARWe critical behav- 
ior. The transition rates and exponents are summarized 
in Table XXX. Above pc one may expect the interface 
to detach from the bottom layer in the same way as the 
interface of monomer models starts to propagate in the 
supercritical phase. However, since dimers are adsorbed 
at neighboring lattice sites, solitary unoccupied sites may 
emerge. These pinning centers prevent the interface from 
moving and lead to the formation of 'droplets'. Due to 
interface fluctuations, the pinning centers can slowly dif- 
fuse to the left and to the right. When two of them meet 
at the same place, they annihilate and a larger droplet 
is formed. Thus, although the interface remains pinned, 
its roughness increases continuously. The width initially 
increases algebraically until it slowly crosses over to a 
logarithmic increase 

The restricted as well as the unrestricted variants un- 
dergo a second phase transition at p = 0.5 (Noh et at, 



2000) where the width increases algebraically with time as 
Wr^t'^. In the RSOS case ordinary Family- Vicsek scal- 
ing (257 occurs with the exponents given in Table XXX. 
The dynamic exponent is Z — a/P 3. This value 
stays the same if one allows dimer digging at the faceting 
transition but other surface exponents a ~ 0.29(4) and 
/3 ~ 0.111(2) will be different (Noh et al, 2000). An ex- 
planation of the latter exponents is given in (Noh et al., 

2001) based on mapping to globally constrained random 
walks. In the SOS cases (variants B,D) large spikes are 
formed, the surface roughens much faster with a growth 
exponent of ~ 0.5. The interface evolves into con- 
figurations with large columns of dimers separated by 
pinning centers. These spikes can grow or shrink almost 
independently. As the columns are spatially decoupled, 
the width does not saturate in finite systems, i.e., the 
dynamic exponents a and Z have no physical meaning. 

For p > 0.5 the restricted models A and C evolve into 
faceted configurations. The width first increases alge- 
braically until the pinning centers become relevant and 
the system crosses over to a logarithmic increase of the 
width. Therefore, the faceted phase may be considered 
as a rough phase. The unrestricted models B and D, 
however, evolve into spiky interface configurations. The 
spikes are separated and grow independently by depo- 
sition of dimers. Therefore, W'^ increases linearly with 
time, defining the free phase of the unrestricted models. 

In the simulations mentioned by now the interface was 
grown from flat initial conditions. It turns out that start- 
ing with random initial conditions /i^ = 0, 1 the den- 
sities rife turn out to decay much slower. For restricted 
variants an algebraic decay of rip with an exponent 



6o ~ 0.13 



(304) 



was observed. Similarly, the critical properties of the 
faceting transition at p = 0.5 are affected by random 
initial conditions. The non-universal behavior for ran- 
dom initial conditions is related to an additional parity 
conservation law. The dynamic rules not only conserve 
parity of the particle number but also conserve the par- 
ity of droplet sizes. Starting with a flat interface the 
lateral size of droplets is always even, allowing them to 
evaporate entirely. However, for a random initial config- 
uration, droplets of odd size may be formed which have 
to recombine in pairs before they can evaporate, slow- 
ing down the dynamics of the system. In the language 
of BARW2 processes the additional parity conservation 
law is due to the absence of nearest-neighbor diffusion. 
Particles can only move by a combination of offspring 
production and annihilation, i.e., by steps of two lattice 
sites. Therefore, particles at even and odd lattice sites 
have to be distinguished. Only particles of different par- 
ity can annihilate. Starting with a fully occupied lattice 
all particles have alternating parity throughout the whole 
temporal evolution, leading to the usual critical behavior 
at the PC transition. For random initial conditions, how- 
ever, particles of equal parity cannot annihilate, slowing 
down the decay of the particle density. Similar sector 
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CLASS ID 


features 


Section 


DP 


time reversal symmetry 


IV.A 


DyP 


long memory 


IV.B 


VM 


Zi, symmetry 


IV.C 


PCP 


coupled frozen field 


V.E 


NDCF 


global conservation 


V.I 


PC 


Z2 symmetry, BARW2 conservation 


IV.D.l 


BP 


DP coupled to ARW 


V.F 


DCF 


coupled diffusive, conserved field 


V.H 


N-BARW2 


N-component BARW2 conservation 


V.K 


N-BARW2S 


symmetric NBARW2 + exclusion 


V.L 


N-BARW2a 


asymmetric NBARW2 + exclusion 


V.L 



TABLE XXXI Summary of known absorbing state universal- 
ity classes in homogeneous isotropic systems. 



decomposition has been observed in diffusion of A;-nier 
models (Barma and Dhar, 1994; Barma et al., 1993). 



spatially isotropic systems. Those which arc below the 
horizontal line exhibit fluctuating absorbing states. The 
necessary and sufRcient conditions for these classes are 
usually unknown. The mean-field classes can also give 
a guide to distinguish classes below dc- In Table XXXII 
I collected the mean-field exponents and upper critical 



CLASS 


/3 




z 




a 


S 


V 


dc 


DP 


1 


1 


2 


1 


1 


1 





4 


DyP 


I 


I 


2 


I 


1 


1 





6 


VM 





I 


2 


I 





1 





2 


PC 


I 





2 


I 


I 





-1/2 


2 


BkARW 


1 

fc-i 





2 


1 


1 

fc-i 








2 


PARWs 


1 





2 


n 


I/n 










PARWa 


1 

m — n 




2 


m— 1 
m — n 


1 

m — 1 








NDCF 


1 


I 


2 


I 


I 


I 





4 


NBARW2 


1 





2 


1 


1 








2 



TABLE XXXII Mean-field classes of known, homogeneous 
absorbing-state transitions 



VII. SUMMARY 



In summary dynamical extensions of classical equilib- 
rium classes were introduced in the first part of this re- 
view. New exponents, concepts, sub-classes, mixing dy- 
namics and some unresolved problems were discussed. 
The common behavior of these models was the strongly 
fluctuating ordered state. In the second part genuine 
nonequilibrium dynamical classes of reaction-diffusion 
systems and interface growth models were overviewed. 
These were related to phase transitions to absorbing 
states of weakly fluctuating ordered states. The class be- 
havior is usually determined by the spatial dimension, 
symmetries, boundary conditions and inhomogeneities 
like in case of equilibrium models but in low dimensions 
hard-core exclusion was found to be a relevant factor too, 
splitting up criticality in fermionic and bosonic models. 
The symmetries are not so evident as in case of equi- 
librium models, they are expressed in terms of the rela- 
tions of fields and response fields most precisely. Further- 
more in recently discovered coupled systems with binary, 
triplet or quadruplet particle production no special sym- 
metry seems to be responsible for a novel type of criti- 
cal behavior. Perhaps a proper field theoretical analysis 
of the coupled DP-I-ARW system could shed some light 
on this mystery. The parity conservation in hard-core 
and in binary spreading models seems to be irrelevant. 
In Table XXXI I summarized the most well known fami- 
lies of absorbing phase transition classes of homogeneous, 



Note that in the general nA — > {n + k)A, mA — > {m—l)A 
type of RD systems the values of n and m determine the 
critical behavior. 

In d > 1 dimension the mapping of spin-systems onto 
RD systems of particles is not straightforward, instead 
one should also deal with the theory of branching inter- 
faces (Cardy, 1998). Preliminary simulations found in- 
teresting critical phenomena by generalized Potts models, 
exhibiting absorbing states (Lipowski and Droz, 2002a). 

Further research should also explore the universality 
classes of such nonequilibrium phase transitions that oc- 
cur by external current driven systems or by other mod- 
els exhibiting fluctuating ordered states (Evans, 2000; 
Evans et ai, 1998). Nonequilibrium phase transitions 
in quantum systems (Racz, 2002) or by irregular graph 
or network based systems are also of current interest of 
research. Finally, having settled the problems arisen by 
fundamental nonequilibrium models one should turn to- 
wards the study of more application motivated systems. 
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Abbreviations 



ABC 


active boundary condition 


(iV.A.4) 


AKW 


annihilating random walk (^AA — > 0) 


(1 V .O.i ) 


Ar 


annihilation-fission process 


(V-i* j 


rSArtW 


branching and annihilation random walk 


(IV.A) 


rSAxlWC 


even-offspring branching and annihilation random walk 


(1 V .U.i j 


rsAKWo 


odd-ofFspring branching and annihilation random walk 


(iV.A.o) 


rsAKWz 


two-offspring branching and annihilation random walk 


(iV.U.i) 


rSKArtW 


branching and k particle annihilating random walk 


(IV .rj ) 


Or 


binary production 


(V.I') 


<^AiVi 


coherent anomaly method 


(iV.A) 


KjDr 


compact directed percolation 


(IV.A./) 




diffusive conserved field 


(V.H) 


Ulv 


Domany-Kinzel cellular automaton 


CWT A 0^ 
(1 V .A.Z j 


L)ir 


directed percolation 


(IV.A) 


Do 


damage spreading 


(i.D.i) 


UyF 


dynamical percolation 


(IV.B) 




Jildwards- W ilkmson 


(VLB) 




density matrix renormalization group 


(IV) 




generalized epidemic process 


(IV.B) 




generalized Domany-Kinzel cellular automaton 


(iV.JJ.oJ 


/-^ TV /TTTi 


generalized mean-field approximation 


/'TT^ 

(II) 




inactive boundary condition 


[iM .A. 4 ) 


xVJr Zj 


Kardar-Parisi- Zhang 


(Vl.Uj 


T TAyf 
ijllVi 


linear interface model 


(VI.C) 


IVlr 


mean-field approximation 


(i.t) 


AT 1I> A O^^/^O 


even-offspring, N-component branching and annihilation random walk 


(V.K) 




non-Qittusive conserved nelu 


(V.I) 


ATTTT^TlV/r 
IN ILlvliVl 


nonequilibrium Ising model 


(^i V .JJ.z ) 




parity conserving 


(IV.D) 




pair contact process 


(V.E) 




pair contact process with particle diffusion 


(V.F.I) 


PAKWs 


symmetric production and m-particle annihilating random walk 


(IV.F) 


PAKWa 


asymmetric production and m-particle annihilating random walk 


(IV.F) 




polimiclear growth models 


(VI.H) 




reflecting boundary condition 


('W7 A A\ 

(iV .A.4 j 


rlJJ 


reaction-diffusion 


(I.E) 




renormalization group 


(I.E) 


RSOS 


restricted solid on solid model 


(VI) 


SCA 


stochastic cellular automaton 


(IV.A.2) 


SOS 


solid on solid model 


(VI) 


TTP 


threshold transfer process 


(V.E) 


UCDP 


unidirectionally coupled directed percolation 


(VI.H) 


VM 


voter model 


(IV.C) 
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